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1  Statement  of  the  problem 


Dramatic  innovations  in  mission  critical  electronic  equipment  have  led  to  an  increased  demand 
for  portable  power  in  the  field.  In  addition,  advances  in  small  mobile  autonomous  robotic 
devices,  e.g.  unmanned  aerial  vehicles  (UAV’s)  or  unmanned  ground  vehicles  (UGV’s),  as  well 
as  biomechatronic  exoskeletons  [Pons,  2008]  are  limited  by  power  demands  that  go  beyond  the 
capabilities  of  current  battery  technologies.  Replacing  batteries  by  fuel  cells  is  a  promising 
approach,  as  fuel  cells  offer  higher  power  densities,  do  not  include  moving  parts,  and  thus  are 
free  of  noise  and  vibrations.  If  deployed,  fuel  cell  based  power  systems  can  greatly  reduce  the 
weight  burden  compared  to  conventional  batteries,  ease  refueling,  and  increase  mission  duration. 

At  the  current  state  of  development,  fuel  cell  systems  face  the  critical  issue  of  hydrogen 
generation  and  storage  [NRC,  2002].  A  potential  solution  to  this  problem  is  fuel  reforming, 
where  a  hydrogen -rich  syngas  is  generated  from  liquid  hydrocarbon  fuels.  Liquid  hydrocarbons 
have  the  advantage  of  high  volumetric  energy  densities  and  an  existing  infrastructure  with  the 
potential  of  using  standard  battlefield  fuels.  Those  fuels,  however,  are  not  tailored  to  the  specific 
needs  for  fuel  reforming  and  are  further  hampered  by  a  high  sulfur  content  of  up  to  0.3  wt% 
[MIL-DTL-83133F,  2008].  Fuel  reforming  based  on  conventional  catalytic  techniques  is 
difficult,  as  operating  conditions  and  fuel  quality  need  to  meet  stringent  requirements  to  prevent 
catalyst  degradation  due  to  excessive  temperatures  and  sulfur  poisoning  [Moon  et  al.,  2004].  An 
alternative  approach  is  non-catalytic  fuel  reforming,  which  promises  a  more  robust  design  due  to 
the  lack  of  damageable  catalytic  surfaces. 

Non-catalytic  fuel  reformers  have  the  advantage  of  not  containing  damageable  catalytic  surfaces 
and  thus  can  operate  at  a  wide  range  of  temperatures  with  no  restrictions  regarding  quality  and 
sulfur  content  of  the  fuel.  Optimal  operating  conditions  for  non-catalytic  fuel  reforming  are  in 
the  very  fuel-rich  regime,  either  near  or  beyond  the  conventional  rich  flammability  limit,  where 
reaction  rates  are  slow  and  some  means  of  enhancing  the  rates  is  necessary  for  practical  reactors. 
One  option  for  enhancing  the  reaction  rates  is  to  recirculate  heat  such  that  local  temperatures  are 
superadiabatic  i.e.  higher  than  those  predicted  by  equilibrium.  The  classic  means  of  producing 
superadiabatic  temperatures  is  by  reaction  of  a  fuel/air  mixture  in  a  filtration  reactor  consisting 
of  a  single  flow  channel  filled  with  a  porous  matrix  [Dhamrat  and  Ellzey,  2006;  Fay,  et  al,  2005, 
Bingue,  et  al.  2002](Fig.  1.1a).  Filtration  combustion  is  characterized  by  a  propagating  reaction 
front  that  is  self-sustaining  once  it  is  established.  In  the  ultra-rich  regime  typical  for  fuel 
reforming,  the  reaction  front  propagates  downstream  (co-flowing)  into  a  region  where  the  solid 
has  already  been  heated  by  the  hot  product  gases,  thus  producing  superadiabatic  temperatures 
locally.  The  high  local  temperatures  are  sufficient  to  drive  the  reactions  of  mixtures  that  are 
otherwise  not  flammable.  One  of  the  major  practical  limitations,  however,  of  the  filtration  reactor 
is  that  it  cannot  be  operated  continuously.  Due  to  the  downstream  propagation,  the  reaction  front 
ultimately  reaches  the  end  of  the  bed  and  the  process  must  be  restarted  or  the  flow  direction 
reversed,  thus  complicating  the  design  of  the  reactor. 
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Outlet 
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channel 


Post-reaction  zone: 
heat  transfer 
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(a)  Filtration  reactor  (b)  Counterflow  reactor 

Figure  1.1:  Reactors  used  in  this  study 

In  order  to  overcome  the  limited  operation  time  of  a  filtration  reactor,  we  have  developed  a  novel 
counterflow  reactor  concept  with  stationary  combustion  zones  [Schoegl  and  Ellzey,  2007; 
Schoegl  and  Ellzey  2009]that  is  capable  of  continuous  operation.  In  this  design,  heat  is 
transferred  from  hot  reaction  products  to  cold  reactants  across  a  wall  dividing  adjacent  flow 
channels  (Fig.  1.1b).  The  operating  principle  is  based  on  preheating  of  the  reactants  by 
counterflow  heat  exchange,  which  enables  combustion  of  otherwise  non-flammable  mixtures 
typical  for  fuel  reforming.  Thus,  superadiabatic  temperatures  are  attained  with  stationary  reaction 
zones  allowing  continuous  operation. 

In  the  research  conducted  with  ARO  support,  both  reactor  designs  were  used  to  study  the 
feasibility  of  non-catalytic  reforming  of  butanol,  heptane  and  Jet-A.  The  more  rugged  filtration 
reactor  (Fig.  1.1a)  was  used  to  demonstrate  the  conversion  of  butanol  and  Jet-A,  and  the 
counterflow  reactor  (Fig.  1.1b)  was  used  for  heptane.  Together  these  results  greatly  increase  our 
understanding  of  the  potential  for  non-catalytic  reforming. 


1.1  Project  Objectives 


To  address  the  challenges  regarding  the  increasing  power  demands  in  portable  or  mobile  devices, 
we  investigated  the  principles  of  a  compact,  non- catalytic  fuel  reformer  to  produce  a  hydrogen- 
rich  syngas  from  heptane,  butanol,  and  Jet-A.  The  purpose  of  this  project  was  to  (1)  demonstrate 
the  principle  of  non-catalytic  reforming  of  butanol  and  Jet-A,  (2)  demonstrate  that  liquid  fuels 
can  be  reformed  in  the  counterflow  reactor,  and  (3)  understand  the  critical  design  parameters 
controlling  the  conversion  efficiency  in  the  counterflow  reactor  and  develop  scaling  laws  for 
down-sizing  the  counterflow  reactor  for  applications  to  systems  with  delivered  power  of  200- 
1000  W. 
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In  Section  2,  a  summary  of  the  major  findings  is  presented.  Sub-sections  2. 1-2.3  discuss  results 
related  to  the  specific  objectives  of  the  project.  Sub-sections  2.4  and  2.5  present  additional 
results  related  to  the  project. 

2  Summary  of  most  important  results 

In  this  section,  the  most  important  results  from  the  study  of  non-catalytic  conversion  of  liquid 
and  gaseous  fuels  are  presented.  More  detailed  results  are  included  in  the  appendixes. 

2.1  Non-catalytic  reforming  of  butanol  and  Jet-A 

Both  butanol  and  Jet-A  were  successfully  converted  in  the  filtration  reactor.  Results  of  the 
conversion  of  butanol  and  jet  fuel  to  syngas  via  non-catalytic  transient  filtration  combustion  are 
summarized:  At  maximum,  about  42%  of  the  hydrogen  in  jet  fuel  was  converted  to  H2  and  56% 
of  carbon  was  converted  to  CO.  The  H2  yield  continued  to  increase  with  (J)  in  the  jet  fuel 
experiments  whereas  in  the  butanol  experiments  the  yields  for  H2  and  CO  both  reached  peaks 
within  the  tested  operating  range.  The  peak  CO  yield  for  experiments  with  butanol  was  72%  and 
the  peak  H2  yield  was  43%.  Most  of  the  chemical  energy  exiting  the  reactor  was  bound  in  H2 
and  CO,  but  CH4,  C2H2  and  C2H4  were  observed  in  substantial  amounts  in  experiments  with 
both  fuels,  especially  in  experiments  with  butanol  at  high  (j).  In  contrast,  equilibrium  predicted 
nearly  negligible  amounts  of  energy-containing  compounds  besides  H2,  CO  and  solid  carbon. 

Soot  production  proved  to  be  substantial  in  experiments  with  both  fuels,  with  the  qualitative 
observation  of  soot  from  both  fuels  exceeding  the  soot  production  observed  in  similar 
experiments  with  methane,  heptane,  and  ethanol.  Soot  production  limited  the  maximum  testable 
(j)  for  jet  fuel  experiments  to  <j)  =  3.15  since  soot  clogged  the  reactor  when  (j)  exceeded  this  value. 
An  analysis  of  the  carbon  balance  for  jet  fuel  experiments  showed  that  up  to  40%  of  the  carbon 
entering  the  system  was  deposited  on  the  porous  medium.  Similarly  for  butanol,  the  carbon 
balance  showed  significant  loss  up  to  (j)  ~  3.  At  the  most  extreme  (|>,  the  reaction  zone  propagates 
downstream  oxidizing  some  of  the  soot  that  has  been  deposited  previously,  resulting  in  an 
exhaust  stream  with  more  carbon  atoms  than  are  contained  in  the  inlet  fuel  at  the  time  of  sample. 

Equilibrium  calculations  and  experimental  results  were  analyzed  to  determine  trends  in  fuel 
conversion  efficiency  as  a  function  of  fuel  type.  The  equilibrium  calculations  predicted  that 
larger  fuel  molecules  would,  generally,  produce  higher  syngas  yields,  but  the  experimental 
results  disagreed  with  this  trend.  Generally,  experimental  results  show  that  smaller  hydrocarbons 
produce  higher  syngas  yields  with  jet  fuel  having  the  lowest  hydrogen  and  carbon  monoxide 
yields.  The  comparison  of  these  equilibrium  data  with  the  experimental  data  available  in  the 
literature  suggests  that  kinetic  effects  are  very  important  in  fuel  conversion  to  syngas  and 
equilibrium  results  should  be  interpreted  carefully  when  predicting  syngas  production  from  the 
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combustion  of  very  rich  mixtures  of  fuel  and  air.  Lastly,  the  demonstration  of  the  conversion  of 
jet  fuel  and  butanol,  along  with  previous  demonstrations  with  methane,  ethanol  and  heptane 
show  the  robustness  of  the  non-catalytic  filtration  reactor. 

Further  details  on  the  non-catalytic  conversion  of  Jet-A  and  butanol  are  presented  in  Appendix 
A. 

2.2  Reforming  of  heptane  in  counterflow  reactor 

Experiments  demonstrated  the  non-catalytic  reforming  of  heptane  to  produce  hydrogen-rich 
syngas  in  a  counterflow  reactor.  A  key  advantage  of  this  reactor  is  that  the  reaction  zones  are 
stationary,  which  permits  continuous  operation  and  practical  applications.  Stable  operation  of 
the  reactor  with  premixed  heptane/air  reactants  was  established  for  equivalence  ratios  between 
2.8  and  3.9  and  inlet  velocities  between  50  cm/s  and  200  cm/s.  Peak  wall  temperatures  were 
monitored  throughout  reactor  operation  and  were  at  all  points  in  excess  of  the  adiabatic 
equilibrium  temperature,  thus  excess-enthalpy  flames  were  achieved  through  internal  heat 
recirculation. 

Exhaust  gas  composition  was  measured  using  gas  chromatography  for  two  sets  of  tests: 
equivalence  ratio  was  varied  from  2. 9-3. 8  with  inlet  velocity  held  constant  at  125  cm/s,  and 
equivalence  ratio  was  held  constant  at  3.0  with  inlet  velocity  varied  from  50-200  cm/s.  H2 
concentration  decreased  with  increasing  equivalence  ratio  from  a  maximum  of  14.5%  to  a 
minimum  of  10.6%,  while  CO  concentration  remained  nearly  constant  around  17.0%.  Both  H2 
and  CO  concentrations  showed  little  dependence  overall  on  inlet  velocity,  with  the  former 
varying  from  12.8-14.4%  and  the  latter  varying  between  16.6-17.5%.  The  largest  effect  of  inlet 
velocity  was  found  at  the  lower  end  of  the  tested  range  where  rates  of  heat  release  and 
temperatures  are  lower. 

In  all  cases,  unburned  hydrocarbons  were  measured  in  excess  of  concentrations  predicted  by 
equilibrium  composition.  Total  unbumed  hydrocarbon  levels  increased  with  increasing 
equivalence  ratio.  Soot  was  noted  with  significant  formation  observed  at  higher  equivalence 
ratios  and  inlet  velocities.  These  results  and  observations  strongly  suggested  that  the  extent  of 
fuel  breakdown  is  reduced  by  insufficient  reaction  rate  and  residence  time. 

The  results  of  this  study  showed  trends  similar  to  previously  published  studies  on  syngas 
production  from  methane  and  propane  using  the  counterflow  reactor  [Schoegl  and  Ellzey,  2007 ; 
Schoegl  and  Ellzey  2009].  Each  fuel  had  H2  species  and  energy  conversion  efficiencies  below 
equilibrium  values,  but  total  conversion  efficiencies  were  comparable  to  equilibrium.  In  addition, 
energy-dense  unburned  hydrocarbon  concentrations  were  in  excess  of  equilibrium  with 
concentrations  an  order  of  magnitude  higher  than  predicted  by  equilibrium  at  the  highest  tested 
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equivalence  ratios  in  this  study.  The  comparison  between  the  different  fuels  indicates  that, 
although  the  stable  operating  range  for  heptane  does  include  the  theoretically  optimal  value  as 
predicted  by  hydrogen  conversion  efficiencies  calculated  from  equilibrium,  the  current  reactor 
design  does  not  permit  complete  fuel  conversion  at  high  equivalence  ratios.  In  addition,  the 
formation  of  soot  at  some  conditions  suggests  the  importance  of  operating  at  more  moderate 
conditions  for  optimal  performance. 

This  study  demonstrated  the  fuel  flexibility  of  the  counterflow  reactor  by  reforming  a  liquid  fuel, 
and  results  supplemented  previous  reforming  studies  on  the  counterflow  reactor.  Peak  H2 
conversion  efficiency  was  significantly  lower  than  that  attained  via  filtration  combustion  in 
porous  media,  an  alternate  method  of  non-catalytic  reforming.  Differences  were  attributed  to 
higher  peak  temperatures  in  the  porous  media  reactor,  which  has  higher  specific  surface  area  for 
heat  transfer  as  compared  to  the  counterflow  reactor  under  investigation.  Furthermore,  although 
the  tested  equivalence  ratio  range  in  this  study  does  include  (|>  =  2.9,  where  peak  conversion 
efficiency  is  predicted  to  occur  by  equilibrium,  previous  studies  predict  peak  conversion 
efficiencies  over  a  range  of  4>  =  1-8  to  cj)  =  3.5.  The  results  of  hydrogen  conversion  in  this  study 
indicate  that  the  optimum  equivalence  ratio  might  be  outside  the  operating  range  of  the 
counterflow  reactor. 

Further  details  of  non-catalytic  reforming  of  heptane  in  the  counterflow  reactor  are  presented  in 
Appendix  B. 

2.3  Critical  design  parameters  for  counterflow  reactor 

A  computational  study  was  undertaken  to  determine  the  effect  of  geometric  scale  on  fuel 
reforming  in  the  meso-scale  parallel-plate  counterflow  reactor.  Length  and  height  of  the  reactor 
channels  were  scaled  relative  to  the  original  dimensions  of  the  counterflow  reactor  that  has  been 
used  in  previous  analytical  and  experimental  investigations.  Results  were  quantified  in  terms  of 
steady  state  reactor  operation,  heat  recirculation  efficiency,  temperature  ratio,  and  hydrogen 
conversion  efficiency. 

Reactor  operating  range  was  found  to  be  strongly  impacted  by  channel  scaling.  Decreasing 
channel  height  resulted  in  a  shift  in  stable  operating  range  towards  higher  equivalence  ratios. 
Decreasing  channel  length  resulted  in  a  significant  narrowing  of  the  stable  equivalence  ratio  and 
inlet  velocity  ranges.  Most  stable  operating  points  were  above  the  conventional  rich  flammability 
limit  and  adiabatic  flame  speed  of  methane  and  air  combustion  based  upon  inlet  conditions, 
which  highlights  the  importance  of  heat  recirculation  in  achieving  operation  of  the  counterflow 
reformer  at  these  conditions. 
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Heat  recirculation  efficiency  and  temperature  ratios  were  found  to  be  significantly  higher  for 
reactors  with  decreased  channel  height.  These  findings  indicate  more  effective  heat  recirculation 
with  smaller  channel  height,  where  channel  wall  surface  area-to-gas  volume  ratio  is  increased 
and  the  transverse  distance  for  mass  and  heat  diffusion  is  reduced.  Channel  length  showed  a  less 
significant  impact  on  heat  recirculation,  with  shorter  channel  lengths  producing  slightly  higher 
heat  recirculation  efficiencies  and  temperature  ratios  than  longer  lengths. 

Hydrogen  conversion  efficiency  decreased  with  increasing  equivalence  ratio  for  all  reactors,  with 
slightly  higher  conversion  efficiencies  found  for  smaller  channel  heights  and  lengths.  Hydrogen 
conversion  efficiency  initially  increased  with  inlet  velocity  at  the  lowest  tested  firing  rates,  and 
decreased  at  the  highest  tested  inlet  velocities,  suggesting  that  an  optimal  firing  rate  exists. 

Scaled  down  reactors  were  achieved  by  decreasing  reactor  channel  height  and  length,  although 
the  latter  resulted  in  a  significantly  narrowed  stable  range.  The  results  of  this  study  highlight  the 
effectiveness  of  small  scale  reactor  channels  in  achieving  heat  recirculation  and  significantly 
superadiabatic  temperatures.  These,  in  turn,  permit  the  effective  extension  of  flammability  limits 
to  stable  operation  at  rich  equivalence  ratios.  While  increasing  heat  recirculation  may  be  required 
to  significantly  extend  the  stable  operating  range,  the  results  of  this  study  suggest  that 
superadiabatic  temperatures  serve  to  extend  the  flammable  range  of  reactor  operating  points  but 
otherwise  do  not  strongly  impact  hydrogen  conversion  efficiencies. 

Further  details  of  the  scaling  study  for  the  counterflow  reactor  are  presented  in  Appendix  C. 

2.4  Operation  of  the  counterflow  reactor  at  lean  conditions 

The  funding  from  ARO  was  leveraged  with  funding  from  additional  sources  which  allowed  the 
investigation  of  related  questions  on  the  operation  of  the  counterflow  reactor.  Specifically,  lean 
premixed  combustion  of  methane,  propane  and  heptane  was  studied  in  the  counterflow  reactor. 
This  reactor  design  was  previously  used  to  study  the  conversion  of  rich  methane,  propane  and 
heptane  reactant  mixtures  to  hydrogen -rich  synthesis  gas.  The  current  study  expanded  the 
established  capability  of  the  counterflow  reactor  to  operate  in  the  lean  combustion  regime  on 
increasingly  complex  fuels.  In  particular,  liquid  fuels  are  generally  easier  and  safer  to  transport 
than  gaseous  fuels,  and  are  studied  because  of  their  logistical  importance. 

Results  included  reactor  stability  maps  for  methane,  propane  and  heptane  combustion  that  show 
stable  operating  points,  as  well  as  points  where  operation  became  unstable  due  to  flashback, 
blow-off  or  extinction,  for  rich  and  lean  reactant  mixtures.  Emissions  of  CO,  NOx  and  UHC 
were  presented  for  lean  methane,  propane  and  heptane  operation.  Additionally,  operating  regions 
were  defined  in  terms  of  peak  reactor  wall  temperatures.  Similarities  in  trends  between  the  three 
tested  fuels  highlight  the  correlation  of  emissions  with  peak  reactor  wall  temperature,  and 
support  previous  findings  that  a  minimum  firing  rate  is  necessary  for  optimal  performance. 
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Further  analysis  highlights  the  importance  of  achieving  peak  wall  temperatures  near,  or  in  excess 
of,  the  calculated  adiabatic  equilibrium  temperature  of  the  fuel  at  inlet  conditions  in  heat¬ 
recirculating  reactors  for  emissions  minimization  and  stable  range  extension. 

Further  details  of  counterflow  reactor  operation  at  lean  conditions  are  presented  in  Appendix  D. 

2.5  Effect  of  preheat  on  premixed  methane/air  flames 

In  the  filtration  reactor  the  analysis  of  the  conversion  of  a  fuel  to  syngas  is  complicated  by  the 
presence  of  the  porous  material.  In  order  to  better  understand  the  effect  of  preheat  on  conversion 
an  experimental  and  modeling  study  of  burner- stabilized  flames  was  conducted.  First 
experiments  were  conducted  with  unheated  reactants  so  that  confidence  in  the  models  could  be 
gained  for  their  prediction  of  combustion  characteristics,  especially  species  yields,  at  high 
equivalence  ratios.  The  effect  of  preheat  was  investigated  by  preheating  reactants  up  to  630  K 
while  holding  the  inlet  velocity  and  equivalence  ratio  constant. 

The  results  showed  that  the  burner  stabilized  flame  model  accurately  predicts  standoff  distance, 
limit  conditions  and  species  yields  for  experiments  with  flat  flame  burners.  Experiments  with 
preheated  reactants  demonstrated  that  the  operating  limits  for  a  flat  flame  burner  increase 
substantially  with  increased  reactant  temperature,  allowing  the  stable  operation  of  the  burner  at 
richer  equivalence  ratios — where  syngas  yields  are  greater — than  those  achievable  with  unheated 
reactants.  The  specific  effect  of  inlet  temperature  was  investigated,  and  it  was  found  that  inlet 
temperature  did  not  affect  the  yield  of  syngas  under  the  tested  range  of  temperatures.  Burner- 
stabilized  flames  with  preheated  reactants  showed  deviations  from  equilibrium  results  in  terms  of 
species  yields,  in  contrast  with  unheated  conditions  where  equilibrium  predicts  species  yields  for 
burner-stabilized  flames  very  well.  The  burner-stabilized  flame  model  was  able  to  predict  the 
deviation  from  equilibrium  when  reactants  were  preheated  to  -630  K.  At  high  preheat 
temperatures,  the  limiting  operating  conditions  of  the  experiment  were  well-predicted  by  the 
burner-stabilized  flame  model  but  not  by  the  laminar  flame  speed. 

These  results  have  important  implications  with  regards  to  heat  recirculating  reactors.  First,  the 
models,  and  therefore  the  kinetics  mechanism,  performed  very  well  up  to  richest  equivalence 
ratio  examined  of  1.75,  which  exceeds  the  conventional  rich  flammability  limit.  The  only 
significant  difference  between  the  models  and  the  experiments  was  the  difference  between  the 
limit  equivalence  ratio  and  the  equivalence  ratio  where  the  laminar  flame  speed  was  25  cm/s 
when  the  reactants  were  preheated  to  617  K.  This  means  that  those  working  on  reactor  models 
for  rich  combustion  should  have  confidence  in  GRI3.0  for  predicting  syngas  production  even 
beyond  the  rich  flammability  limit.  In  heat-recirculating  reactors,  preheat  is  accomplished  by 
transport  enhancement,  which  occurs  in  proximity  to  the  reaction  zone,  as  part  of  the  reactor 
design.  In  the  experiments  described  here,  the  preheat  occurred  far  upstream  of  the  reaction  zone 
and  was  input  from  an  external  source.  This  experimental  design  provided  a  means  to  study  the 
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effect  of  preheat  temperature  independent  of  enhanced  transport  mechanisms  and  interaction 
with  the  reaction  zone.  The  results  show  that  preheat  is  necessary  to  attain  stable  operation  at 
extreme  equivalence  ratios.  This  was  the  expected  result  since  conditions  (such  as  downstream 
propagation  of  the  reaction  zone  in  a  filtration  reactor)  that  enhance  preheat  have  been  observed 
in  heat-recirculating  reactors  when  mixtures  expected  to  be  unstable  are  burned.  It  was  also 
found  that  further  preheat  does  not  enhance  syngas  production.  This  result  aids  in  the  design  of 
heat  recirculating  reactors  because  it  indicates  that  maximal  heat  recirculation  may  not  be 
necessary  for  conversion  of  fuel  to  syngas. 

Further  details  of  the  study  of  the  effect  of  preheat  on  premixed  methane/air  flames  are  presented 
in  Appendix  E. 
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Appendix  A:  Non-catalytic  reforming  of  butanol  and  Jet-A 


The  first  objective  of  this  ARO-funded  project  was  to  demonstrate  that  Jet-A  and  butanol  could 
be  converted  to  syngas  through  non-catalytic  means.  This  work  is  described  below  and  was 
published  in  C.H.  Smith,  D.I.  Pineda,  C.D.  Zak,  and  J.L.  Ellzey,  J.  Hydrogen  Energy  38  (2013) 
879-889. 

In  a  filtration  reactor,  which  was  used  in  this  study,  a  reaction  zone  may  propagate  upstream  or 
downstream  relative  to  the  porous  solid  at  relatively  low  speed  (-1-10  cm/min).  This  type  of 
reactor  differs  substantially  from  other  types  of  porous  media  and  heat  recirculating  reactors  in 
which  the  front  is  stabilized  because  there  are  no  theoretical  flammability  limits  in  filtration 
reactors  [13].  Filtration  reactors  and  other  types  of  heat  recirculating  reactors  have  been  used  to 
convert  various  hydrocarbons  to  syngas,  particularly  gaseous  hydrocarbons  [14-21].  Conversion 
to  syngas  has  been  successfully  accomplished  with  liquid  fuels  including  ethanol  [22]  and 
heptane  [23]  in  a  filtration  reactor.  In  other  types  of  heat-recirculating  reactors,  diesel  [11,  24], 
octane  [25],  heptane  [26]  and  kerosene  (jet  fuel)  [27]  have  been  converted  to  syngas. 

Both  butanol  and  jet  fuel  have  been  studied  extensively  at  moderate  equivalence  ratios  (4>). 
Studies  on  butanol  include  basic  investigations  on  flames  [28],  kinetic  modeling  [29,  30], 
combustion  in  engines  [31],  and  catalytic  combustion  [32].  At  rich  conditions  where  the  interest 
is  on  the  production  of  syngas,  butanol  research  has  focused  on  equilibrium  thermodynamic 
analysis  [33-36],  and  conversion  by  catalytic  methods  [37-39]. 

At  moderate  (j),  there  has  been  significant  work  on  jet  fuel  and  other  high  molecular  weight  fuels 
over  the  last  two  decades.  Fundamental  work  on  diffusion  flames  [40,  41]  and  kinetic  modeling 
[42,  43]  has  significantly  advanced  our  understanding  of  these  complex  fuels,  as  has  work  on 
cracking  [44]  and  fuel  stability  [45].  Equilibrium  analyses  and  computational  studies  [46]  have 
shown  the  potential  for  syngas  production  from  jet  fuel  using  various  methods  including 
catalysts  [4,  47-54]  and  plasmas  [55].  Investigations  of  noncatalytic  reforming  in  porous  media 
have  generally  focused  on  reactors  with  stationary  reaction  zones  [25,  27]. 

Soot  formation  is  often  a  consideration  when  using  these  high  molecular  weight  fuels,  and  their 
propensity  to  soot  has  been  studied  in  various  experimental  configurations  including  diffusion 
flames  [42,  56,  57],  pool  fires  [58,  59],  and  engine  conditions  [60-62].  Soot  production  in 
premixed  flames  has  also  been  investigated  at  moderately  rich  (j)  [63,  64],  though,  as  noted  in  the 
literature  [65],  very  high  equivalence  ratios  are  difficult  to  study  because  of  flame  instability.  In 
filtration  combustion  soot  production  from  heptane  [23]  has  been  investigated.  In  other  porous 
media  reactors  soot  production  from  iso-octane  and  petrol  [25]  was  investigated  for  4>  up  to  4.25 
and  heptane,  diesel  oil,  kerosene  (jet  fuel)  and  bio-diesel  [27]  for  (J)  up  to  3. 
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As  a  demonstration  and  initial  investigation  of  the  principle  of  noncatalytic  conversion  of  jet  fuel 
and  butanol  to  syngas  by  transient  filtration  combustion,  we  performed  experiments  over  a  range 
of  rich  equivalence  ratios  and  inlet  velocities  (cold  gas).  Both  fuels  are  compared  for  syngas 
production  potential  and  soot  production,  and  a  comparison  to  other  fuels  is  performed  based  on 
data  found  in  the  literature. 

A.l  Principle  of  Operation  of  Filtration  Reactors 

In  filtration  combustion,  a  gaseous  mixture  flows  through  a  porous  solid  and  a  reaction  front  is 
established  within  the  pores.  The  fuel  may  be  either  the  solid  or  the  gas,  and  in  the  case  presented 
here,  the  solid  is  inert  and  the  combustible  gas  mixture  flows  through  the  porous  matrix.  Heat  is 
transferred  upstream  from  the  hot  products  to  the  reactants  through  the  porous  solid  by 
conduction  and  radiation.  Depending  on  the  operating  conditions,  the  reaction  front  may 
propagate  upstream  or  downstream  relative  to  the  solid.  This  propagation  speed  is  typically  much 
less  than  the  inlet  velocity  (V)  of  the  fuel/air  mixture  and  is  dependent  upon  the  properties  of  the 
solid,  the  properties  of  the  fuel  and  the  operating  conditions  [66]. 

When  the  front  is  propagating  upstream,  it  is  moving  into  a  relatively  cool  solid  resulting  in 
subadiabatic  temperatures,  i.e,  temperatures  lower  than  those  predicted  by  equilibrium. 
Conversely,  when  the  front  is  propagating  downstream  it  is  moving  into  a  preheated  solid, 
resulting  in  additional  heat  gain  to  the  reactants  and  local  superadiabatic  temperatures,  which  are 
above  those  predicted  by  equilibrium.  Figure  A.l,  which  is  a  computational  model  solution  for 
ethanol/air  [67],  describes  these  phenomena  by  showing  the  gas  temperature  as  a  function  of 
axial  position  for  four  different  times.  The  equilibrium  temperature  is  -1700  K,  and  the 
superadiabatic  temperature  spikes  can  be  clearly  seen  in  each  curve.  The  figure  also  shows  that 
the  preheat  region  expands  as  the  front  propagates  downstream. 
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Figure  A.l.  Representative  temperature  vs.  axial  position  for  a  reaction  zone  propagating 

downstream  in  a  filtration  reactor  [67] 

The  strongly  superadiabatic  temperatures  created  in  filtration  reactors  make  possible  the 
combustion  of  mixtures  far  outside  the  conventional  flammability  limits  [68].  Theoretical 
analysis  predicted  no  flammability  limit  for  filtration  combustion  [13],  and  experimental 
observations  have  confirmed  that  self-sustaining  combustion  is  achieved  at  both  very  rich  and 
very  lean  conditions  [69].  Under  conditions  when  the  reaction  zone  propagates  upstream  or 
downstream,  the  reactor  cannot  be  run  continuously  because  the  reaction  zone  eventually  either 
blows  off  or  flashes  back.  This  has  prompted  the  design  and  development  of  reactors  with 
stationary  reaction  zones  [9,  70].  In  contrast  to  reactors  with  propagating  reaction  zones,  reactors 
with  stationary  reaction  zones  have  narrower  operating  ranges,  only  producing  peak  temperatures 
that  are  slightly  superadiabatic.  Another  design,  called  a  reciprocal  flow  burner,  retains  the 
advantage  of  strongly  superadiabatic  temperatures  but  can  be  run  continuously  because  the 
reaction  zone  is  constrained  to  a  fixed  volume  in  a  porous  medium  by  reversing  the  flow  of  the 
reactants  periodically.  Filtration  combustion  has  been  the  topic  of  multiple  review  papers  [11, 
71-73],  and  the  detailed  principles  of  operation  can  be  found  in  the  literature  [66,  74,  75]. 

A.2  Experimental  Method 

Experimental  Apparatus 

The  experimental  apparatus  (Figure  )  consists  of  the  reactor,  the  fuel  vaporization  system,  the 
reactant  delivery  system,  and  the  data  acquisition  system.  A  quartz  cylinder  5  cm  in  diameter  and 
30  cm  in  length  includes  both  the  reactor  and  vaporization  chamber.  The  reactor  section  of  the 
quartz  cylinder  is  filled  with  99.5%  pure  alumina  (AL203)  spheres  of  5  mm  in  diameter  (Union 
Process)  and  insulated  on  the  outside  with  2  cm  thick  alumina.  The  alumina  insulation  has  holes 
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every  2  cm  along  the  length  to  allow  visual  observation  of  the  reaction  zone  location  and  soot 
deposition. 


Quartz  Probe 


The  vaporization  system  consists  of  an  air-atomizing  nozzle  and  the  quartz  chamber  for  mixing. 
The  atomized  fuel  droplets,  approximately  20  microns  in  size,  are  vaporized  in  a  separate  heated 
air  stream  in  the  quartz  chamber.  The  base  of  the  chamber,  which  holds  the  nozzle  and  air  heater, 
is  machined  and  welded  stainless  steel.  The  outer  surface  of  the  quartz  chamber  is  wrapped  with 
resistance  band  heaters  and  insulation  to  prevent  condensation  of  the  fuel.  Portions  of  the  band 
heaters  and  insulation  are  removable  so  that  any  potential  condensation  of  the  jet  fuel  can  be 
observed.  Some  operating  conditions  were  not  tested  (high  V  with  high  cj))  because  the 
vaporization  system  was  unable  to  produce  a  gaseous  stream  of  reactants.  Data  are  reported 
under  conditions  when  accumulation  of  liquid  fuel  was  not  observed,  and  when  the  reactant 
stream  was  always  transparent  (ie:  no  observable  droplets  of  liquid  fuel)  before  entry  into  the 
porous  media.  The  temperature  of  the  inlet  mixture  was  always  between  150°C  and  200°C  as 
measured  by  the  thermocouple  placed  before  entry  to  the  reactor. 


The  fuel  delivery  system  includes  a  fuel  tank,  a  pump  and  rotameters  for  fuel  regulation.  The 
stream  of  dry  laboratory  air  entering  the  nozzle  was  unheated,  but  an  inline  air  heater  heated  the 
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other  stream  of  air.  The  air  entering  the  nozzle  and  the  air  entering  the  inline  air  heater  were 
controlled  with  separate  mass  flow  controllers. 

For  these  experiments,  jet  fuel  from  Berry  Aviation  at  the  San  Marcos  Municipal  Airport  was 
used.  AVFuel  Corporation  supplies  the  fuel  to  Berry  Aviation.  Since  jet  fuel  can  have  a  varying 
composition,  an  average  chemical  formula  must  be  assumed.  In  this  work,  a  chemical  formula  of 
C11H21  was  used  [76].  The  density,  as  measured  by  a  laboratory  scale,  was  797  kg/m3,  and  the 
lower  heating  value  was  43.2  MJ/kg  as  given  by  [77].  The  butanol  that  was  used  in  this  study 
was  obtained  from  Superior  Solvents  and  Chemicals.  The  density  was  measured  as  810  kg/m3 
and  the  lower  heating  value  was  33.1  MJ/kg  [77]. 

Experimental  Procedure 

When  a  specific  operating  condition  (cj),  V)  was  to  be  tested,  the  reaction  front  was  first 
positioned  5-10  cm  from  the  bottom  of  the  porous  medium  part  of  the  reactor  by  igniting  a 
slightly  lean  flame  at  the  exit  of  the  porous  medium  and  allowing  it  to  propagate  upstream.  This 
procedure  ensured  that  any  soot  that  was  deposited  within  the  reactor  during  operation  at  the 
previous  condition  was  completely  burned  off  before  the  start  of  the  next  condition  to  be  tested. 
Then  the  flows  of  fuel  and  air  were  set  to  the  desired  values,  and  sampling  of  the  exhaust  gases 
by  a  Varian,  Inc.  gas  chromatograph  (GC)  began.  Each  GC  sample  required  160  seconds,  so 
samples  were  taken  when  the  reaction  zone  was  at  different  positions  along  the  length  of  the 
reactor  if  the  reaction  zone  was  propagating.  The  first  sample,  which  was  discarded  along  with 
the  second  sample  to  clear  the  line  between  the  reactor  and  the  GC,  was  taken  when  the  reaction 
zone  was  between  5  and  10  cm  from  the  bottom  of  the  porous  media.  The  location  of  the  reaction 
zone  was  monitored  visually  through  the  holes  in  the  insulation  surrounding  the  reactor,  and  no 
samples  were  taken  when  the  reaction  zone  was  less  than  ~5  cm  from  the  exit  of  the  reactor.  At 
least  4  samples  and  most  often  5  were  taken  for  each  operating  condition,  and  the  data  reported 
in  this  paper  are  averages  of  the  last  2  or  3  samples.  Given  this  procedure  all  reported 
measurements  were  taken  when  the  reaction  zone  was  between  10  and  15  cm  from  the  bottom  of 
the  reactor. 

The  GC  was  calibrated  to  measure  the  following  species:  hydrogen,  oxygen,  nitrogen,  methane, 
carbon  dioxide,  carbon  monoxide,  ethylene,  ethane,  acetylene  and  select  higher  hydrocarbons. 
The  GC  sampled  the  exhaust  gas  for  measurement  through  a  water-cooled  quartz  probe. 
Measurement  uncertainty  was  calculated  using  a  Student-t  distribution.  The  contributions  to  the 
uncertainty  were  uncertainty  in  the  GC  calibration  gases,  uncertainty  in  the  GC  calibration,  and 
uncertainty  in  the  flow  rates  of  air  and  fuel.  Uncertainties  for  all  calculations  were  based  on 
sequential  perturbation.  For  calculations  of  yield  and  carbon  balance,  the  molar  flow  rates  of 
species  in  the  exhaust  were  required.  These  molar  flow  rates  were  calculated  by  assuming  the 
flow  rate  of  nitrogen  into  the  reactor  equaled  the  flow  rate  of  nitrogen  out  of  the  reactor. 
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For  every  tested  condition,  the  first  two  measurements  were  discarded  and  the  rest  of  the 
measurements  were  averaged. 

A.3  Results  and  Discussion 

In  the  experiments  for  each  fuel,  we  obtained  one  set  of  data  by  varying  (|)  and  holding  V 
constant  and  another  by  varying  V  while  holding  (j)  constant.  The  experimental  operating 
conditions  are  described  in  Table  A.l. 


Fuel 

Equivalence  Ratio  (cf>) 

Inlet  Velocity  (V) 
[cm/s] 

Butanol 

1-5 

30 

Butanol 

3 

30-60 

Jet  A 

1-3.15 

34 

Jet  A 

2.7 

34-46 

Table  A.l:  Experimental  operating  conditions 

The  tested  ranges  for  (j)  and  V  were  determined  by  multiple  factors.  For  both  fuels,  maximum  V 
was  limited  by  the  rate  of  downstream  reaction  zone  propagation.  Species  data  are  reported  from 
conditions  when  the  GC  could  sample  the  exhaust  gases  at  least  4  times,  and  this  was  possible 
only  when  the  downstream  propagation  rate  was  less  than  ~2.5  cm/min.  For  butanol  and  jet  fuel, 
the  reaction  zone  remained  nearly  stationary  (the  reaction  zone  remained  within  5-10  cm  from 
the  bottom  of  the  reactor  in  the  time  that  5  samples  were  taken)  relative  to  the  porous  media 
when  the  operating  conditions  were  <j)  ~  2,  V  ~  32  so  the  GC  could  sample  the  exhaust  gases 
many  times.  When  (j)  or  V  increased  from  these  values,  the  reaction  zone  propagated  downstream 
up  to  ~3  cm/min,  and  the  speed  of  propagation  was  higher  for  jet  fuel  than  for  butanol  under  the 
near-overlapping  tested  conditions.  The  -2.5  cm/min  limit  on  the  propagation  speed  determined 
the  limiting  V  for  both  fuels  and  the  limiting  (|>  for  butanol,  however  soot  formation  limited  the 
maximum  testable  (j)  for  jet  fuel.  The  minimum  testable  V  was  determined  by  the  minimum  flow 
rates  required  for  safe  operation  of  the  inline  air  heater. 

Important  metrics  for  syngas  production  are  the  hydrogen  yield  and  the  carbon  monoxide  yield. 
These  metrics,  defined  below,  describe  how  effectively  the  reactor  converts  hydrogen  and  carbon 
bound  in  the  fuel  to  diatomic  hydrogen  and  carbon  monoxide,  respectively. 

hydro qen  yield  —  100  x  2xWhz  (A.l) 

J  ^  J  NHxNFuel  y 
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carbon  monoxide  yield  =  100  x - : - 

NcXNpuel 


(A. 2) 


where  the  units  of  NFuei,  NHz and  Nco  are  moles  per  second.  NH  =  10  and  Nc  =  4  for  butanol, 
Nh  =  21  and  Nc  =  11  for  jet  fuel  [76]. 


We  also  calculated  values  for  the  chemical  energy  conversion  efficiency.  This  metric  describes 
the  percentage  of  chemical  energy  in  the  original  fuel  that  was  contained  in  particular  exhaust 
species.  The  energy  conversion  efficiencies  were  calculated  using  the  LHV  of  the  individual 
species  as  described  by  the  following  representative  equation  for  hydrogen: 


hydrogen  energy  conversion  efficiency  —  100  x 


Nh2xLHVh2 
LHVFueix  N puei 


(A. 3) 


where  LHV  values  are  in  units  of  kJ  per  mole. 


Similarly,  we  calculated  the  total  energy  conversion  efficiency,  the  percentage  of  chemical 
energy  in  the  original  fuel  contained  in  all  fuel  exhaust  species.  If  all  species  are  at  ambient 
temperature  and  no  energy  is  lost  through  heat  transfer,  then  this  quantity  should  be  100%. 
Values  lower  than  100%  indicated  how  much  energy  is  associated  with  heat  losses  and  the 
sensible  energy  of  the  products.  The  total  energy  conversion  efficiency  was  calculated  as 
follows: 


total  energy  conversion  efficiency  =  100  x 


If NjXLHVj 
LH VpueiX  Npuei 


(A.  4) 


The  summation  is  performed  over  all  gaseous  species  in  the  exhaust  with  positive  LHV  and  with 
non-negligible  concentration  (concentration  >0.1%).  LHV  values  were  taken  from  [78]. 

For  comparison  with  experimental  results,  equilibrium  values  are  also  presented.  We  calculated 
these  values  for  a  constant  pressure/enthalpy  process  using  the  Cantera  software  suite  [79].  In 
heat-recirculating  reactors  peak  temperatures  can  significantly  exceed  the  adiabatic  flame 
temperature.  However,  there  is  no  external  energy  input  so  under  ideal  conditions  the  products 
exit  at  the  adiabatic  flame  temperature;  for  this  reason  it  is  appropriate  to  compare  the 
experimental  results  to  equilibrium  calculations  for  a  constant  enthalpy  process.  Previous  work 
[23,  80]  has  shown  that  equilibrium  calculations  predict  trends  in  exhaust  species  as  a  function  of 
reactant  composition  and  can  provide  insight  into  the  thermodynamic  character  of  the  conversion 
process.  All  equilibrium  calculations  included  the  formation  of  solid  carbon. 
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Soot  Observations 


One  of  the  challenges  of  investigating  butanol  and  jet  fuel  at  rich  conditions  was  the  propensity 
of  these  fuels  to  produce  soot.  For  both  fuels,  when  c|>  was  greater  than  approximately  1.5,  soot 
was  observed  in  the  reactor  through  the  holes  in  the  insulation  surrounding  the  reactor  and  at  the 
exit  of  the  reactor.  When  the  reaction  zone  propagated  downstream  soot  was  observed  both 
upstream  and  downstream  of  the  reaction  zone.  The  highest  (j)  tested  with  jet  fuel  was  3.58 
because  soot  laydown  at  this  equivalence  ratio  and  higher  was  so  great  that  it  clogged  the  pores 
of  the  reactor,  and  the  porous  media  was  ejected  from  the  reactor.  With  butanol  the  reactor  was 
not  clogged  at  any  tested  equivalence  ratio,  although  significant  soot  was  observed.  The 
quantities  of  soot  increased,  based  on  qualitative  observation,  with  increasing  (j). 

In  order  to  investigate  further  the  magnitude  of  soot  deposition,  a  carbon-flow  balance  was 
performed  from  the  experimental  measurements  and  equilibrium  calculations.  Figure  A. 2  shows 
this  carbon  balance  in  terms  of  the  ratio: 


carbon  out  _  measured  flow  rate  carbon  atoms  in  gaseous  exhaust 
carbon  in  flow  rate  of  carbon  atoms  in  fuel 


(A. 5) 


Figure  A.2:  Ratio  of  carbon  out  to  carbon  in  vs.  equivalence  ratio  with  inlet  velocity  =  32  + 
2  cm/s.  Lines  on  the  figure  indicate  approximate  equivalence  ratios  at  which  upstream  and 

downstream  propagation  occurred. 
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As  shown  in  the  figure,  equilibrium  predicts  that  all  carbon  is  in  gaseous  species  until  (|>  ~  3, 
above  which  significant  solid  carbon  formation  occurs.  This  carbon  formation  results  in  the 
relatively  sharp  decreases  in  CO  yield  and  total  energy  conversion  efficiency  seen  in  Figures 
A.4,  A. 5,  A.7,  A. 8,  A.10  and  A.ll.  In  the  experiments  near  stoichiometric  conditions,  carbon  is 
in  gaseous  compounds  for  both  butanol  and  jet  fuel.  As  4>  increases,  however,  the  ratio  of  carbon 
leaving  to  carbon  entering  decreases  steadily  for  jet  fuel  until  the  maximum  (})  tested  (3.15)  is 
reached.  This  indicates  that  carbon  is  being  converted  to  forms  not  measured  by  the  GC, 
including  large  hydrocarbons,  hydrocarbons  condensed  on  the  porous  media  and  solid  carbon 
deposited  on  the  porous  media.  For  jet  fuel  at  4>  =  3.15,  about  40%  of  the  carbon  entering  the 
system  is  not  measured  by  the  GC.  These  data  are  consistent  with  the  qualitative  observation  of 
heavy  soot  deposition  and  reactor  reactor  clogging  in  the  experiments  with  jet  fuel.  For  butanol 
the  ratio  decreases  until  4>  =  3,  and  then  rises  again  to  exceed  1  at  (j>  =  5. 

The  carbon  ratio  greater  than  1  is  counter-intuitive  and  initially  suggests  that  carbon  is  being 
“created.”  On  closer  inspection,  however,  this  is  a  consequence  of  the  transient  nature  of 
filtration  combustion.  As  shown  in  Figure  A.l,  when  the  front  propagates  downstream  three 
important  parameters  change  with  time:  the  residence  time  of  the  gas  upstream  of  the  reaction 
increases,  the  residence  time  of  the  gas  downstream  of  the  reaction  zone  decreases,  and  the 
preheat  region  expands.  These  effects,  coupled  with  the  soot  deposition  and  potential  large 
hydrocarbon  condensation  complicate  the  analysis  of  the  experimental  results  for  downstream- 
propagating  reaction  zones.  As  indicated  in  Figure  A. 2,  the  reaction  zone  propagated 
downstream  when  (j>  exceeded  2.5,  and  this  downstream  propagation  rate  increased  with 
equivalence  ratio. 

Soot  is  complex,  containing  a  variety  of  chemical  species  and  structures,  including  adsorbed 
hydrocarbons  [81].  Since  soot  is  composed  mainly  of  carbon  and  hydrogen,  it,  along  with 
adsorbed  hydrocarbons,  can  be  oxidized  to  other  species  that  are  measured  in  the  exhaust  [64, 
82].  One  would  expect  this  effect  to  be  particularly  pronounced  when  02  is  in  proximity  to 
surfaces  coated  with  soot  at  high  temperature,  and  this  occurs  significantly  when  the  reaction 
zone  propagates  downstream  as  soot  was  observed  upstream  of  the  reaction  zone.  Soot  particles 
can  grow  on  timescales  similar  to  the  residence  times  in  the  reactor  [83,  84].  Therefore  an 
additional  effect  of  downstream  propagation  is  that  hydrocarbon  species,  particularly  acetylene, 
that  contribute  to  soot  growth  have  less  time  to  do  so  because  their  residence  time  between  the 
reaction  zone  and  the  sampling  point  decreases  with  downstream  propagation.  The  carbon 
balance  exceeding  1  at  the  highest  (]>  is  suspected  to  be  the  result  of  the  soot 
consumption/oxidation  rate  exceeding  the  soot  deposition  rate  at  the  time  those  samples  were 
taken. 

The  data  that  were  taken  when  the  reaction  zone  was  nearly  stationary  (c[>  ~  2)  or  upstream 
propagating  (<|)  <  1.5)  are  considered  as  appropriate  representations  of  the  conversion  potential  of 
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these  fuels.  In  these  cases,  soot  is  not  observed  upstream  of  the  reaction  zone,  and  therefore  there 
is  no  potential  for  additional  fuel  species  to  be  added  to  the  gaseous  stream.  Results  in  the 
following  sections  are  presented  in  terms  of  yield  and  energy  conversion  efficiency,  which  are 
defined  relative  to  the  known  fuel  input  at  the  time  of  sample.  The  data  shown  in  Figure  A. 2 
imply  that  the  yields  and  energy  conversion  efficiencies  must  be  interpreted  with  the  carbon  flow 
balance  in  mind.  If  there  is  a  net  deposition  of  soot  at  the  time  a  given  sample  is  taken,  this  will 
be  reflected  in  lower  yields  and  energy  conversion  efficiencies.  If  the  rate  of 
oxidation/consumption  exceeds  the  deposition  rate  at  the  time  a  sample  is  taken  the  yields  and 
energy  conversion  efficiencies  will  be  artificially  high. 

Not  only  do  the  results  shown  in  Figure  A. 2  and  our  qualitative  observation  of  soot  (including 
complete  pore  blockage  in  the  jet  fuel  experiments)  have  implications  for  the  work  presented 
here,  but  these  data  also  have  important  implications  for  any  experimental  work  with  butanol  and 
jet  fuel  at  high  (j)  in  premixed  combustion.  Burning  the  soot  out  of  the  reactor  at  all  conditions 
was  accomplished  by  burning  a  lean  mixture  that  propagated  upstream  though  the  reactor.  For  a 
filtration  reactor  to  be  operated  in  any  practical  application  this  procedure  would  be  required.  In 
addition,  such  a  cleaning  procedure  would  likely  be  required  for  other  reactors  that  operate  with 
these  fuels  at  (|>  greater  than  about  2. 

Conversion  of  Jet  Fuel 

Ideal  fuel  conversion  results  in  the  complete  conversion  of  carbon  and  hydrogen  in  the  reactants 
to  CO  and  H2  in  the  products,  respectively.  The  yields  of  these  species,  then,  describe  reactor 
performance  relative  to  the  ideal  case.  The  effect  of  (j)  on  syngas  production  in  terms  of  H2  and 
CO  yields  is  shown  in  Figure  A. 3. 


Figure  A.3:  Yield  vs.  equivalence  ratio  for  jet  fuel  with  inlet  velocity  =  34  cm/s. 
Downstream  propagation  observed  for  equivalence  ratio  >  2.5 
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The  yields  for  CO  and  H2  are  zero  at  stoichiometric,  and  very  small  for  (j)  =  1.25.  At  higher  the 
equilibrium  values  are  all  higher  than  the  corresponding  experimental  measurements.  The  CO 
yield  increases  rapidly  with  (j)  until  (j)  =  2.25  while  the  H2  yield  continues  to  increase  over  the 
entire  tested  range.  The  maximum  yields  for  H2  and  CO  occur  at  different  (j)  and  are  42%  at  (j)  = 
3.15  and  56%  at  (j)  =  2.25,  respectively.  In  contrast,  equilibrium  calculations  predict  that  the 
peaks  in  conversion  of  H2  and  CO  occur  at  the  same  (j)  ((|)  ~  2.9). 

As  mentioned  previously,  the  energy  conversion  efficiency  is  the  percentage  of  chemical  energy 
in  the  fuel  that  is  contained  in  specific  species  in  the  exhaust.  The  total  energy  conversion 
efficiency  is  the  percent  of  the  chemical  energy  in  the  fuel  that  is  contained  in  the  exhaust.  Figure 
A.4  and 

Figure  A.5  show  the  total  energy  conversion  efficiency  and  species  energy  conversion  efficiency 
as  a  function  of  and  V  for  jet  fuel. 


Figure  A.4:  Energy  conversion  efficiency  vs.  equivalence  ratio  for  jet  fuel  with  inlet  velocity 
=  34  cm/s.  Downstream  propagation  observed  for  equivalence  ratio  >  2.5 


Figure  A.5:  Energy  conversion  efficiency  vs.  inlet  velocity  for  jet  fuel  with  equivalence  ratio 
=  2.7.  Downstream  propagation  observed  for  all  conditions 
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Similar  to  the  trend  observed  in  Figure  A. 3,  the  experimental  energy  conversion  efficiencies  for 
CO,  H2  and  the  total  are  significantly  lower  than  equilibrium  (Figure  A. 4).  The  total  energy 
conversion  efficiency  increases  with  (j)  up  to  a  maximum  of  61%  at  4>  =  3.15,  meaning  that  39% 
of  the  energy  in  the  original  fuel  is  lost  to  heat.  Besides  H2  and  CO,  acetylene,  and  to  a  lesser 
degree  methane,  are  the  only  energy-containing  species  that  contributes  substantially  to  the  total 
energy  conversion  efficiency.  This  is  in  contrast  to  equilibrium  which  predicts  almost  no 
hydrocarbons  in  the  products. 


Figure  A. 5  shows  equilibrium  and  experimental  energy  conversion  efficiency  as  a  function  of  V 
for  (j)  =  2.7.  The  experimental  data  show  that  V  does  not  affect  the  energy  conversion  efficiencies 
(species  specific  or  total)  to  a  large  degree  over  the  tested  range.  Previous  studies  have  shown 
that  the  energy  conversion  efficiency  can  depend  strongly  on  inlet  velocity;  however,  the  tested 
range  in  this  study  was  not  as  large  as  the  tested  ranges  reported  in  previous  studies. 

Conversion  of  Butanol 

The  potential  of  butanol  for  conversion  to  syngas  was  investigated  through  a  similar  set  of 
experiments  to  those  performed  with  jet  fuel.  Butanol  differs  substantially  from  jet  fuel  because 
it  has  a  smaller  molecular  weight  and  is  an  alcohol,  containing  an  oxygen  atom  in  the  molecule. 
In  Figure  A. 6,  the  yields  of  H2  and  CO  are  presented  as  a  function  of  (j)  with  V  held  constant  at 
30  cm/s. 


Figure  A.6:  Hydrogen  and  carbon  monoxide  vs.  equivalence  ratio  for  butanol  with  inlet 
velocity  =  30  cm/s.  Downstream  propagation  observed  for  equivalence  ratio  >  2.5 
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The  yields  for  H2  and  CO  increase  rapidly  as  c|>  is  increased  from  stoichiometric.  The  peak  CO 
yield  (72  %)  occurs  at  (j)  =  2,  and  the  peak  H2  yield  (43  %)  occurs  at  (j>  =  2.5.  As  cj>  increases  from 
these  values,  the  yields  for  CO  and  H2  decrease  slowly.  For  both  CO  and  H2,  the  yields  are  very 
similar  to  the  yields  predicted  by  equilibrium  until  cj>  ~  2.  At  higher  (j),  the  equilibrium  yields  are 
higher  than  the  experimentally  measured  yields,  though  at  (j)  =  5,  the  experimental  value  and 
equilibrium  values  for  CO  yield  are  nearly  identical. 

Figure  A. 7  and  Figure  A. 8  shows  the  energy  conversion  efficiency  for  the  measured  species  and 
the  total  energy  conversion  efficiency  as  a  function  of  (j)  and  V,  respectively.  Both  experimental 
data  and  equilibrium  results  are  shown.  In  Figure  A.7  the  total  energy  conversion  efficiency 
determined  experimentally  follows  the  equilibrium  value  until  <|>  >  3  (where  relatively  fast 
downstream  propagation  was  observed).  As  (|>  increases  beyond  3,  the  experimental  value  for  the 
total  energy  conversion  efficiency  increases  until  it  slightly  exceeds  100%.  As  explained  above, 
it  is  suspected  that  an  imbalance  in  soot  deposition  and  consumption  at  the  time  of  exhaust  gas 
sampling  is  the  reason  for  the  unrealistically  high  energy  conversion  efficiencies  observed  at  4*  > 
3.5. 


Figure  A.7:  Energy  conversion  efficiency  vs.  equivalence  ratio  for  butanol  with  inlet 
velocity  =  30  cm/s.  Downstream  propagation  observed  for  equivalence  ratio  >  2 
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Inlet  Velocity  (cm/s) 

Figure  A.8:  Energy  conversion  efficiency  vs.  inlet  velocity  for  butanol  with  equivalence 
ratio  =  3,  with  downstream  flame  front  propagation  observed  for  all  conditions 

For  conversion  of  butanol  to  H2  and  CO,  lower  (j)  are  optimal  since  the  conversion  of  energy 
from  butanol  to  H2  and  CO  peaks  at  (j)  ~  2.  At  (j)  ~  2,  a  significant  amount  (10%)  of  the  chemical 
energy  is  in  acetylene,  and  the  contribution  of  acetylene  to  the  chemical  energy  outflow  increases 
to  a  maximum  of  20%  at  (j)  =  3.25  and  then  decreases  as  c|)  increases.  Methane  and  ethylene  also 
contributed  significantly  to  the  chemical  energy  outflow  as  (j)  exceeded  2.  The  energy  in  methane 
peaks  at  (j)  =  4  with  a  maximum  of  23%,  while  the  energy  in  ethylene  is  nonzero  at  (j)  =  3  and 
rises  steadily  with  (j)  up  to  a  maximum  of  36%  at  (j)  =  5.  These  results  contrast  with  those  for  jet 
fuel  (Fig.  A. 5)  which  shows  significantly  less  contribution  from  hydrocarbons  to  the  energy 
content  of  the  exhaust. 

Figure  A.8  shows  the  energy  conversion  efficiencies  as  a  function  of  V  for  butanol  with  (j)  =  3. 
Similar  to  the  results  for  jet  fuel,  the  energy  conversion  efficiencies  for  H2  and  CO  do  not 
change  significantly  as  V  is  increased,  but  the  total  energy  conversion  efficiency  does  increase 
somewhat  as  V  exceeds  30  cm/s.  Exceeding  the  equilibrium  value  for  total  energy  conversion 
efficiency  and  having  near  100%  conversion  when  the  V  was  40  cm/s  and  50  cm/s  is,  again, 
likely  the  result  of  soot  deposition/consumption  imbalance  and  reaction  zone  propagation 
(discussed  below). 
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Fuels  Comparison 


Thermodynamic  analysis  is  often  used  to  investigate  fuel  conversion  to  syngas  [34,  70,  85].  In 
order  to  provide  insight  into  the  relative  syngas  production  efficiencies  of  a  variety  of  fuels, 
equilibrium  calculations  were  completed  for  a  set  of  seven  fuels  with  a  range  of  molecular 
weights.  The  computed  hydrogen  yield  and  carbon  monoxide  yield  as  a  function  of  (j)  for  these 
fuels  is  shown  in  Figure  A. 9  and  Figure  A.  10,  respectively. 
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Figure  A.9:  Hydrogen  yield  vs.  equivalence  ratio  (equilibrium  values) 


Q) 

o 

o 


3°1.5 


*-  octane 
«■  jet  fuel 
*■  ethanol 
o  butanol 


2  2.5  3  3.5  4  4.5 

Equivalence  Ratio 

Figure  A. 10:  Carbon  monoxide  yield  vs.  equivalence  ratio  (equilibrium  values) 


The  equilibrium  data  show  a  trend  of  increasing  H2  yield  with  increasing  fuel  molecular  weight. 
The  trend  is  the  same  for  CO  yield  excluding  the  behavior  of  the  oxygenated  fuels  at  (j)  >  3. 
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Based  on  equilibrium,  jet  fuel  is  expected  to  produce  the  highest  yields  of  H2  and  CO,  while 
ethanol  should  produce  the  lowest  yields.  Amongst  the  non-oxygenated  fuels,  methane  should 
produce  the  lowest  yields. 

In  earlier  publications,  the  authors’  group  reported  results  of  experiments  converting  methane 
[86],  heptane  [23],  and  ethanol  [22]  to  syngas  by  filtration  combustion.  The  filtration  reactors 
that  were  used  in  these  studies  were  similar  to  that  used  in  the  current  study,  except  that  for 
methane  where  the  reactor  consisted  of  ceramic  foam  rather  than  a  bed  of  alumina  spheres.  For 
each  fuel,  the  operating  conditions  were  different,  so  direct  comparison  across  a  range  of  V  and  (j) 
is  not  possible.  However,  there  are  (j),  V  pairs  that  overlapped  or  near-overlapped  between  each 
of  the  fuels  studied,  so  a  limited  comparison  can  be  made  amongst  the  fuels  in  terms  of  hydrogen 
and  carbon  monoxide  yields.  These  data  can  be  used  to  determine  the  accuracy  of  using 
equilibrium  data  to  predict  relative  fuel  performance.  The  hydrogen  and  carbon  monoxide  yields 
for  given  (j),  V  pairs  are  listed  in  Table  A. 2. 
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Table  A.2:  Comparison  of  hydrogen  and  carbon  monoxide  yields  for  various  fuels 
(Representative  uncertainty  in  yield  values:  ±5%,  representative  uncertainty  in  equivalence 
ratio:  0.1,  representative  uncertainty  in  inlet  velocity:  ±2  cm/s.  Data  for  ethanol  from  [22], 
data  for  methane  from  [86],  and  data  for  heptane  from  [23]) 

CO  yield  was  not  reported  in  the  experiments  with  methane.  However,  H2  yield  can  be 
compared  with  the  various  fuels.  As  the  tables  show,  the  H2  yield  for  methane  exceeded  the  H2 
yield  that  was  observed  with  each  of  the  other  fuels.  In  terms  of  H2  and  CO  yield,  ethanol  is  next 
best,  having  a  slightly  higher  yield  of  both  H2  and  CO  than  butanol  and  heptane  under  the  same 
operating  conditions.  Heptane  has  the  next  highest  conversion  with  yields  higher  than  butanol. 
Jet  fuel  is,  by  far,  the  least  efficient  producer  of  syngas  with  both  H2  and  CO  yields  less  than 
those  of  all  the  other  fuels  except  for  the  H2  yield  compared  with  butanol.  Ten  binary 
comparisons  can  be  made  amongst  the  fuels  in  terms  of  H2  yield,  and  six  can  be  made  in  terms 
of  CO  yield  with  these  data.  Among  the  ten  H2  yield  comparisons  and  six  CO  yield  comparisons 
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only  in  two  cases  does  a  fuel  with  a  larger  molecular  weight  have  a  greater  hydrogen  yield  than 
one  with  lower  molecular  weight  (heptane  vs  butanol  for  both  H2  and  CO).  These  data,  then, 
suggest  that  syngas  production  efficiency  generally  decreases  with  molecular  weight. 
Interestingly,  this  is  the  opposite  of  what  is  predicted  by  thermodynamic  equilibrium  (Figure  A. 9 
and  Figure  A.  10). 

Experimental  investigations  of  conversion  of  rich  hydrocarbon  mixtures  have  been  conducted  by 
other  groups  in  the  last  few  years,  and  these  data  allow  a  further  comparison  of  the  relative 
potential  of  fuels  to  be  converted  to  syngas.  Pedersen-Mjaanes,  Chan,  and  Mastorakos  [87] 
investigated  the  conversion  of  methanol,  methane,  octane  and  petrol  to  syngas  in  a  2-section  heat 
recirculating  reactor.  They  found  methanol  to  have  the  highest  maximum  yield  followed  by 
methane  and  then  octane  and  petrol  [87].  Pastore  and  Mastorakos  performed  experiments  with 
heptane  and  diesel  in  a  2-section  reactor  [24],  and,  in  contrast  to  the  trends  observed  in  the  work 
of  the  authors’  group,  they  reported  diesel  to  have  a  significantly  higher  rate  of  conversion  than 
heptane  in  terms  of  the  conversion  of  energy  in  fuel  to  energy  in  both  H2  and  CO.  Kennedy  and 
coworkers  at  the  University  of  Illinois  at  Chicago  performed  experiments  with  methane,  ethane, 
and  propane  in  a  filtration  reactor  very  similar  to  the  one  used  in  the  present  experiments  [16]. 
Their  experiments  showed  nearly  identical  CO  yields  for  each  of  the  three  fuels,  but  maximum 
H2  yields  increased  in  order  of  decreasing  molecular  weight.  In  a  recent  paper,  the  relative 
conversion  efficiencies  for  experiments  with  methane,  propane  and  heptane  in  a  counterflow 
reactor  were  compared  [26].  In  those  experiments,  propane  had  the  highest  efficiency  of 
conversion  to  H2,  followed  by  methane  and  then  heptane.  With  these  data  from  the  literature  ten 
binary  comparisons  can  be  made  in  terms  of  H2  and  CO  yield.  Of  these  ten  comparisons,  three 
are  in  agreement  with  equilibrium  prediction  (yield  increase  with  molecular  weight),  five  suggest 
that  yield  decreases  with  molecular  weight,  and  two  support  neither  trend.  As  with  the  filtration 
combustion  experiments,  these  data  suggest  that  the  trend  in  syngas  production  efficiency  as  a 
function  of  fuel  size  is  not  predicted  by  equilibrium. 

Based  on  the  data  reported  in  all  of  these  studies,  including  the  present  study,  the  potential  of 
fuels  to  be  converted  to  syngas  roughly  decreases  with  increasing  molecular  weight.  Currently, 
the  available  literature  does  not  provide  enough  data  to  draw  strong  conclusions  about  fuel 
optimization  for  conversion  to  syngas,  and  there  are  multiple  exceptions  to  this  general  trend 
indicating  that  there  is  a  need  for  a  comprehensive  study  of  the  effect  of  fuel  on  conversion 
where  fuels  are  tested  under  identical  conditions.  The  major  conclusion  of  this  fuels  comparison 
is  that  equilibrium  predictions  of  syngas  production  as  a  function  of  fuel  size  are  inaccurate 
under  most  conditions.  Since  equilibrium  is  often  used  to  predict  general  trends  in  syngas 
production,  these  data  suggest  that  care  must  be  taken  when  interpreting  the  results  of 
equilibrium  computations. 
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A.4  Conclusions 


Results  of  the  conversion  of  butanol  and  jet  fuel  to  syngas  via  non-catalytic  transient  filtration 
combustion  are  presented.  At  maximum,  about  42%  of  the  hydrogen  in  jet  fuel  was  converted  to 
H2  and  56%  of  carbon  was  converted  to  CO.  The  H2  yield  continued  to  increase  with  (j)  in  the  jet 
fuel  experiments  whereas  in  the  butanol  experiments  the  yields  for  H2  and  CO  both  reached 
peaks  within  the  tested  operating  range.  The  peak  CO  yield  for  experiments  with  butanol  was 
72%  and  the  peak  H2  yield  was  43%.  Most  of  the  chemical  energy  exiting  the  reactor  was  bound 
in  H2  and  CO,  but  CH4,  C2H2  and  C2H4  were  observed  in  substantial  amounts  in  experiments 
with  both  fuels,  especially  in  experiments  with  butanol  at  high  <j).  In  contrast,  equilibrium 
predicted  nearly  negligible  amounts  of  energy-containing  compounds  besides  H2,  CO  and  solid 
carbon. 

Soot  production  proved  to  be  substantial  in  experiments  with  both  fuels,  with  the  qualitative 
observation  of  soot  from  both  fuels  exceeding  the  soot  production  observed  in  similar 
experiments  with  methane,  heptane,  and  ethanol.  Soot  production  limited  the  maximum  testable 
cj)  for  jet  fuel  experiments  to  <j)  =  3.15  since  soot  clogged  the  reactor  when  (j)  exceeded  this  value. 
An  analysis  of  the  carbon  balance  for  jet  fuel  experiments  showed  that  up  to  40%  of  the  carbon 
entering  the  system  was  deposited  on  the  porous  medium.  Similarly  for  butanol,  the  carbon 
balance  showed  significant  loss  up  to  (j)  ~  3.  At  the  most  extreme  4>,  the  reaction  zone  propagates 
downstream  oxidizing  some  of  the  soot  that  has  been  deposited  previously,  resulting  in  an 
exhaust  stream  with  more  carbon  atoms  than  are  contained  in  the  inlet  fuel  at  the  time  of  sample. 

Equilibrium  calculations  and  experimental  results  were  analyzed  to  determine  trends  in  fuel 
conversion  efficiency  as  a  function  of  fuel  type.  The  equilibrium  calculations  predicted  that 
larger  fuel  molecules  would,  generally,  produce  higher  syngas  yields,  but  the  experimental 
results  disagreed  with  this  trend.  Generally,  experimental  results  show  that  smaller  hydrocarbons 
produce  higher  syngas  yields  with  jet  fuel  having  the  lowest  hydrogen  and  carbon  monoxide 
yields.  The  comparison  of  these  equilibrium  data  with  the  experimental  data  available  in  the 
literature  suggests  that  kinetic  effects  are  very  important  in  fuel  conversion  to  syngas  and 
equilibrium  results  should  be  interpreted  carefully  when  predicting  syngas  production  from  the 
combustion  of  very  rich  mixtures  of  fuel  and  air.  Lastly,  the  demonstration  of  the  conversion  of 
jet  fuel  and  butanol,  along  with  previous  demonstrations  with  methane,  ethanol  and  heptane 
show  the  robustness  of  the  non-catalytic  filtration  reactor. 
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Appendix  B:  Reforming  of  heptane  in  counterflow  reactor 


The  second  objective  of  this  ARO-funded  project  was  to  demonstrate  liquid  fuel  reforming 
capability  in  the  counterflow  reactor.  This  work  is  described  below,  and  was  published  as  E.L. 
Belmont,  S.M.  Solomon,  J.L.  Ellzey,  Combust.  Flame  159  (2012)  3624-3631 . 

One  of  the  largest  impediments  to  the  implementation  of  hydrogen-based  portable  power 
technologies,  such  as  fuel  cells,  is  the  lack  of  infrastructure  for  the  production  and  distribution  of 
hydrogen.  [1],  Distribution  of  hydrogen  entails  significant  challenges,  many  of  which  can  be 
overcome  if  hydrogen  is  produced  at  the  site  of  use  from  fuels  that  are  more  readily  transported. 
These  include  natural  gas  and  liquid  hydrocarbon  fuels,  the  latter  of  which  have  high  volumetric 
energy  density  [2,3].  The  ability  to  reform  liquid  fuels  on  a  small  scale  and  at  the  site  of  use 
therefore  has  the  potential  to  replace  or  supplement  batteries  for  energy  storage. 

Various  reactors  have  employed  the  concept  of  heat  recirculation  through  a  solid  matrix,  with 
some  designs  having  a  means  of  stabilizing  the  reaction  front.  A  two  section  porous  media 
reactor  can  support  a  stationary  combustion  zone  and  this  design  has  been  used  to  reform  rich 
mixtures  of  methanol,  methane,  octane  and  automotive-grade  petrol  [4].  Alternatively,  a 
diverging  section  of  porous  media  permits  stabilized  reaction  zones  for  fuel  reforming  in  porous 
media  and  has  been  used  successfully  in  several  studies  [5,6,7].  In  some  heat  recirculating 
reactors,  the  reaction  front  is  not  stabilized  and  propagates  through  the  solid  matrix.  When  the 
front  propagates  downstream,  it  passes  over  preheated  solid  and  the  gas  gains  additional  energy, 
resulting  in  superadiabatic  temperatures.  These  higher  temperatures  increase  reaction  rates  and 
potentially  result  in  high  conversion  efficiencies  [8,9].  A  disadvantage  of  this  design  is  that  the 
wave  ultimately  reaches  the  end  of  the  reactor  requiring  periodic  restarting  or  reversal  of  the 
flow  [10]. 

The  advantage  of  superadiabatic  temperatures  has  been  realized  in  other  reactor  designs  as  well. 
A  Swiss  roll  reactor  [11,12,13]  has  been  used  for  both  reforming  and  thermal  oxidation  [14].  In 
this  reactor,  heat  from  the  exhaust  is  transferred  through  the  walls  to  the  incoming  reactants. 
This  basic  principle  was  applied  in  the  study  of  the  counterflow  reactor  consisting  of  straight 
channels  which  simplify  the  geometry  and  permit  a  wide  range  of  stable  operating  conditions. 
Figure  B.l  is  a  schematic  of  the  reactor.  The  flow  direction  alternates  in  adjacent  channels  so 
that  energy  from  hot  combustion  products  is  transferred  to  cold  incoming  reactants.  The 
preheated  reactants  may  reach  temperatures  in  the  combustion  zone  that  are  superadiabatic.  The 
counterflow  reactor  design  has  been  proposed  analytically  [15,16,17]  and  validated  in  two 
experimental  studies  on  the  reforming  of  methane  [18]  and  propane  [19]. 
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Figure  B.l:  Counterflow  reactor  utilizes  heat  transfer  from  hot  combustion  products 
through  channel  walls  to  preheat  incoming  reactants 


One  of  the  compelling  characteristics  of  non-catalytic  reactors  is  their  ability  to  operate  on 
various  fuels.  There  have  been  numerous  studies  on  the  production  of  syngas  from  gaseous  fuels 
[4,6,20,21,22],  but  fewer  have  done  so  from  liquid  fuels.  In  this  study  we  established  the  first 
liquid  fuel  reforming,  using  heptane  fuel,  in  the  counterflow  reactor.  The  results  augment  those 
obtained  on  an  identical  reactor  using  methane  [18]  and  propane  [19].  Heptane  is  a  particularly 
interesting  liquid  fuel  as  a  surrogate  for  commercial  fuels.  Dixon  et  al.  [23]  presented 
experimental  and  numerical  analyses  of  syngas  production  from  heptane  in  a  porous  media 
reactor.  Equivalence  ratio  (<j))  and  inlet  velocity  (u)  were  varied,  and  exhaust  gas  hydrogen 
concentrations  were  found  to  increase  with  both  of  these  parameters.  These  trends  are  similar  to 
those  found  in  previous  filtration  combustion  studies  [4,21].  Pastore  and  Mastorakos  [24] 
reformed  heptane  to  syngas  in  a  two-section  porous  burner  and  obtained  comparable  results  to 
Dixon  et  al.  [23]  at  some  operating  conditions  but  with  a  more  limited  stable  burner  range. 


The  earlier  study  of  the  porous  media  reactor  operating  on  heptane  provided  guidance  for  this 
work.  In  that  study  [20],  the  reactor  operated  near  the  conventional  flammability  limits  of 
heptane  without  prohibitive  soot  formation  causing  any  observed  degradation  in  performance.  At 
more  extreme  conditions,  soot  buildup  did  pose  a  challenge,  and  Dixon  et  al.  [23]  developed  a 
regeneration  procedure  that  is  utilized  in  this  study  as  well.  The  non-catalytic  surfaces  of  the 
counterflow  reactor  imply  that  small  amounts  of  deposits  should  not  degrade  the  performance  as 
long  as  the  channel  dimensions  and  wall  properties  are  not  significantly  changed  due  to  buildup, 
thus  opening  the  possibility  of  using  fuels  that  may  produce  more  soot. 


In  this  work,  we  focused  on  conversion  of  heptane  to  syngas  in  the  counterflow  reactor.  The 
operating  limits  and  conversion  efficiencies  were  determined,  and  these  quantities  were 
compared  to  those  for  the  same  reactor  operating  on  methane  and  propane. 

B.l  Experimental  method 

Experimental  apparatus 

Figure  B.2  shows  the  counterflow  reactor,  which  is  the  same  design  used  in  previous  studies  on 
reforming  of  methane  and  propane  [18,19].  It  consists  of  four  4  mm  high  parallel  channels 
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constructed  of  1  mm  thick  silicon  carbide  (SiC)  walls.  The  reactor  is  17.3  cm  long,  with  a  91.5 
mm  long  main  section  and  two  40.8  mm  long  reactor  heads  where  inlets  and  outlets  are  located. 
Also  located  in  the  reactor  heads  are  end  plugs,  which  prevent  the  mixing  of  combustion 
products  and  unreacted  fuel/air  mixture,  and  two  6.4  mm  wide  sections  of  SiC  porous  foam  (17.7 
pores/cm,  9%  density)  which  act  as  flow  straighteners  and  flame  arresters.  Two  alumina  walls 
spaced  33.6  mm  apart  enclose  the  SiC  channels.  Fine  wire  B-type  thermocouples  are  inserted 
into  the  channels  through  the  alumina  insulation  10  mm  apart.  They  sit  flush  with  the  channel 
wall  in  order  to  prevent  flame  holding.  The  whole  reactor  is  shrouded  in  alumina  insulation  to 
minimize  heat  losses.  Further  details  about  the  reactor  construction  can  be  found  in  Schoegl  and 
Ellzey  [18]. 


6.4  mm  wide  10  mm  spaced 

Plugs  Inlets  Ceramic  Foam  Thermocouples  Insulation  Ceramic  Walls 


■IX! 


Outlets 


4  mm 
Channel 
Height 


40.8  mm  I  91.5  mm  I  40.8  mm 

Reactor  Head  Main  Reactor  Section  Reactor  Head 

Figure  B.2:  Counterflow  reactor 


n-Heptane  (Fisher  Scientific,  HPLC  grade)  is  pumped  to  an  atomizing  nozzle  mounted  inside  a 
mixing  chamber.  Air  is  fed  to  the  mixing  chamber  using  a  Hastings  mass  flow  controller. 
Downstream  of  the  controller  the  total  air  is  split  between  the  atomizing  nozzle  and  heater  air  by 
a  control  valve.  The  heated  air  is  introduced  into  the  annular  space  surrounding  the  nozzle  in 
order  to  vaporize  the  atomized  heptane  and  prevent  condensation  of  heptane  on  the  mixing 
chamber  walls.  The  vaporized  fuel/air  mixture  is  fed  into  a  manifold  that  splits  the  flow  among 
the  four  channels.  The  delivery  system,  from  the  mixing  chamber  to  the  inlet  ports,  is  wrapped  in 
heat  tape  and  monitored  using  K-type  thermocouples  to  maintain  the  system  above  150  °C. 


Exhaust  gas  composition  is  analyzed  using  a  Varian  CP  4900  gas  chromatograph  (GC)  with  three 
columns:  a  molecular  sieve  (Molsieve)  measures  diatomic  hydrogen  (H2),  nitrogen  (N2),  oxygen 
(02),  carbon  monoxide  (CO)  and  methane  (CH4);  a  porous  polymer  unit  (PPU)  gives 
concentrations  of  carbon  dioxide  (C02),  ethylene  (C2H4),  ethane  (C2H6),  acetylene  (C2H2)  and 
propane  (C3H8);  and  a  CP-SIL  column  detects  larger  hydrocarbons  (i-,n~)  butane,  (i-,n-) 
pentane,  n -hexane  and  n- heptane.  In  order  to  quench  reactions  outside  of  the  reactor  and  prevent 
moisture  accumulation  in  the  GC,  exhaust  samples  are  drawn  from  the  reactor  exhaust  through  a 
quartz  probe.  The  sample  then  passes  through  inert  Silco  Steel  tubing  and  a  filter  before  entering 
the  GC.  Previous  findings  [19]  showed  slightly  lower  H2  and  CO  product  concentrations  from 
outside  channels  and  minimal  variation  of  composition  with  position  within  a  single  channel 
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except  near  the  alumina  walls.  Exhaust  samples  in  this  study  were  taken  from  the  center 
channels,  which  are  expected  to  be  most  representative  of  interior  channels  in  scaled  up  reactors, 
and  the  sample  probe  was  placed  at  the  middle  of  the  channel  exit. 

Experimental  method 

The  start-up  procedure  began  with  warming  the  system  to  above  150  °C  using  both  heating  tape 
and  heated  air  flowing  through  the  reactant  delivery  system  and  the  reactor.  A  slightly  rich, 
vaporized  premixed  heptane/air  mixture  was  then  introduced  to  the  system  and  ignited  at  the 
channel  outlets  at  an  inlet  velocity  of  50  cm/s.  Flame  fronts  propagated  upstream  in  the  reactor 
channels  until  they  stabilized  downstream  of  the  porous  SiC  flow  straighteners.  Operating 
conditions  were  gradually  adjusted  to  (j)  =  3.0  and  u  =  125  cm/s,  where  the  reactor  was  operated 
until  temperatures  stabilized.  The  initial  warm-up  phase  took  approximately  30-45  minutes. 

Experiments  included  determination  of  stable  operating  conditions,  and  studying  the  effect  of 
inlet  velocity  and  equivalence  ratio  variation  on  exhaust  gas  composition.  Test  conditions  were 
considered  stable  if  combustion  could  be  sustained  in  the  main  reactor  section  for  ten  minutes. 
Alumina  reactor  wall  temperatures  in  each  of  the  four  channels  were  monitored  throughout  the 
experiments.  Gas  temperatures  were  not  directly  measured  due  to  the  tendency  of  thermocouples 
inserted  into  the  channels  to  act  as  flame  holders;  however,  an  analytical  model  of  the 
counterflow  reactor  predicts  peak  gas  temperatures  significantly  above  the  maximum  wall 
temperatures  [20].  Peak  temperatures  were  recorded  after  the  reactor  stabilized  at  a  given  set 
point,  and  showed  less  than  ±  5  °C  variation  during  the  time  of  continuous  operation  at  that  set 
point  and  while  exhaust  gas  samples  were  taken.  All  data  presented  in  this  study  were  taken  at 
operating  points  that  met  the  above  criterion  for  stability,  and  no  directional  drift  in  measured 
temperatures  or  species  concentrations  was  observed  with  repeated  sampling. 

Inlet  velocity  and  equivalence  ratio  were  changed  in  steps  of  Au  =  15  cm/s  and  Acj)  =  0.1, 
respectively,  and  the  inlet  velocity  was  specified  at  standard  conditions  of  25  °C  and  1  atm. 
Results  were  obtained  for  variation  of  inlet  velocity  from  50-200  cm/s  with  equivalence  ratio 
held  constant  at  3.0,  and  variation  of  equivalence  ratio  from  2. 9-3. 8  with  inlet  velocity  held 
constant  at  125  cm/s.  The  inlet  velocity  range  was  restricted  at  the  upper  limit  in  order  to  avoid 
damage  to  the  reactor  due  to  wall  temperatures  approaching  1300  °C,  while  the  lower  limit  was 
determined  by  flame  instability  and  extinction.  The  equivalence  ratio  range  was  limited  by 
flashback  on  the  lean  end  of  the  tested  range,  in  which  the  flame  propagated  upstream  into  the 
porous  SiC  foam,  and  blowoff  on  the  rich  end,  where  the  flame  front  propagated  out  the  end  of 
the  reactor.  Significant  deposits  of  soot  and  pyrolytic  graphite  [25]  were  observed  during 
extended  periods  of  operation  at  the  highest  tested  equivalence  ratios  and  inlet  velocities.  For  all 
experimental  conditions  reported  in  this  paper,  the  soot  build-up  was  prevented  by  briefly 
interrupting  the  fuel  and  purging  the  reactor  with  air  for  10-20  seconds  before  continuing  the 
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experiments.  At  high  equivalence  ratio  and  inlet  velocity  operating  points,  where  soot  buildup 
was  most  likely,  the  reactor  was  purged  approximately  every  20  minutes. 

At  each  operating  point,  five  consecutive  GC  measurements  were  taken  and  the  first  two  were 
discarded  to  ensure  that  any  residual  gases  in  the  sample  line  had  been  eliminated.  Each 
operating  point  was  tested  twice,  giving  a  total  of  six  data  points.  Total  uncertainty  was 
calculated  at  each  operating  point  as  the  root-sum-square  of  the  contributing  uncertainties. 
Therefore  the  uncertainty  of  GC  results  is  due  to  repeatability  of  measurements,  calculated  using 
a  student-t  distribution,  calibration  gases,  with  specified  uncertainties  of  1-5%,  and  the  GC  itself 
with  uncertainty  of  1%  of  the  maximum  calibrated  value  for  each  species.  The  error  bars  shown 
in  plotted  results  represent  average  uncertainties  for  each  set  of  operating  conditions.  Inlet 
velocity  uncertainty  is  due  to  uncertainties  of  1%  of  the  air  flow  controller  maximum  flow  rate 
and  2%  of  the  rotameter  maximum  flow  rate,  as  well  as  an  estimated  5%  tolerance  of  the  channel 
cross-sectional  area.  Equivalence  ratio  uncertainty  is  due  to  contributions  from  the  air  flow 
controller  and  rotameter.  Inlet  velocity  uncertainty  increases  with  increasing  u ,  where  Su  =  ±4.0 
cm/s  at  u  =  50  cm/s  and  du  =  ±10.3  cm/s  at  u  =  200  cm/s.  Equivalence  ratio  uncertainty  increases 
slightly  with  increasing  (j>  but  more  significantly  with  decreasing  u,  with  <-)'(|)  =  ±0.05  at  u  =  200 
cm/s  and  =  ±0.20  at  u  =  50  cm/s. 

B.2  Results  and  discussion 

Operating  Range 

One  of  the  critical  aspects  of  the  counterflow  reactor  is  the  ability  to  recirculate  heat  from  the  hot 
products  to  the  incoming  reactants  in  the  neighboring  channel.  The  preheating  increases  the 
burning  rate  for  a  reactant  mixture  compared  to  standard  conditions.  The  combustion  zone 
stabilizes  at  a  position  in  the  reactor  where  the  burning  rate  is  balanced  by  the  inlet  velocity.  For 
example,  increasing  the  equivalence  ratio  from  3.2  to  3.6  at  a  fixed  velocity  results  in  lower  heat 
release  and  the  stable  position  for  the  flame  is  further  downstream  in  the  channel  where  the  heat 
recirculation  is  more  significant.  Ultimately  as  the  equivalence  ratio  is  increased  further,  heat 
release  and  recirculation  are  insufficient  compared  to  heat  losses  and  the  flame  blows  off.  If 
instead  heat  release  is  increased  compared  to  heat  losses,  for  example  by  lowering  the 
equivalence  ratio  from  3.2  to  2.8  at  a  constant  inlet  velocity,  then  the  combustion  zone  will 
migrate  upstream  to  a  stable  position  where  heat  recirculation  is  reduced.  With  further  reduction 
in  equivalence  ratio  the  flame  will  propagate  into  the  porous  flow  straighteners  and  flash  back. 

Figure  B.3  shows  the  stable  operating  conditions  for  premixed  heptane/air  reactants.  As  the 
velocity  increases,  flashback  occurs  at  leaner  conditions.  For  example,  flashback  occurs  at  4>  = 
2.9  at  an  inlet  velocity  of  50  cm/s.  At  u  =  200  cm/s,  where  heat  release  and  the  burning  rate  are 
greater,  flashback  occurs  at  the  leaner  condition  of  4>  =  2-7.  Similarly,  blowoff  occurs  at  (j)  =  4.0 
when  the  inlet  velocity  is  50  cm/s  and  <|>  =  3.6  when  the  inlet  velocity  is  200  cm/s.  Extinction 
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occurs  at  inlet  velocities  below  50  cm/s  due  to  excessive  heat  losses  relative  to  heat  release  rates. 
The  two  sets  of  test  conditions  at  which  exhaust  gas  composition  were  evaluated,  (|>  =  3.0  with 
inlet  velocity  varied  from  50-200  cm/s  and  u  =  125  cm/s  with  equivalence  ratio  varied  from  2.9- 
3.8,  are  highlighted  in  the  figure.  These  values  were  selected  because  the  reactor  operating  range 
is  relatively  large  and  equilibrium  predicts  that  the  conversion  of  heptane  to  hydrogen  is 
maximum  at  an  equivalence  ratio  of  2.9. 
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Figure  B.3:  Stable  and  unstable  operating  conditions,  with  highlighted  operating  points 

where  exhaust  gas  composition  was  analyzed 


Combustion  temperature  and  exhaust  gas  composition 

Although  equilibrium  calculations  should  not  necessarily  predict  the  products  of  the  reactor,  they 
do  form  a  standard  for  comparison  of  trends  as  well  as  degree  of  reaction  progression.  In  Table 
B.l,  equilibrium  values  are  determined  for  a  constant  pressure,  adiabatic  system.  Wet  and  dry 
equilibrium  values  are  tabulated  for  (j)  =  3.0  and  (j)  =  3.5,  where  dry  values  do  not  include  water  in 
the  equilibrium  composition  and  will  be  used  for  comparison  against  experimental  results. 
Oxygen,  heptane,  and  most  other  hydrocarbons  including  C2H2,  C2H6  and  C3H8  are  not  present 
in  significant  amounts  in  the  equilibrium  composition.  Methane  is  the  only  unbumed 
hydrocarbon  present  in  non-negligible  amounts,  and  its  concentration  increases  with  increasing 
equivalence  ratio.  Another  compositional  change  that  correlates  to  an  increasing  (|)  is  the 
appearance  of  solid  carbon.  The  adiabatic  equilibrium  temperature  is  predicted  to  decrease  from 
808.9  °C  to  774.6  °C  as  equivalence  ratio  increases  from  (j)  =  3.0  to  (j)  =  3.5. 
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3.0 

3.0 

3.5 

3.5 

WET 

DRY 

WET 

DRY 

T(°C) 

808.9 

808.9 

774.6 

774.6 

H2  [%] 

26.61 

26.82 

28.00 

28.33 

H20 

0.75 

— 

1.15 

— 

CO 

23.46 

23.64 

20.51 

20.75 

CH4 

0.29 

0.29 

0.49 

0.49 

C2H2 

<lE-5 

<lE-5 

<lE-5 

<lE-5 

C2H6 

<lE-5 

<lE-5 

<lE-5 

<lE-5 

C3H8 

<lE-5 

<lE-5 

<lE-5 

<lE-5 

N2 

48.17 

48.54 

44.50 

45.01 

C(s) 

0.00 

0.00 

4.35 

4.40 

Table  B.l:  Equilibrium  values  of  adiabatic  equilibrium  temperature  and  product 
composition  for  premixed  heptane/air  reactant  mixtures  of  <)>  =  3.0  and  <j)  =  3.5 


Figures  B.4-B.5  and  B.6-B.7  show  experimental  measurements  of  exhaust  gas  composition  and 
peak  reactor  wall  temperatures  as  functions  of  equivalence  ratio  and  inlet  velocity,  respectively. 
A  key  feature  of  heat  recirculating  reactors  is  their  ability  to  produce  superadiabatic 
temperatures,  which  are  necessary  to  increase  reaction  rates  and  achieve  high  conversion 
efficiency  [23,26].  Figures  B.4  and  B.6  illustrate  that  the  peak  wall  temperature  for  all  cases  is 
superadiabatic.  Temperatures  follow  the  trend  of  equilibrium  and  exceed  the  adiabatic 
equilibrium  temperature  by  190-460  °C  in  constant  equivalence  ratio  tests,  and  290-380  °C  in 
constant  inlet  velocity  tests. 


Figure  B.4  also  shows  the  impact  of  varying  equivalence  ratio  on  major  exhaust  components  H2 
and  CO.  Experimental  measurements  of  CO  concentration  show  little  change  over  the  range  of 
tested  equivalence  ratios,  and  results  are  within  the  average  uncertainty.  In  contrast  to 
equilibrium  values,  experimental  values  of  H2  concentrations  decrease  with  increasing 
equivalence  ratio.  Despite  peak  reactor  temperatures  in  excess  of  theoretical  adiabatic 
temperatures,  both  H2  and  CO  concentrations  are  below  equilibrium  values.  This  is  attributed  to 
a  significant  amount  of  unconverted  hydrocarbons  and  is  indicative  that  the  reactions  did  not 
proceed  to  equilibrium  within  the  residence  time  of  the  reactor.  Such  results  have  been  suggested 
in  calculations  done  previously  by  Al-Hamamre  et  al.  [6]  for  methane,  in  which  an  initial  sharp 
increase  in  product  concentration  is  observed  followed  by  a  more  gradual  increase  towards 
equilibrium  over  time  scales  that  are  significantly  longer  than  those  that  exist  in  the  present 
experiments. 


Figure  B.5  shows  experimental  values  of  intermediate  hydrocarbons  in  the  reactor.  Total 
unbumed  hydrocarbons  (UHC)  increase  monotonically  from  5.2%  to  9.8%  with  increasing 
equivalence  ratio.  Methane  and  ethylene,  with  H/C  ratios  of  4  and  2  respectively,  both  increase 
with  equivalence  ratio  while  acetylene,  with  an  H/C  ratio  of  1,  decreases.  This  shows  that 
hydrogen  is  increasingly  bound  in  hydrocarbons  rather  than  in  diatomic  hydrogen.  Ethane, 
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propane  and  /-butane  concentrations  increase  with  increasing  equivalence  ratio,  but  are  at  all 
points  present  in  amounts  less  than  0.2%  and  are  not  shown.  At  high  equivalence  ratios,  the 
results  show  carbon  monoxide  near  equilibrium  and  measured  UHC  in  excess  of  equilibrium. 
The  amount  of  carbon  in  exhaust  species  that  are  not  measured  by  GC  can  be  estimated  using  a 
carbon  balance.  Based  on  the  species  that  are  measured  by  GC,  with  carbon  numbers  up  to  C7,  it 
is  expected  that  the  unmeasured  species  are  larger  hydrocarbons  and  possibly  soot.  At  the  highest 
equivalence  ratio  of  3.8,  the  experimental  carbon  balance  shows  that  -11%  of  incoming  carbon 
is  contained  in  species  not  measured  by  the  GC.  In  contrast,  equilibrium  predicts  that  -25%  of 
incoming  carbon  goes  to  form  solid  carbon  product.  Therefore,  in  experiments  there  are  higher 
values  of  measured  UHC  than  equilibrium  predicts  but  lower  values  of  higher  hydrocarbons  and 
soot. 

Combustion  chemistry  literature  provides  insight  into  oxidation  of  rich  hydrocarbon  mixtures 
and  aids  in  interpretation  of  these  trends.  //-Heptane  oxidation  is  initiated  via  H  abstraction  by  H, 
O  or  OH  radicals.  The  resultant  heptyl  radical  is  then  susceptible  to  breakdown  by  //-scission  to 
produce  ethylene  and  pentyl  radical.  The  latter  undergoes  repeated  further  //-scission  to  produce 
ethylene  and  a  methyl  radical.  Ethylene  may  react  with  an  O  radical  to  form  a  methyl  radical  and 
HCO,  or  in  fuel  rich  mixtures  ethylene  undergoes  further  H  abstraction  by  H  and  OH  radicals  to 
produce  a  vinyl  radical,  which  quickly  reacts  to  form  acetylene.  Therefore  acetylene  appears  as  a 
product  further  along  in  the  progression  of  the  complex  hydrocarbon  breakdown  process. 
Methane  can  form  through  multiple  pathways,  both  as  an  early  and  a  late  breakdown  product 
[27].  Constant  inlet  velocity  results  show  via  the  correlation  of  methane  with  ethylene  trends  that 
its  formation  is  not  predominantly  due  to  further  reduction  of  acetylene. 

Although  hydrocarbons  larger  than  C7  were  not  measured,  soot  formation  was  observed  at  high 
equivalence  ratios  and  inlet  velocities.  The  data  are  too  limited  to  draw  conclusions  about  the 
pathway  of  polycyclic  aromatic  hydrocarbons  (PAH)  formation,  but  it  is  interesting  to  note  that 
acetylene,  which  is  present  at  all  operating  points  in  excess  of  equilibrium  concentrations, 
participates  in  the  chemistry  of  formation  of  PAH,  a  precursor  to  soot  [27]. 
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Figure  B.4:  Exhaust  gas  and  equilibrium  concentrations  of  H2  and  CO,  and  peak  wall  and 
adiabatic  equilibrium  temperatures  for  u  -  125  cm/s  and  varying  <|) 


Figure  B.5:  Exhaust  gas  concentrations  of  CH4,  C2H2,  C2H4  and  total  UHC,  and 
equilibrium  concentrations  of  CH4  (equilibrium  concentrations  of  C2H2  and  C2H4  are 
negligible  and  not  shown)  for  u  -  125  cm/s  and  varying  <|> 

Figure  B.6  shows  the  impact  of  inlet  velocity  on  peak  reactor  wall  temperature  and  on  the 
concentrations  of  H2  and  CO.  Equilibrium  values  are  shown  for  comparison.  In  experiments 
with  increasing  velocity,  the  volumetric  heat  release  rate  increases  resulting  in  higher  peak 
temperatures.  The  H2  concentration  is  influenced  by  two  competing  effects  in  this  system. 
Higher  temperatures  increase  the  reaction  rate  and  result  in  faster  fuel  conversion  [18,23], 
however,  increased  velocities  imply  shorter  residence  times  which  may  counteract  the  influence 
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of  the  higher  temperatures  [19].  These  competing  effects  are  reflected  in  experimental  values  of 
H2  concentrations,  which  increase  initially  at  velocities  up  to  80  cm/s  and  then  remain  nearly 
constant.  CO  concentrations  show  minimal  change  with  inlet  velocity.  The  increase  in  H2  and 
UHC  (Fig.  B.7)  with  no  corresponding  increase  in  carbon  monoxide  at  lower  inlet  velocities 
suggests  that  larger  intermediate  hydrocarbons  not  measured  by  the  GC  may  be  present  at  the 
lowest  u.  Molecular  hydrogen  is  released  as  those  species  are  further  broken  down,  while  carbon 
remains  bound  in  the  smaller  intermediate  hydrocarbon  species. 


Figure  B.7  shows  exhaust  gas  concentrations  of  CH4,  C2H2  and  C2H4,  as  well  as  total  UHC,  for 
(j)  =  3.0  and  varying  u.  All  other  measured  hydrocarbons  are  present  in  levels  less  than  0.1%. 
Total  unburned  hydrocarbons  remain  constant  within  the  margin  of  uncertainty  at  an  average 
value  of  5.5%  up  to  an  inlet  velocity  of  155  cm/s,  above  which  the  total  decreases  monotonically 
to  4.4%  at  200  cm/s.  Acetylene,  which  is  a  species  formed  late  in  hydrocarbon  breakdown, 
initially  increases  in  concentration  with  increasing  inlet  velocity  up  to  80  cm/s  and  subsequently 
decreases  at  higher  inlet  velocities.  Methane  and  ethylene  increase  with  inlet  velocity  up  to  u  = 
155  cm/s,  above  which  concentrations  decrease  with  further  increase  in  velocity.  Ethylene  is  less 
reduced  than  acetylene  and  forms  earlier  in  hydrocarbon  breakdown.  The  existence  of  a  peak 
acetylene  value  followed  by  ethylene  and  methane  maxima  at  higher  inlet  velocities  suggests  that 
fuel  breakdown  is  diminished  by  decreased  residence  time.  At  the  highest  inlet  velocities, 
ethylene  and  methane  concentrations  as  well  as  total  measured  unburned  hydrocarbons  decrease. 
In  conjunction  with  no  increase  in  H2  levels  and  visual  observation  of  significant  soot  formation, 
the  results  suggest  that  fuel  breakdown  is  promoted  by  high  temperatures  but  intermediate 
species  may  form  larger  soot  precursors  instead  of  further  reduction  at  the  highest  inlet 
velocities. 
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Figure  B.6:  Exhaust  gas  and  equilibrium  concentrations  of  H2  and  CO,  and  peak  wall  and 
adiabatic  equilibrium  temperatures  for  <|)  =  3.0  and  varying  u 
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Figure  B.7:  Exhaust  gas  concentrations  of  CH4,  C2H2,  C2H4  and  UHC,  and  equilibrium 
concentrations  of  CH4  (equilibrium  concentrations  of  C2H2  and  C2H4  are  negligible  and 

not  shown)  for  <|)  =  3.0  and  varying  u 


Fuels  comparisons 

The  conversions  of  methane  [18]  and  propane  [19]  to  syngas  were  investigated  in  earlier  studies 
using  the  counterflow  reactor  design.  We  compared  results  for  heptane  to  those  from  the 
previous  publications  using  three  metrics:  species  conversion  efficiency,  and  species  and  total 
energy  conversion  efficiencies. 

Species  conversion  efficiency,  also  called  species  yield,  compares  actual  species  production  to 
the  theoretical  maximum  based  on  the  amount  of  that  species’  components  bound  in  the  reactant 
fuel  stream.  Equation  B.l  gives  H2  conversion  efficiency,  which  is  the  number  of  moles 
generated  compared  to  the  number  of  moles  of  hydrogen  bound  in  the  heptane  feedstock. 


t1h% 


AC 


AC 


N  8  N 

1  v  H2  .max  oly  C7Hl6 


(B.l) 


Species  energy  conversion  efficiency  represents  the  amount  of  chemical  energy  bound  in  a 
product  species  compared  to  the  total  chemical  energy  in  the  reactant  fuel.  Chemical  energy  is 
represented  by  the  lower  heating  value  (LHV)  in  kJ/kmol  at  standard  conditions.  Equation  B.2 
gives  H2  energy  conversion  efficiency. 


46 


(B.2) 


V  En  ergy,l 


nhlhvh 2 

Nch  inLHVCH 


Total  energy  conversion  efficiency  compares  the  chemical  energy  bound  in  all  of  the  energy 
containing  product  species  to  the  initial  chemical  energy  of  the  reactant  fuel.  The  difference 
between  the  two  is  the  sum  of  energy  lost  to  surroundings  and  sensible  enthalpy  of  the  product 
gas.  Equation  B.3  gives  total  energy  conversion  efficiency. 
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Figures  B.8-B.10  show  experimental  and  equilibrium  results  of  methane,  propane  and  heptane 
reforming  in  terms  of  the  above  three  metrics  for  experiments  in  which  inlet  velocity  was  held 
constant  at  125  cm/s  and  equivalence  ratio  was  varied.  Experiments  were  conducted  over 
different  equivalence  ratio  ranges  because  of  differences  in  stability  when  operating  the  reactor 
with  the  three  fuels.  As  previously  discussed,  the  combustion  zone  is  stable  within  the  reactor 
when  the  inlet  velocity  equals  the  burning  rate,  which  is  a  function  of  temperature.  Variation  in 
temperature  between  the  three  studies  is  due  to  differences  in  the  adiabatic  equilibrium 
temperatures  of  different  reactant  species,  equivalence  ratios,  and  levels  of  preheat.  Heptane  and 
propane  have  higher  adiabatic  equilibrium  temperatures  than  methane.  In  addition,  heptane  is 
preheated  to  maintain  the  vapor  phase;  therefore  preheated  heptane  has  a  higher  adiabatic 
equilibrium  temperature  than  propane.  For  example,  the  adiabatic  equilibrium  temperatures  of 
propane  and  heptane  at  25  °C  for  (|>  =  2.9  are  850  °C  and  847  °C,  whereas  the  equilibrium 
temperature  of  heptane  increases  to  934  °C  when  the  reactants  are  preheated  to  150  °C.  Due  to 
these  differences,  an  equivalence  ratio  and  velocity  combination  that  is  in  the  stable  range  but 
near  the  flashback  limit  for  methane  will  be  outside  the  stable  range  for  propane.  The  same  is 
true  for  propane  and  heptane.  A  comparison  of  burner  maps  for  the  three  studies  [18,19]  further 
illustrates  these  stability  differences. 

Methane  is  the  only  one  of  the  three  fuels  that  shows  a  maximum  value  in  the  experimental 
results  of  hydrogen  species  (Fig.  B.8)  and  energy  (Fig.  B.9)  conversion  efficiencies  under  the 
tested  conditions.  Methane,  propane  and  heptane  efficiencies  decrease  with  increasing 
equivalence  ratio  above  (j)  =  2.2,  and  experimental  conversion  efficiencies  are  significantly  less 
than  predicted  by  equilibrium  at  these  points,  despite  a  peak  hydrogen  conversion  efficiency  for 
heptane  predicted  at  t()  =  2.9  by  equilibrium.  The  occurrence  of  peak  conversion  efficiency  at 
equivalence  ratios  below  those  predicted  by  equilibrium  is  consistent  with  the  results  of  Al- 
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Hamamre  et  al.  [6]  for  methane,  in  which  incomplete  conversion  occurs  on  counterflow  reactor 
timescales.  This  suggests  that  propane  and  heptane  might  also  produce  peak  H2  at  a  leaner 
operating  point  if  such  stable  operation  were  possible  but,  as  discussed  above,  leaner  operation  is 
limited  by  flashback  and  therefore  these  points  could  not  be  tested.  The  results  also  indicate  that 
the  current  reactor  design  does  not  permit  complete  fuel  conversion  at  high  equivalence  ratios. 


Equivalence  Ratio 


Figure  B.8:  Experimental  and  equilibrium  values  of  H2  conversion  efficiencies  of  methane 
[21],  propane  [22]  and  heptane  reforming  with  constant  u  =  125  cm/s  and  varying 

equivalence  ratio 
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Figure  B.9:  Experimental  and  equilibrium  H2  energy  conversion  efficiencies  of  methane 
[21],  propane  [22]  and  heptane  reforming  with  constant  u  =  125  cm/s  and  varying 

equivalence  ratio 


Experimental  values  of  total  energy  conversion  efficiencies  show  trends  and  values  similar  to 
equilibrium  (Fig.  B.10)  for  the  three  data  sets.  In  each  case  unbumed  hydrocarbon  concentrations 
increase  with  increasing  equivalence  ratio  [18,19],  and  the  chemical  energy  bound  in  those 
hydrocarbons  contributes  significantly  to  the  total.  Peak  reactor  temperatures  decrease  with 
increasing  equivalence  ratio  for  all  three  fuels  [18,19].  This  lessens  convective  and  radiative 
losses  to  the  environment  as  well  as  energy  required  to  increase  product  gas  enthalpy. 
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Figure  B.10:  Experimental  and  equilibrium  total  energy  conversion  efficiencies  of  methane 
[21],  propane  [22]  and  heptane  reforming  with  constant  u  =  125  cm/s  and  varying 

equivalence  ratio 

Figure  B.ll  shows  the  results  of  methane,  propane  and  heptane  reforming  in  terms  of  species 
and  energy  conversion  efficiencies  for  experiments  in  which  inlet  velocity  was  varied  from  50- 
200  cm/s  and  equivalence  ratio  was  held  constant  at  <])  =  2.2  for  methane,  (()  =  2.4  for  propane  and 
cj)  =  3.0  for  heptane.  All  conversion  efficiencies,  with  the  exception  of  methane  total  energy 
conversion  efficiency,  increase  initially  with  inlet  velocity  and  show  little  change  at  higher  u. 
The  initial  increase  supports  the  discussion  above  for  heptane  that  conversion  improves  with 
increased  energy  release  rates  and  resultant  higher  peak  temperatures.  The  marginal  change  over 
the  remainder  of  the  tested  inlet  velocity  range  suggests  that  product  composition  is  affected  by 
competing  effects  of  increased  temperature,  which  promotes  conversion,  and  reduced  residence 
time,  which  inhibits  conversion.  The  decrease  in  total  energy  conversion  efficiency  of  heptane  at 
the  highest  u  is  attributed  to  the  formation  of  soot,  which  consists  of  high  chemical  energy 
species  that  are  not  accounted  for  by  GC  measurements. 

Heptane  conversion  has  also  been  studied  in  filtration  reactors  consisting  of  a  column  of  porous 
media  which  showed  a  broader  operating  range  of  1.4-3. 8  in  comparison  to  the  stable  range  of 
2. 8-3. 9  in  the  counterflow  reactor  [23];  however  both  reactors  were  tested  at  or  near  an 
equivalence  ratio  of  2.9,  where  peak  hydrogen  conversion  efficiency  is  estimated  from 
equilibrium  calculations  to  occur.  The  tested  range  in  the  current  study  also  partially  overlaps  the 
equivalence  ratio  range  of  2. 5-3. 5  that  was  estimated  by  Dixon  et  al.  [23]  to  be  favorable 
conditions  for  heptane  reforming  in  a  filtration  reactor.  Pastore  and  Mastorakos  [24]  reported 
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that  the  peak  conversion  efficiencies  for  heptane  occurred  at  equivalence  ratios  of  1.8-  2.5 
depending  on  the  reactor.  These  differing  values  of  the  optimum  equivalence  ratio  indicate  that 
conversion  is  sensitive  to  the  details  of  the  reactor  and  operating  conditions. 

Peak  H2  conversion  efficiencies  from  the  counterflow  reactor  were  less  than  those  obtained  in 
porous  media,  from  which  efficiencies  of  greater  than  80%  were  obtained.  Hydrocarbon  product 
concentrations  in  the  porous  reactor  were  significantly  lower  than  those  measured  in  the 
counterflow  reactor,  at  less  than  1%  compared  to  5.8%,  indicating  more  complete  fuel 
breakdown.  Peak  temperatures  in  porous  media  were  in  excess  of  1400  °C  at  (j)  =  3.0  and  u  =  60 
cm/s,  compared  to  1190  °C  at  (j)  =  3.0  and  u  =  125  cm/s  in  the  counterflow  reactor.  Residence 
times  for  the  porous  media  and  counterflow  reactors  were  approximately  the  same  at  these 
operating  conditions,  with  both  on  the  order  of  100  ms.  The  difference  in  conversion  efficiency 
is  attributed  to  higher  peak  temperatures  in  the  porous  media  reactor  [19,28]. 
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Figure  B.ll:  H2,  H2  energy  and  total  energy  conversion  efficiencies  of  methane  (<j)  =  2.2) 
[21],  propane  (<|)  =  2.4)  [22]  and  heptane  (<j)  =  3.0)  with  varying  inlet  velocity 

B.3  Conclusions 

Experimental  results  of  non-catalytic  reforming  of  heptane  to  produce  hydrogen-rich  syngas  in  a 
counterflow  reactor  were  presented.  A  key  advantage  of  this  reactor  is  that  the  reaction  zones  are 
stationary,  which  permits  continuous  operation  and  practical  applications.  Stable  operation  of 
the  reactor  with  premixed  heptane/air  reactants  was  established  for  equivalence  ratios  between 
2.8  and  3.9  and  inlet  velocities  between  50  cm/s  and  200  cm/s.  Peak  wall  temperatures  were 
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monitored  throughout  reactor  operation  and  were  at  all  points  in  excess  of  the  adiabatic 
equilibrium  temperature,  thus  excess-enthalpy  flames  were  achieved  through  internal  heat 
recirculation. 

Exhaust  gas  composition  was  measured  using  gas  chromatography  for  two  sets  of  tests: 
equivalence  ratio  was  varied  from  2. 9-3. 8  with  inlet  velocity  held  constant  at  125  cm/s,  and 
equivalence  ratio  was  held  constant  at  3.0  with  inlet  velocity  varied  from  50-200  cm/s.  H2 
concentration  decreased  with  increasing  equivalence  ratio  from  a  maximum  of  14.5%  to  a 
minimum  of  10.6%,  while  CO  concentration  remained  nearly  constant  around  17.0%.  Both  H2 
and  CO  concentrations  showed  little  dependence  overall  on  inlet  velocity,  with  the  former 
varying  from  12.8-14.4%  and  the  latter  varying  between  16.6-17.5%.  The  largest  effect  of  inlet 
velocity  was  found  at  the  lower  end  of  the  tested  range  where  rates  of  heat  release  and 
temperatures  are  lower. 

In  all  cases,  unburned  hydrocarbons  were  measured  in  excess  of  concentrations  predicted  by 
equilibrium  composition.  Total  unburned  hydrocarbon  levels  increased  with  increasing 
equivalence  ratio.  Soot  was  noted  with  significant  formation  observed  at  higher  equivalence 
ratios  and  inlet  velocities.  These  results  and  observations  strongly  suggested  that  the  extent  of 
fuel  breakdown  is  reduced  by  insufficient  reaction  rate  and  residence  time. 

The  results  of  this  study  showed  trends  similar  to  previously  published  studies  on  syngas 
production  from  methane  and  propane  using  the  counterflow  reactor  [18,19].  Each  fuel  had  H2 
species  and  energy  conversion  efficiencies  below  equilibrium  values,  but  total  conversion 
efficiencies  were  comparable  to  equilibrium.  In  addition,  energy-dense  unburned  hydrocarbon 
concentrations  were  in  excess  of  equilibrium  with  concentrations  an  order  of  magnitude  higher 
than  predicted  by  equilibrium  at  the  highest  tested  equivalence  ratios  in  this  study.  The 
comparison  between  the  different  fuels  indicates  that,  although  the  stable  operating  range  for 
heptane  does  include  the  theoretically  optimal  value  as  predicted  by  hydrogen  conversion 
efficiencies  calculated  from  equilibrium,  the  current  reactor  design  does  not  permit  complete  fuel 
conversion  at  high  equivalence  ratios.  In  addition,  the  formation  of  soot  at  some  conditions 
suggests  the  importance  of  operating  at  more  moderate  conditions  for  optimal  performance. 

This  study  demonstrated  the  fuel  flexibility  of  the  counterflow  reactor  by  reforming  a  liquid  fuel, 
and  results  supplemented  previous  reforming  studies  on  the  counterflow  reactor.  Peak  H2 
conversion  efficiency  was  significantly  lower  than  that  attained  via  filtration  combustion  in 
porous  media,  an  alternate  method  of  non-catalytic  reforming.  Differences  were  attributed  to 
higher  peak  temperatures  in  the  porous  media  reactor,  which  has  higher  specific  surface  area  for 
heat  transfer  as  compared  to  the  counterflow  reactor  under  investigation.  Furthermore,  although 
the  tested  equivalence  ratio  range  in  this  study  does  include  (j)  =  2.9,  where  peak  conversion 
efficiency  is  predicted  to  occur  by  equilibrium,  previous  studies  predict  peak  conversion 
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efficiencies  over  a  range  of  (j)  =  1.8  to  (|>  =  3.5.  The  results  of  hydrogen  conversion  in  this  study 
indicate  that  the  optimum  equivalence  ratio  might  be  outside  the  operating  range  of  the 
counterflow  reactor. 
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Appendix  C:  Critical  design  parameters  for  counterflow  reactor 

The  third  objective  of  this  ARO-funded  project  was  to  computationally  investigate  the  critical 
design  parameters  for  reforming  in  the  counterflow  reactor.  This  work  is  described  below,  and  is 
in  preparation  for  publication. 

The  counterflow  reactor  design  was  proposed  analytically  [1,2,3]  and  validated  in  experimental 
studies  on  the  reforming  of  methane  [4],  propane  [5]  and  heptane  [6]  to  produce  syngas.  A 
computational  model  of  the  counterflow  reactor  was  developed  to  study  the  reaction  zone 
characteristics  of  rich  methane  combustion,  in  which  the  reactor  channels  were  modeled  in  two 
dimensions,  and  the  detailed  GRI  2.11  kinetics  mechanism  was  used.  The  model  was  validated 
against  experimental  results,  where  the  impact  of  inlet  velocity  and  equivalence  ratio  operating 
conditions  on  combustion  zone  characteristics  was  successfully  modeled  [7]. 

This  study  utilized  the  previously  developed  computational  model  to  study  the  impact  of 
geometric  scaling  on  rich  methane  reforming  in  the  counterflow  reactor.  Scalability  is  an 
important  consideration  for  portable  power  applications,  as  well  as  for  large-scale,  stationary 
applications.  In  either  case,  reactor  scaling  may  be  achieved  by  changing  the  geometry  of  the 
individual  channels  or  by  altering  the  number  of  parallel  channels  [8].  This  study  focused  on 
geometry  at  the  individual  channel  scale. 

Several  studies  have  examined  the  effects  of  geometric  scaling  on  the  stability  and  performance 
of  meso-scale,  heat-recirculating  combustors  [9].  The  parameters  that  have  been  examined 
include  channel  wall  thickness  and  material  in  Swiss-roll  and  parallel-channel  reactors 
[10,11,12,13],  as  well  as  channel  length  [14]  and  channel  width  [12,15,16,17,18,19].  Kaisare  and 
Vlachos  [20]  numerically  studied  the  effect  of  reactor  length,  wall  thickness  and  reactor  opening 
size  on  flame  stability  in  micro-  and  meso-scale  channel  reactors. 

The  effect  of  scaling  the  counterflow  reactor  on  non-catalytic  fuel  reforming  was  the  focus  of 
this  computational  study,  where  reactor  channel  length  or  height  were  varied  from  the  original 
dimensions  examined  in  previous  experimental  [4-6,21]  and  analytical  investigations  [3,21].  The 
impact  of  scaling  on  methane  reforming  was  quantified  in  a  number  of  ways:  combustion 
stability,  efficiency  of  reactant  conversion  to  hydrogen,  efficiency  of  heat  recirculation  from  hot 
combustion  products  to  incoming  cold  reactants,  and  degree  of  superadiabicity,  defined  as  the 
ratio  of  peak  to  equilibrium  temperatures,  achieved  in  the  reactor.  Computational  investigation  of 
reactor  operating  range  and  hydrogen  conversion  provided  an  efficient  means  of  testing 
numerous  reactor  geometries.  Quantification  of  heat  recirculation  efficiency  and  degree  of 
superadiabicity  provided  insights  into  reactor  operation  that  are  not  readily  measured  during 
experimentation.  Heat  recirculation  efficiency  in  the  counterflow  reactor  was  also  compared  to 
that  achieved  in  a  porous  reactor  [22] . 
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C.l  Numerical  method 


Numerical  model 


The  data  presented  in  this  study  were  obtained  from  simulations  performed  using  the 
computational  fluid  dynamics  (CFD)  code  Fluent  14.5.  The  CFD  code  solves  conservation  of 
mass,  energy  and  species  equations  and  the  incompressible  Navier-Stokes  equations.  The 
conservation  of  mass  is  given  by 

+  V  ‘  (Pfv)  =  o  (C.l) 

and  the  conservation  of  momentum  is  given  by 

d-^p-  +  V  ■  (pfm)  =  -Vp  +  V  ■  f  +  F  (C.2) 

where  Pf  is  fluid  density,  v  is  the  velocity  vector,  and  p  is  the  static  pressure.  The  stress  tensor, 
r,  and  F  represent  volumetric  forces  including  viscous  forces.  The  fluid  is  approximated  by  the 
ideal  gas  law,  thus  fluid  density  is  given  by  Pf  =  p/(RT /Mw)  where  R  is  the  ideal  gas  constant, 
T  is  temperature  and  Mw  is  the  molecular  weight  of  the  gas. 

Conservation  of  energy  for  the  fluid  phase  is  given  by 

+  V  ■  (v(pfE  +  p))  =  V  ■  (/C/VT  -  Si  hJi)  +  ShJ  (C.3) 

where  E  —  h  —  p / Pf  +  v2  and  h  is  sensible  enthalpy,  defined  for  incompressible  fluid  as 
h  =  Si  Tj/ij  with  Y[  as  the  mass  fraction  of  species  i.  Thermal  conductivity  is  given  by  k,  J 
represents  Fickian  species  diffusive  flux,  and  Shj  is  a  volumetric  source  term  that  includes  the 
heat  of  chemical  reactions. 

Conservation  of  energy  for  the  solid  phase  is  given  by 

=  V  ■  (/csVT)  +  SKs  (C.4) 

where  Sh  s  represents  user-defined  volumetric  heat  sources. 

Conservation  of  species  is  given  by 
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d(PfYj) 


+  V-(pfvYi)  =  ~V-Ji  +  Ri 


(C.5) 


dt 


where  the  net  rate  of  production  of  species  i,  Ri,  is  the  sum  of  Arrhenius  reaction  rates  for  N 
reactions.  Combustion  chemistry  is  modeled  using  detailed  reaction  chemistry  mechanism  GRI 
2.1 1.  Further  details  about  the  numerical  model  can  be  found  in  [7]. 


Model  Geometry 

Figure  C.la  is  a  schematic  of  the  four-channel  reactor  used  in  previous  experimental 
investigations  [4-6,21].  Figure  C.lb  is  a  schematic  of  the  computational  domain  of  the  original 
reactor  geometry,  with  dimensions  of  the  counterflow  reactor  used  in  previous  experimental 
reforming  and  thermal  oxidation  studies.  The  domain  is  reduced  to  one  half-channel,  which 
lessens  the  computational  load  significantly,  through  the  use  of  symmetry  and  a  user-defined 
function  (UDF).  Symmetry  in  the  y-direction  is  imposed  at  the  channel  center  line,  and 
temperature  symmetry  is  imposed  by  a  UDF  at  the  wall  center  line.  External  losses  are  accounted 
for  by  radiative  wall  ends  and  convective  losses  from  the  wall.  Therefore,  the  computational 
model  approximates  the  reactor  as  an  infinite  number  of  parallel  channels. 
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Figure  C.l:  (a)  A  schematic  of  the  four-channel  reactor  and  (b)  the  computational  domain 
of  the  original  reactor  geometry  (VR  1.0) 


The  reactors  modeled  in  this  study  are  scaled  relative  to  the  volume  of  the  original  reactor  (122 
cm3)  by  a  volume  ratio,  defined  as 


Volume  Ratio(VR )  = 


Volume  of  scaled  reactor 
Volume  of  original  reactor 


(C.6) 
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where  the  scaled  reactor  volume  is  achieved  by  varying  either  channel  height  or  main  reactor 
section  length,  while  maintaining  the  other  dimension  at  the  reference,  original  reactor  value.  In 
addition,  geometry  is  limited  by  maintaining  the  original  reactor  head  dimensions,  except  as 
channel  height  is  varied,  and  channel  wall  thickness  is  maintained  at  the  original  value  of  1  mm. 
Volume  ratios  of  0.75,  1.0  and  1.25  are  investigated,  where  volume  ratio  1.0  (VR  1.0) 
corresponds  to  the  original  reactor  design  and  dimensions.  Table  C.l  shows  the  dimensions  of 
each  scaled  reactor. 

Table  C.l:  Counterflow  reactor  dimensions  for  volume  ratios  (VR)  0.75,  1.0  and  1.25, 
obtained  by  varying  channel  height  or  channel  length  with  the  other  held  constant  at  its 
reference  value 


Volume  Ratio 
(VR) 

Reactor  ID 

Channel  height 
(mm) 

Total  channel  length 
(mm) 

1.0 

VR  1.0 

4.0 

173.1 

0.75 

VR(h)  0.75 

2.7 

173.1 

1.25 

VR(h)  1.25 

5.3 

173.1 

0.75 

VR(F)  0.75 

4.0 

129.8 

1.25 

VR(F)  1.25 

4.0 

216.4 

Solution  Approach 

Scaled  reactors  were  modeled  in  ANSYS  ICEM  CFD.  The  initial  model  grids  had  axial  and 
transverse  spacing  equal  to  that  used  in  Schoegl  et  al  [7].  A  base  case  of  each  reactor  model  was 
run  for  20k  time  steps  of  size  2e-5  seconds  at  inlet  conditions  of  4>=2.2  and  u=  1 25  cm/s  in  order 
to  achieve  ignition  and  initial  stabilization  of  the  flame  front  in  the  main  reactor  section  of  the 
channel.  The  solution  of  the  base  case  was  then  used  as  an  initial  solution  for  each  tested 
equivalence  ratio  and  inlet  velocity  operating  point.  Solutions  were  continued  until  near-steady 
state  was  achieved,  where  steady  state  was  determined  to  occur  when  the  flame  front  velocity 
within  the  main  section  of  the  reactor  model  was  less  than  0.5  mm/s.  Flame  front  location  was 
determined  by  the  location  of  peak  H  radical  concentration. 

After  near-steady  state  was  achieved,  the  mesh  for  each  case  was  refined  in  a  three  phase 
process.  Each  mesh  was  locally  refined  in  regions  of  high  concentrations  of  specific  chemical 
species,  following  [7].  This  process  effectively  increased  the  mesh  resolution  up  to  64-fold 
around  the  combustion  zone  by  partitioning  quadrilateral  elements  into  four  sub-elements. 
Following  each  stage  of  mesh  refinement,  simulations  were  continued  for  a  total  of  5k  additional 
time  steps  of  size  le-5  seconds.  After  the  third  refinement,  the  spatial  discretization  solution 
technique  was  changed  from  first-order  upwind  to  second-order  upwind,  and  the  final  result  was 
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obtained  after  an  additional  3k  time  steps  of  size  le-5  s.  The  Semi-Implicit  Method  for  Pressure 
Linked  Equations  (SIMPLE)  algorithm  was  used  for  pressure-velocity  coupling  during  the 
transient  calculations. 

C.2  Results  and  Discussion 

Results  are  presented  below  for  each  of  the  five  scaled  reactors  over  the  range  of  tested  inlet 
velocities  and  equivalence  ratios.  Data  include  operating  conditions  where  steady  state 
combustion  was  achieved  in  the  main  section  of  the  reactor,  as  well  as  heat  recirculation 
efficiency,  gas  and  wall  temperature  ratios,  and  hydrogen  conversion  efficiencies  at  those  steady 
state  points. 

Steady  State  Combustion 

Each  of  the  five  reactor  geometries  were  tested  over  a  range  of  inlet  velocities  and  equivalence 
ratios  to  determine  which  operating  points  produce  a  steady  state  flame  front  solution  in  the  main 
section  of  the  reactor  channel,  which  corresponds  to  stable  reactor  operating  points  in 
experimental  studies.  Inlet  velocity  was  varied  from  50-200  cm/s  in  increments  of  25  cm/s  at 
(j)=2.2,  and  equivalence  ratio  was  varied  from  c|)=l .4-2.6  in  increments  of  0.2  at  «=125  cm/s. 
Figure  C.2  shows  the  operating  limits  within  these  ranges  for  the  three  reactor  geometries  that 
were  scaled  by  channel  height.  Beyond  these  limits,  the  steady  state  flame  front  solution  is  either 
upstream  of  the  main  reactor  section  at  the  porous  media  flow  straightener,  which  corresponds  to 
flashback  in  experimental  investigations,  or  the  flame  front  propagated  out  of  the  channel  exit, 
which  corresponds  to  blow-off  in  experimentation. 

In  the  case  of  equivalence  ratio  variation,  the  computations  indicate  flashback  and  blow-off  at 
equivalence  ratios  below  and  above  the  lowest  and  highest  points,  respectively,  for  which  steady 
state  solutions  are  obtained  in  the  main  reactor  section.  When  the  volume  ratio  is  increased  by 
increasing  channel  height,  the  equivalence  ratio  at  which  flashback  occurs  decreases.  For  VR(h) 
0.75,  flashback  occurred  at  c|>=1.8.  For  VR  1.0  and  VR(h)  1.25,  flashback  occurred  at  c|)=1.6  and 
(j)=1.4,  respectively.  The  results  show  a  shift  in  stable  equivalence  ratio  range  towards 
stoichiometric  as  channel  height  is  increased.  Nearly  all  of  the  tested  equivalence  ratios  are 
above  the  flammability  limit  of  (|)=1.67  for  methane  at  inlet  conditions,  and  preheating  of 
reactants  is  required  to  achieve  the  observed  extension  of  flammability  limits  [23].  Therefore,  the 
operating  range  results  suggest  that  preheating  becomes  more  effective  with  decreasing  channel 
height.  As  channel  height  increases,  axial  convection  begins  to  dominate  over  transverse  heat 
and  mass  diffusion,  where  the  latter  are  required  for  conversion  of  fuel  at  the  channel  centerline. 

In  the  case  of  inlet  velocity  variation,  the  two  lower  volume  ratio  reactors  (VR(h)  0.75  and  VR 
1.0)  produced  steady  state  flame  fronts  in  the  main  reactor  section  over  the  entire  tested  inlet 
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velocity  range,  while  VR(h)  1.25  became  unsteady  due  to  blow-off  at  w=150  cm/s.  Therefore  the 
largest  volume  ratio  reactor,  VR(h)  1.25,  is  the  only  one  of  the  three  reactors  scaled  by  height 
that  is  limited  over  the  tested  inlet  velocity  range.  Blow-off  occurred  at  «=150  cm/s  because  the 
reactants  were  insufficiently  preheated  to  raise  the  burning  velocity  to  equal  the  inlet  velocity, 
which  is  the  criterion  for  steady  state  operation  in  the  counterflow  reactor.  Similar  results  were 
observed  in  Oh  et  al  [10]  and  Norton  et  al  [12],  where  larger  channel  widths  resulted  in  lower 
reaction  rates  due  to  slow  heat  transfer  from  the  wall  where  ignition  starts  [7]  towards  the 
reactants  at  the  channel  centerline.  In  contrast,  smaller  separation  distances  between  channel 
walls  produce  more  localized  reaction  zones  with  steeper  temperature  and  composition  gradients. 
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Figure  C.2:  Steady  state  operating  limits  for  counterflow  reactors  scaled  to  smaller  and 
larger  volume  ratios  by  varying  channel  height  (VR(h)  0.75  and  VR(h)  1.25)  and  original 
reactor  geometry  (VR  1.0)  over  the  tested  range  of  inlet  velocities,  m=50-200  cm/s  at  <|>=2.2, 
and  the  tested  range  of  equivalence  ratios,  <j)=1.4-2.6  at  «=125  cm/s. 

Figure  C.3  shows  the  steady  state  operating  limits  within  the  ranges  of  tested  inlet  velocities  and 
equivalence  ratios  for  the  three  reactor  geometries  scaled  by  channel  length.  In  the  case  of 
equivalence  ratio  variation,  flashback  occurred  at  decreasing  equivalence  ratios  as  volume  ratio 
increased:  (j)=1.8  for  VR(L)  0.75,  (j)=1.6  for  VR  1.0  and  (j)=1.4  for  VR(L)  1.25.  Blow-off  occurred 
at  (j)=2.4  for  the  reactor  with  the  shortest  length,  VR(L)  0.75,  whereas  steady  state  operation  was 
achieved  at  the  upper  end  of  the  tested  equivalence  ratio  range  for  the  longer  reactors  VR  1 .0  and 
VR(L)  1.25.  In  the  case  of  inlet  velocity  variation,  the  two  longer  reactors,  VR(L)  1.25  and  VR 
1.0,  produced  steady  state  solutions  throughout  the  tested  range  of  50-200  cm/s.  The  lowest 
volume  ratio  reactor  with  shortest  channel  length,  VR(L)  0.75,  had  a  steady  state  range  of  100- 
150  cm/s,  due  to  flashback  at  75  cm/s  and  blow-off  at  175  cm/s. 

The  results  suggest  that  steady  state  range  increases  with  channel  length,  however  VR(L)  0.75 
and  VR  1.0  produced  steady  state  results  throughout  the  tested  inlet  velocity  range  and  therefore 
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the  extents  of  inlet  velocity  range  were  not  tested.  The  observed  steady  state  solutions  are 
consistent  with  experimental  observations  of  the  stable  operating  range.  It  is  likely  that  there  is 
an  optimal  channel  length  for  maximum  stability  range  due  to  the  competing  effects  of  heat  loss 
and  stabilization  space;  as  channel  length  increases,  surface  area  and  thus  heat  loss  also  increase, 
however  there  is  also  increased  axial  space  for  the  reaction  front  to  stabilize.  Since  there  is  no 
predetermined  combustion  zone  location  in  the  main  section  of  the  reactor,  the  flame  front  is  free 
to  stabilize  at  the  point  where  the  burning  velocity  equals  the  local  inlet  velocity.  For  example, 
increasing  the  equivalence  ratio  from  (j)=2.0  to  <J)=2.2  at  a  fixed  inlet  velocity  results  in  lower  heat 
release.  The  flame  front  will  then  stabilize  further  downstream,  where  the  heat  recirculation  is 
more  significant  and  the  reactants  can  be  further  preheated  prior  to  combustion.  Likewise,  the 
flame  front  location  will  adjust  to  changes  in  inlet  velocity  at  a  fixed  equivalence  ratio;  the  flame 
front  will  shift  downstream  in  response  to  an  increase  in  inlet  velocity,  where  additional  heat 
recirculation  will  provide  sufficient  preheating  for  the  burning  velocity  to  equal  the  inlet 
velocity.  The  reactor  scaled  down  by  channel  length,  VR(L)  0.75,  has  the  most  limited  axial 
distance  within  the  main  section  for  the  reaction  front  to  adjust  to  changes  in  operating 
conditions  before  blow-off  or  flashback  will  occur,  thus  its  stable  range  is  significantly  narrower 
than  longer  reactors  VR  1.0  and  VR(L)  1.25. 
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Figure  C.3:  Steady  state  operating  limits  for  counterflow  reactors  scaled  to  smaller  and 
larger  volume  ratios  by  varying  channel  length,  VR(L)  0.75  and  VR(L)  1.25,  and  original 
reactor  geometry,  VR  1.0,  over  the  tested  range  of  inlet  velocities,  n=50-200  cm/s  at  <j)=2.2, 
and  the  tested  range  of  equivalence  ratios,  <|)=1.4-2.6  at  n=125  cm/s. 

Heat  Recirculation  Efficiency 

The  counterflow  reactor  transfers  heat  by  conduction  and  radiation  upstream  from  hot 
combustion  products  in  the  post-flame  region  to  cold  reactants  in  the  preheat  region.  Heat  is  also 
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AVR  1.0 
XVR(L)  1.25 
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transferred  through  the  channel  walls  to  the  unburned  reactants  counterflowing  in  the  adjacent 
channels,  thereby  preheating  the  reactants.  Figure  C.4  shows  temperature  profiles  for  gas  at  the 
reactor  channel  centerline  and  the  channel  wall,  obtained  for  the  original  reactor  geometry  (VR 
1.0)  at  operating  conditions  of  4>=2.2  and  u=  125  cm/s.  The  preheat  zone  is  identified  as  the 
region  between  the  channel  inlet  and  the  point  where  wall  and  gas  temperatures  are  equal. 
Upstream  of  this  point,  the  wall  temperature  exceeds  the  gas  temperature  and  heat  is  transferred 
to  the  gas.  The  two  sections  of  porous  media  flow  straighteners  (PM)  are  identified,  and  the 
impact  of  these  conductive  inserts  on  increasing  the  gas  temperature  can  be  seen.  Downstream  of 
the  preheat  zone,  the  reaction  zone  occurs  and  gas  temperatures  exceed  wall  temperatures.  Heat 
is  then  transferred  from  the  gas  to  the  wall,  where  it  conducts  and  radiates  upstream  and  preheats 
reactants  in  adjacent  channels. 
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Figure  C.4:  Gas  (Tgas)  and  wall  (Twall)  temperature  profiles  along  the  axial  length  of  a 
counterflow  reactor  channel,  obtained  from  the  original  reactor  geometry  (VR  1.0)  with 
inlet  conditions  of  <(>=2.2  and  u= 125  cm/s.  The  preheat  zone  is  identified,  as  well  as  the 
locations  of  the  two  sections  of  porous  media  flow  straighteners  (PM),  the  high  thermal 
conductivity  of  which  contributes  to  preheating  the  reactants. 

The  amount  of  heat  transferred  to  the  reactants  in  the  preheat  zone  is  obtained  by  integrating  the 
heat  flux  into  the  channel  upstream  of  the  reaction  zone,  from  the  inlet  to  the  point  where  the  gas 
and  solid  temperatures  are  equal  and  the  flux  equals  zero.  Heat  recirculation  efficiency  is  a 
measure  of  this  heat  transfer  and  is  defined  as 

Heat  recirculation  efficiency  — -  (C.7) 

1  1  '  Firing  rate 
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where  firing  rate  is  calculated  based  upon  the  lower  heating  value  of  methane. 


Figure  C.5  shows  the  calculated  heat  recirculation  efficiency  for  each  of  the  five  scaled  reactors 
as  a  function  of  equivalence  ratio.  All  reactors  show  a  decrease  in  heat  recirculation  efficiency 
with  increasing  equivalence  ratio.  A  significant  increase  in  efficiency  is  found  with  decreasing 
channel  height,  where  the  smallest  channel  height  VR(h)  0.75  shows  approximately  30%  and 
26%  heat  recirculation  efficiency  at  (j)=2.0  and  (j)=2.2,  respectively,  versus  19%  and  18% 
efficiency  produced  by  the  largest  channel  height  VR(h)  1.25.  The  increase  in  heat  recirculation 
efficiency  with  decrease  in  volume  ratio  is  attributed  to  the  increase  in  wall  surface  area  for  heat 
transfer  per  unit  of  gas  volume  as  channel  height  decreases.  The  shift  of  stable  equivalence  ratio 
range  towards  progressively  richer  equivalence  ratios  with  decrease  in  channel  height,  shown  in 
Figure  C.3,  is  attributed  to  this  increased  heat  recirculation  efficiency;  the  more  effective 
preheating  provided  by  increased  heat  recirculation  efficiency  permits  the  combustion  of 
progressively  richer  reactant  mixtures. 
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Figure  C.5:  Heat  recirculation  efficiency  is  shown  for  each  of  the  scaled  reactors,  including 
those  scaled  by  channel  height  (VR(h)  0.75  and  VR(h)  1.25)  and  those  scaled  by  channel 
length  (VR(L)  0.75  and  VR(L)  1.25),  as  a  function  of  equivalence  ratio. 

Heat  recirculation  efficiency  decreases  slightly  with  increasing  channel  length,  where  the 
shortest  reactor  length  VR(L)  0.75  produced  efficiencies  of  25%  and  22%  versus  22%  and  20% 
observed  for  the  longest  reactor  length  VR(L)  1.25  at  (j)=2.0  and  (j)=2.2.  This  trend  is  attributed  to 
greater  amounts  of  heat  loss  from  longer  reactors  with  larger  external  surface  area.  Additionally, 
flame  fronts  stabilize  further  downstream  with  increase  in  equivalence  ratio.  This  brings  the 
locations  of  peak  gas  temperatures  in  adjacent  channels  closer  to  the  axial  center,  which 
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minimizes  heat  loss  and  temperature  decrease  of  hot  combustion  gases  prior  to  preheating 
reactants  in  adjacent  channels.  Thus  the  heat  recirculation  efficiency  of  the  longest  reactors, 
VR(L)  1.25  and  VR  1.0,  appear  to  converge  with  increasing  equivalence  ratios,  where  the  flame 
fronts  in  both  reactors  are  approaching  the  channel  axial  length  centerlines.  The  shortest  reactor 
channel  length,  VR(L)  0.75,  has  limited  equivalence  ratio  stability  range,  but  the  peak  reactor 
temperatures  are  in  close  proximity  to  the  preheat  zones  of  adjacent  channels  at  the  points  where 
stability  is  attained. 

Figure  C.6  shows  the  calculated  heat  recirculation  efficiency  for  each  of  the  five  scaled  reactors 
as  a  function  of  inlet  velocity.  A  decrease  in  efficiency  is  observed  for  all  reactors  except  that 
with  the  smallest  channel  height,  VR(h)  0.75,  which  shows  a  peak  efficiency  at  w=150  cm/s.  The 
decrease  in  efficiency  exhibited  by  most  reactors  is  attributed  to  the  decrease  in  residence  time 
with  increasing  inlet  velocity,  which  impacts  fuel  conversion  and  heat  release.  A  similar  trend 
was  observed  in  Barra  et  al  [22].  The  different  trend  exhibited  by  VR(h)  0.75  is  likewise 
attributed  to  fuel  conversion  efficiency  and  peak  temperatures  in  this  reactor. 


As  was  observed  in  Figure  C.5  for  the  variation  of  equivalence  ratio,  heat  recirculation  efficiency 
is  observed  to  increase  with  decreasing  channel  height  and  channel  length,  with  the  largest 
difference  observed  between  VR(h)  0.75  and  VR(h)  1.25. 
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Figure  C.6:  Heat  recirculation  efficiency  is  shown  for  each  of  the  scaled  reactors,  including 
those  scaled  by  height  (VR(h)  0.75  and  VR(h)  1.25)  and  those  scaled  by  length  (VR(L)  0.75 
and  VR(L)  1.25),  as  a  function  of  inlet  velocity. 
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Temperature  Ratio 


Excess  enthalpy  reactors  recirculate  energy  to  increase  the  enthalpy  of  the  reaction  zone,  which 
permits  the  combustion  of  mixtures  at  equivalence  ratios  beyond  the  conventional  limits  of  the 
inlet  mixture.  As  a  result,  flame  temperatures  may  exceed  the  adiabatic  flame  temperature 
predicted  by  equilibrium  based  upon  initial  reactant  temperature.  The  extent  to  which 
temperatures  become  superadiabatic  depends  upon  the  preheat  temperature  achieved  by  the 
reactants,  and  therefore  the  amount  of  heat  recirculated  as  well  as  the  amount  of  heat  loss  from 
the  reactor.  Peak  reactor  temperatures  are  therefore  reflective  of  the  effectiveness  of  heat 
recirculation.  A  temperature  ratio  is  defined  as  a  measure  of  superadiabicity,  or  extent  of 
temperature  increase  above  the  adiabatic  equilibrium  temperature,  as 


_  .  Peak  gas  or  wall  temperature 

Temperature  ratio  — - - -  .  (C.8) 

Adiabatic  equilibrium  temperature 


Temperature  ratio  is  defined  in  terms  of  both  peak  gas  and  peak  wall  temperatures,  where  the 
latter  is  more  readily  measured  during  experimentation. 
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Figure  C.7:  Separated  by  a  dashed  line,  gas  (above)  and  wall  (below)  temperature  ratios 
are  shown  for  each  of  the  scaled  reactors,  including  those  scaled  by  height  (VR(h)  0.75  and 
VR(h)  1.25),  those  scaled  by  length  (VR(L)  0.75  and  VR(L)  1.25),  and  original  reactor 
geometry  (VR  1.0),  as  a  function  of  equivalence  ratio. 


Figure  C.7  shows  gas  and  wall  temperature  ratios  for  each  of  the  five  scaled  reactors  as  a 
function  of  equivalence  ratio.  All  temperature  ratios  are  above  unity,  indicating  superadiabatic 
peak  gas  and  wall  temperatures  at  all  equivalence  ratios.  Gas  temperature  ratios  are  also 
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significantly  higher  than  wall  temperature  ratios  at  all  points,  where  the  latter  are  the 
temperatures  typically  measured  during  experimentation.  Both  gas  and  wall  temperature  ratios 
increase  with  equivalence  ratio,  suggesting  that  higher  levels  of  preheat,  and  therefore 
increasingly  superadiabatic  conditions,  are  required  to  achieve  reactor  stability  at  equivalence 
ratios  increasingly  above  the  conventional  flammability  limit.  Gas  temperature  ratios  are  highest 
in  the  reactors  with  decreased  channel  height  and  length,  VR(h)  0.75  and  VR(L)  0.75,  and  lowest 
in  the  reactors  with  increased  channel  height  and  length,  VR(h)  1.25  and  VR(L)  1.25. 
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Figure  C.8:  Separated  by  a  dashed  line,  gas  (above)  and  wall  (below)  temperature  ratios 
are  shown  for  each  of  the  scaled  reactors,  including  those  scaled  by  height  (VR(h)  0.75  and 
VR(h)  1.25),  those  scaled  by  length  (VR(L)  0.75  and  VR(L)  1.25),  and  original  reactor 
geometry  (VR  1.0),  as  a  function  of  inlet  velocity. 

Figure  C.8  shows  gas  and  wall  temperature  ratios  for  each  of  the  five  scaled  reactors  as  a 
function  of  inlet  velocity.  All  temperature  ratios  are  above  unity,  indicating  superadiabatic  peak 
gas  and  wall  temperatures  at  all  inlet  velocities.  Gas  temperature  ratios  are  also  significantly 
higher  than  wall  temperature  ratios  at  all  points.  Gas  temperature  ratios  are  highest  in  the  reactors 
with  decreased  channel  height  and  length,  VR(h)  0.75  and  VR(L)  0.75,  and  lowest  in  the  reactors 
with  increased  channel  height  and  length,  VR(h)  1.25  and  VR(L)  1.25.  This  correlates  with  heat 
recirculation  efficiency  trends  observed  in  Figure  C.6.  Gas  temperature  ratios  generally  show  a 
slight  increase  with  inlet  velocity  at  the  lowest  tested  inlet  velocities,  and  a  decrease  with  inlet 
velocity  at  the  highest  tested  values.  This  is  attributed  to  a  low  rate  of  heat  release  at  the  lowest 
firing  rates,  where  heat  losses  may  have  a  significant  impact  on  reactor  temperatures,  and 
incomplete  fuel  conversion  at  the  highest  firing  rates  due  to  decreased  residence  time. 
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Hydrogen  Conversion  Efficiency 


The  goal  of  hydrocarbon  reforming  is  the  production  of  syngas  consisting  primarily  of  carbon 
monoxide  and  hydrogen,  of  which  hydrogen  is  often  the  higher  valued  product.  Hydrocarbon 
reforming  is  carried  out  at  rich  conditions,  where  partial  oxidation  will  occur  due  to  oxidizer 
deficiency  of  the  reactants.  Equilibrium  predicts  an  increase  in  hydrogen  product  concentration 
with  increasing  reactant  equivalence  ratio.  The  equilibrium  yield  is  not  necessarily  achievable  in 
a  practical  reactor,  however,  due  to  limited  residence  time  for  reaction.  Hydrogen  production  in 
the  counterflow  reactor  is  defined  in  terms  of  hydrogen  conversion  efficiency  as 


rr  .  rr-  ■  H7  product  concentration 

Hyconversion  efficiency  — -  (C.9) 

H2  equilibrium  concentration 


where  thermodynamic  equilibrium  provides  a  theoretical  maximum  concentration,  and  hydrogen 
conversion  efficiency  quantifies  the  actual  hydrogen  yield  relative  to  this  value. 


no 
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Figure  C.9:  Hydrogen  conversion  efficiencies  are  shown  for  each  of  the  scaled  reactors, 
including  those  scaled  by  height  (VR(h)  0.75  and  VR(h)  1.25)  and  those  scaled  by  length 
(VR(L)  0.75  and  VR(L)  1.25),  as  a  function  of  equivalence  ratio. 


Figure  C.9  shows  hydrogen  conversion  efficiencies  for  each  of  the  five  scaled  reactors  as  a 
function  of  equivalence  ratio.  All  reactors  show  a  decrease  in  conversion  efficiency  with 
increasing  equivalence  ratio.  The  results  for  the  original  reactor  geometry,  VR  1.0,  are  in 
agreement  with  experiment  findings  that  hydrogen  conversion  efficiency,  as  defined  above, 
decreases  from  approximately  95%  to  75%  as  equivalence  ratio  is  increased  from  (j)=1.8  to  (j)=2.4 
[4].  Computational  reactor  conversion  efficiencies  exceed  100%  at  the  lowest  tested  equivalence 
ratios,  where  hydrogen  product  concentration  exceeds  the  yield  predicted  by  equilibrium  for 
reactants  at  300K.  Equilibrium  predicts  a  shift  in  primary  hydrogen-containing  products  from 
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water  towards  diatomic  hydrogen  as  reactant  temperature  is  increased.  Since  reactants  are 
preheated  within  the  counterflow  reactor,  hydrogen  production  in  excess  of  that  predicted  by 
equilibrium  for  300K  products,  and  therefore  conversion  efficiencies  above  100%,  are 
reasonable. 


Trends  observed  in  hydrogen  conversion  efficiencies  are  similar  to  those  observed  in  temperature 
ratio  and  heat  recirculation  efficiency  variation  with  height  and  length  scaling:  hydrogen 
conversion  efficiencies  are  slightly  higher  in  the  reactors  with  decreased  channel  height  and 
length,  VR(h)  0.75  and  VR(L)  0.75,  as  compared  to  original  reactor  geometry,  VR  1.0,  and 
lowest  in  the  reactors  with  increased  channel  height  and  length,  VR(h)  1.25  and  VR(L)  1.25. 
Overall,  however,  there  is  little  difference  in  hydrogen  conversion  efficiency  between  the 
reactors  at  most  tested  equivalence  ratios.  This  suggests  that  heat  recirculation  and  reactant 
preheating  are  necessary  for  flammability  limit  extension,  but  do  not  significantly  impact  product 
composition  in  many  cases.  This  finding  is  in  agreement  with  a  study  of  burner-stabilized  flames 
of  preheated,  premixed  methane  and  air  by  Smith  et  al  [23].  The  lack  of  significant  impact  of 
reactor  length  on  product  composition  is  also  suggested  by  the  previous  findings  of  Kaisare  et  al 
[20],  who  concluded  that  burner  length  is  unlikely  to  impact  reactor  operation  significantly 
because  homogeneous  combustion  is  largely  localized. 


Figure  C.10  shows  hydrogen  conversion  efficiencies  for  each  of  the  five  scaled  reactors  as  a 
function  of  inlet  velocity.  All  reactors  except  for  VR(h)  0.75  show  an  increase  in  conversion 
efficiency  with  inlet  velocity  at  low  firing  rates,  and  a  decrease  with  inlet  velocity  at  high  firing 
rates.  This  suggests  that  optimal  operation  occurs  around  the  middle  of  the  tested  inlet  velocity 
range,  where  firing  rate  is  sufficiently  high  to  dominate  over  external  losses  but  not  high  enough 
to  limit  conversion  by  decreased  residence  time.  The  results  for  VR(h)  0.75  suggest  that  this 
reactor  may  be  able  to  stabilize  at  higher  inlet  velocities  beyond  the  tested  range. 
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Figure  C.10:  Hydrogen  conversion  efficiencies  are  shown  for  each  of  the  scaled  reactors, 
including  those  scaled  by  height  (VR(h)  0.75  and  VR(h)  1.25)  and  those  scaled  by  length 
(VR(L)  0.75  and  VR(L)  1.25),  as  a  function  of  inlet  velocity. 
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C.3  Conclusions 


A  computational  study  was  undertaken  to  determine  the  effect  of  geometric  scale  on  fuel 
reforming  in  the  meso-scale  parallel-plate  counterflow  reactor.  Length  and  height  of  the  reactor 
channels  were  scaled  relative  to  the  original  dimensions  of  the  counterflow  reactor  that  has  been 
used  in  previous  analytical  and  experimental  investigations.  Results  were  quantified  in  terms  of 
steady  state  reactor  operation,  heat  recirculation  efficiency,  temperature  ratio,  and  hydrogen 
conversion  efficiency. 

Reactor  operating  range  was  found  to  be  strongly  impacted  by  channel  scaling.  Decreasing 
channel  height  resulted  in  a  shift  in  stable  operating  range  towards  higher  equivalence  ratios. 
Decreasing  channel  length  resulted  in  a  significant  narrowing  of  the  stable  equivalence  ratio  and 
inlet  velocity  ranges.  Most  stable  operating  points  were  above  the  conventional  rich  flammability 
limit  and  adiabatic  flame  speed  of  methane  and  air  combustion  based  upon  inlet  conditions, 
which  highlights  the  importance  of  heat  recirculation  in  achieving  operation  of  the  counterflow 
reformer  at  these  conditions. 

Heat  recirculation  efficiency  and  temperature  ratios  were  found  to  be  significantly  higher  for 
reactors  with  decreased  channel  height.  These  findings  indicate  more  effective  heat  recirculation 
with  smaller  channel  height,  where  channel  wall  surface  area-to-gas  volume  ratio  is  increased 
and  the  transverse  distance  for  mass  and  heat  diffusion  is  reduced.  Channel  length  showed  a  less 
significant  impact  on  heat  recirculation,  with  shorter  channel  lengths  producing  slightly  higher 
heat  recirculation  efficiencies  and  temperature  ratios  than  longer  lengths. 

Hydrogen  conversion  efficiency  decreased  with  increasing  equivalence  ratio  for  all  reactors,  with 
slightly  higher  conversion  efficiencies  found  for  smaller  channel  heights  and  lengths.  Hydrogen 
conversion  efficiency  initially  increased  with  inlet  velocity  at  the  lowest  tested  firing  rates,  and 
decreased  at  the  highest  tested  inlet  velocities,  suggesting  that  an  optimal  firing  rate  exists. 

Scaled  down  reactors  were  achieved  by  decreasing  reactor  channel  height  and  length,  although 
the  latter  resulted  in  a  significantly  narrowed  stable  range.  The  results  of  this  study  highlight  the 
effectiveness  of  small  scale  reactor  channels  in  achieving  heat  recirculation  and  significantly 
superadiabatic  temperatures.  These,  in  turn,  permit  the  effective  extension  of  flammability  limits 
to  stable  operation  at  rich  equivalence  ratios.  While  increasing  heat  recirculation  may  be  required 
to  significantly  extend  the  stable  operating  range,  the  results  of  this  study  suggest  that 
superadiabatic  temperatures  serve  to  extend  the  flammable  range  of  reactor  operating  points  but 
otherwise  do  not  strongly  impact  hydrogen  conversion  efficiencies. 
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Appendix  D:  Operation  of  the  counterflow  reactor  at  lean  conditions 

Although  the  focus  of  the  ARO  project  was  to  investigate  conversion  of  rich  mixtures,  additional 
funding  was  available  to  explore  other  aspects  of  the  counterflow  reactor.  The  work  reported  in 
this  section  establishes  that  the  reactor  may  also  be  an  effective  lean  combustor  for  gaseous  and 
liquid  fuels.  Analytical  modeling  indicates  the  mode  of  stabilization  of  lean  flames  in  the 
counterflow  reactor.  This  work  has  been  published  as  E.L.  Belmont,  I.  Schoegl,  J.L.  Ellzey,  Proc. 
Combust.  Inst.  34  (2013)  3361-3367,  and  a  second  publication  entitled  “Lean  heptane  and 
propane  combustion  in  a  non-catalytic  parallel-plate  counterflow  reactor”  has  been  accepted  for 
publication  in  Combustion  and  Flame. 

The  use  of  small  scale  combustors  in  conjunction  with  mechanical  or  electrical  conversion 
devices,  such  as  microturbines  or  thermoelectrics,  has  been  proposed  to  address  the  increasing 
need  for  portable  power  [1-3].  Lean  reactant  mixtures  are  favorable  because  they  can  permit 
increased  levels  of  fuel  conversion  and  thermal  efficiencies  as  compared  to  reactant  mixtures 
nearer  to  stoichiometric.  In  order  to  achieve  a  given  firing  rate,  however,  increased  volumetric 
throughput  is  required  as  reactants  become  progressively  leaner.  Therefore,  the  optimal 
equivalence  ratio  for  combustor  operation  is  dependent  upon  the  application.  In  order  for  a 
combustor  to  be  widely  applicable,  a  broad  operating  range  is  desirable. 

Heat  recirculation  can  be  used  to  extend  conventional  lean  flammability  limits  and  increase 
burning  speed  by  increasing  the  enthalpy  in  the  reaction  zone  above  that  of  adiabatic  conditions. 
Various  heat-recirculating  lean  combustor  designs  have  been  proposed  and  studied,  and  all 
operate  on  the  principal  of  transferring  heat  through  a  solid.  In  some  cases,  the  combustor 
consists  of  a  solid  matrix,  such  as  a  reticulated  foam  or  packed  bed  of  spheres  [4-8].  In  others  it 
consists  of  parallel  channels  [9],  such  as  the  counterflow  reactor  used  in  this  study  [10]. 
Theoretical  analyses  and  experimental  studies  of  channel  reactors  have  shown  that 
superadiabatic  burning  velocities  and  broadened  flammability  limits  can  be  achieved  [11-15]. 

Other  studies  of  the  counterflow  reactor  have  demonstrated  its  ability  to  convert  fuel-rich 
reactants  of  methane  [16],  propane  [17]  and  heptane  [18]  to  hydro  gen -rich  syngas.  The  studies 
presented  here  used  the  same  counterflow  reactor  for  the  study  of  fuel-lean  mixtures  of  methane, 
as  well  as  increasingly  complex  fuels:  propane  and  heptane.  The  latter,  being  a  liquid  fuel,  has 
importance  as  a  single  component  surrogate  for  logistical  fuels  [19].  Furthermore,  the 
demonstration  of  fuel-lean  liquid  combustion  in  the  counterflow  reactor  is  an  important 
demonstration  of  its  fuel-flexibility  for  field  applications  because  liquid  fuels  are  generally  safer 
and  easier  to  transport  than  gaseous  fuels. 

Counterflow  combustor  stability  and  emissions  measurements  of  carbon  monoxide  (CO), 
nitrogen  oxides  (NOx)  and  unbumed  hydrocarbons  (UHCs)  were  analyzed  in  this  study  for  lean 
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operation  on  propane  and  heptane.  Stability  was  compared  for  reactor  operation  in  the  lean  and 
rich  regimes  for  both  fuels.  Reactor  channel  wall  temperatures  were  compared  to  adiabatic 
equilibrium  temperatures  to  determine  superadiabicity.  The  results  of  lean  methane,  propane  and 
heptane  combustion  in  the  counterflow  reactor  were  examined  to  understand  the  influence  of 
superadiabatic  operation  on  operating  range  and  emissions. 

D.l  Analytical  modeling  of  lean  combustion  in  the  counterflow  reactor 

An  analytical  model  was  developed  by  Schoegl  and  Ellzey  [10]  to  investigate  the  operating 
principle  of  the  counterflow  reactor  and  obtain  qualitative  predictions  of  its  behavior  in  response 
to  changes  in  operating  conditions.  A  complete  description  of  the  model  can  be  found  in  [10] 
where  it  was  developed  for  analysis  of  the  counterflow  reactor  as  a  fuel  reformer.  A  summary  of 
the  model  and  description  of  its  recalibration  for  counterflow  reactor  operation  on  lean 
methane/air  reactants  is  described  below. 

The  purpose  of  the  one-dimensional  model  is  to  provide  insight  into  the  general  performance  of 
the  counterflow  reactor.  The  model  assumes  constant  thermophysical  properties,  and  combustion 
is  modeled  using  activation  energy  asymptotics.  Radiative  heat  losses  are  included  through  a 
boundary  condition  in  the  wall  temperature  solution.  The  model  solves  for  non-dimensional  wall 
temperature  (Tw)  and  gas  temperature  (77)  along  a  normalized  x  coordinate.  The  location  of  the 
reaction  zone  is  determined  indirectly  to  be  where  the  local  reaction  rate  produces  a  burning 
speed  that  is  equal  to  the  inlet  velocity,  which  is  the  condition  for  stability.  Operating  conditions 
are  specified  by  equivalence  ratio,  inlet  temperature  To,  and  ambient  temperature  7oo.  Properties 
of  the  reactor  are  specified  using  non-dimensional  parameters  which  can  be  found  in  [10]. 

The  gas  temperatures  in  two  adjacent  channels  are  coupled  by  heat  transfer  through  a  common 
conducting  wall.  Equation  D.l  gives  the  normalized  steady-state  energy  conservation  equations 
and  includes  conduction  and  convection,  transport,  external  losses  and  chemical  reactions,  given 
by 


Sf-^[a+T)rw-i(T1  +  r2)-^ 


(D.la) 


d2Ti 

dx2 


=  ±a^i  +  jL (r  r  )  + 


W; 


(D.  lb) 


where  indices  i  =  1,2  denote  individual  channels  and  Tad  is  the  normalized  adiabatic  flame 
temperature.  Equation  D.2  is  the  species  conservation  equation  which  includes  diffusion, 
transport  and  reactions,  and  gives  the  concentration  of  the  limiting  reactant  species  yi  as 
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(D.2) 


1  d2V;  ,  U  dVi  1 

- —  —  + - — - Wi 

Le  dx2  e  dx  e2 


where  Le  is  the  Lewis  number  and  it  is  the  inlet  velocity  non-dimensionalized  by  the  adiabatic 
flame  speed  of  the  reference  case.  Equation  D.3  represents  the  reaction  terms  wi,  which  are 
modeled  as  point  sources  using  activation  energy  asymptotics  [20]  by 

Wi  W  eS(x-  xci)exp  (D'3) 

where  S(x  —  xc  i)  are  Dirac  delta  functions  at  the  combustion  locations  xc  i.  The  normalized 
activation  energy  [1  is  equal  to  Ze/y,  where  the  Zeldovich  number  Ze  is  given  by  E'aAT'ad/R'T'gd, 
the  prime  (')  denotes  a  dimensional  quantity,  E’a  is  the  activation  energy,  AT'ad  is  the  adiabatic 
temperature  rise  of  a  mixture  over  the  unburned  mixture  7^,  and  R'  is  the  universal  gas  constant. 
The  scaling  factor  y  is  equal  to  AT'ld/T^,e^  where  is  the  adiabatic  equilibrium  temperature  of 
the  reference  case. 


Chemical  reactions  are  modeled  using  single-step,  first-order,  irreversible  Arrhenius  kinetics. 
The  activation  energy  is  calibrated  for  lean  methane/air  mixtures  using  laminar  flame  speeds  for 
various  preheating  levels  at  the  equivalence  ratios  of  interest  calculated  with  Cantera  [21]  and 
GRI  3.0  [22].  Adiabatic  flame  temperature  and  flame  speed  at  the  reference  condition  (|>  =  0.5 
and  Tg  =  177.0  °C  are  used  for  the  non-dimensionalization  of  the  model.  Equation  D.4  is  utilized 
to  obtain  the  activation  energy  by  fitting  [3  at  each  equivalence  ratio 


\ui\  =  exp 


P(Tc,i  Tad) 

2Tc>i/T  ad 


(D.4) 


where  Tc,i  are  adiabatic  flame  temperatures  normalized  by  the  reference  temperature.  The 
normalized  adiabatic  flame  temperature  Tad  is  equal  to  1  at  tj)  =  0.50;  at  other  equivalence  ratios, 
it  is  the  adiabatic  flame  temperature  resulting  in  a  mass  flux  equal  to  the  reference  case  in 
detailed  simulations,  which  matches  flame  speeds  within  the  context  of  the  constant  density 
approximation  used  for  this  model.  Table  D.l  gives  the  conditions  resulting  in  constant  mass  flux 
and  the  resultant  values  for  the  calibrated  activation  energies. 


Due  to  the  simplifications  used  in  the  analytical  model,  its  use  is  limited  to  qualitative 
predictions  of  reactor  behavior.  It  is,  however,  able  to  illustrate  the  stabilization  mechanism  of 
the  reactor  and  the  trends  of  temperature  and  reaction  zone  location  in  response  to  variations  in 
equivalence  ratio  and  inlet  velocity. 
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Table  D.l:  Calibrated  activation  energies  derived  from  calculated  temperatures  and  flame 
speeds 


<t> 

n 

(°C) 

T'ad 

(°C) 

A T'ad 
(°C) 

u'  E'a 

(cm/s)  (kJ/mol) 

0.42 

316 

1285 

696 

20.2 

249.7 

0.44 

281 

1295 

741 

19.1 

247.2 

0.46 

247 

1305 

786 

17.9 

244.2 

D.2  Experimental  method 


Experimental  Apparatus 


Figure  D.l  shows  the  counterflow  reactor,  which  is  the  same  design  used  in  previous  reforming 
studies  [16-18].  It  contains  four  4  mm  high  parallel  channels  constructed  of  1  mm  thick  silicon 
carbide  (SiC)  walls.  The  reactor  is  17.3  cm  long,  with  a  91.5  mm  long  main  section  and  two  40.8 
mm  long  reactor  heads.  Inlets,  outlets,  end  plugs  which  prevent  the  mixing  of  combustion 
products  and  reactants,  and  two  6.4  mm  wide  sections  of  SiC  porous  foam  which  act  as  flow 
straighteners  and  flame  arresters  are  located  in  the  reactor  heads.  Two  alumina  walls  spaced  33.6 
mm  apart  enclose  the  channels.  B-type  thermocouples  (TCs)  are  inserted  into  the  channels  and 
sit  flush  with  the  channel  wall  in  order  to  prevent  flame  holding.  Premixed  fuel  and  air  are 
provided  to  the  reactor  using  mass  flow  controllers  for  methane,  propane  and  air,  and  a  pump  for 
heptane.  Heptane  is  vaporized  prior  to  entering  the  reactor.  A  detailed  description  of  the  heptane 
vaporization  system  used  in  this  study  is  available  in  Belmont  et  al.  [18]. 
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Figure  D.l:  Four-channel  counterflow  reactor  with  dimensions  and  internal  components 

indicated 


Concentrations  of  carbon  monoxide  (CO),  nitrogen  oxides  (NOx),  and  unbumed  hydrocarbons 
(UHCs)  were  measured  in  the  reactor  exhaust  using  Rosemount  Analytical  analyzers.  Samples 
were  drawn  through  an  uncooled  quartz  probe  with  a  2  mm  inner  diameter  and  inert  tubing,  and 
dried  before  analysis.  A  previous  study  [17]  showed  slightly  lower  product  species 
concentrations  from  outside  channels,  and  minimal  variation  of  composition  within  a  single 
channel  except  near  the  alumina  walls.  In  consideration  of  these  findings,  exhaust  samples  in  this 
study  were  taken  from  the  center  channels  at  the  middle  of  the  channel  exit.  Therefore,  the  data 
produced  in  this  study  is  expected  to  be  representative  of  interior  channels  in  scaled  up  reactors. 
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Experimental  Method 


The  start-up  procedure  began  by  burning  a  near- stoichiometric  mixture  of  the  fuel  of  interest  at 
the  reactor  outlets  at  an  inlet  velocity  of  50  cm/s.  The  reaction  fronts  propagated  upstream  in  the 
channels  once  the  reactor  heads  had  warmed  sufficiently,  and  stabilized  downstream  of  the 
porous  SiC  flow  straighteners.  Operating  conditions  were  then  gradually  adjusted  to  the 
reference  conditions  for  each  fuel,  and  were  maintained  at  those  values  until  temperatures 
stabilized.  References  conditions  for  methane  were  (f>  =  0.44  and  u  =  125  cm/s,  and  reference 
conditions  for  propane  and  heptane  were  <f>  =  0.41  and  u  =  125  cm/s. 

The  warm-up  phase  took  approximately  45  minutes.  Experiments  included  the  determination  of 
stable  operating  conditions  and  the  exhaust  gas  concentrations  of  CO,  NOx  and  UHCs  for  each 
fuel  over  a  range  of  inlet  velocities  and  equivalence  ratios.  Test  conditions  were  considered 
stable  if  combustion  was  sustained  in  the  main  reactor  section  for  ten  minutes  at  constant 
temperatures.  Reactor  wall  temperatures  were  monitored  in  each  of  the  four  channels  throughout 
the  experiments.  Gas  temperatures  were  not  measured  directly  due  to  the  tendency  of 
thermocouples  inserted  into  the  channel  to  act  as  flame  holders.  Peak  wall  temperatures  were 
recorded  at  stable  reactor  conditions  and  showed  less  than  ±5°C  variation  during  that  time  and 
while  exhaust  gas  samples  were  taken.  Reactor  temperatures  were  limited  to  less  than  1300°C  in 
order  to  avoid  damage  to  the  reactor. 

Following  the  warm-up  phase,  equivalence  ratio  or  inlet  velocity  was  adjusted  to  the  desired  test 
values.  Inlet  velocity  was  specified  at  standard  conditions  of  25°C  and  1  atm.  Emissions  results 
were  obtained  over  a  range  of  inlet  velocities  with  equivalence  ratio  held  constant  at  (/>  =  0.44  for 
methane  and  (j)  =  0.41  for  propane  and  heptane,  and  over  a  range  of  equivalence  ratios  with  inlet 
velocity  held  constant  at  125  cm/s.  Emissions  concentrations  in  parts  per  million  (ppm)  are 
reported  as  measured  on  a  dry  basis. 

Each  operating  point  was  tested  twice.  Total  uncertainty  was  calculated  at  each  point  as  the  root- 
sum- square  of  the  contributing  uncertainties.  The  uncertainty  of  emissions  results  is  due  to 
repeatability  of  measurements  as  calculated  using  a  student-t  distribution,  the  calibration  gases 
and  the  analyzers.  Average  uncertainties  of  CO,  NOx  and  UHC  measurements  for  methane, 
propane  and  heptane  are  ±10,  2  and  16,  ±13,  2  and  16  ppm  and  ±28,  2  and  1  ppm,  respectively. 
Inlet  velocity  uncertainty  is  due  to  uncertainties  of  the  mass  flow  controllers  for  air,  methane  and 
propane,  the  pump  for  heptane,  and  the  channel  cross-sectional  area.  Equivalence  ratio 
uncertainty  is  due  to  contributions  from  the  flow  controllers  and  pump.  Average  uncertainties  of 
inlet  velocity  and  equivalence  ratio  are  ±10  cm/s  and  ±0.03  for  methane  and  propane,  and  are 
±10  cm/s  and  ±0.04  for  heptane.  Temperature  uncertainty  is  attributed  to  contributions  from 
repeatability  of  measurements,  the  thermocouples,  and  the  data  acquisition  system.  Average 
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temperature  uncertainties  for  methane,  propane  and  heptane  are  estimated  to  be  ±12°C,  ±16°C 
and  ±20°C,  respectively. 


D.3  Results  and  Discussion 
D.3.1  Analytical  model  results 

Equivalence  ratio  and  inlet  velocity  can  be  changed  independently  while  stability  is  maintained 
within  the  reactor.  There  are  no  predetermined  flame  stabilization  points  within  the  channel; 
therefore  the  combustion  zone  freely  adjusts  to  a  new  position  within  the  channel  following  a 
change  in  equivalence  ratio  or  inlet  velocity.  The  combustion  zone  moves  to  a  position  where 
sufficient  heat  is  gained  from  the  adjacent  channel  and  reactants  are  preheated  to  the  point  where 
the  burning  velocity  equals  the  inlet  velocity,  and  stability  is  maintained. 
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Figure  D.2:  Analytical  model  (AM)  gas  (Tg)  and  wall  (Tw)  temperatures  and  experimental 
(EXP)  temperatures  measured  along  axial  position  in  channel  for  varying  (|)  and  u  =  125 

cm/s 


Figure  D.2  shows  the  results  of  the  analytical  model  for  gas  and  wall  temperature  profiles  for 
varying  equivalence  ratio  from  0.42  to  0.46  with  inlet  velocity  held  constant  at  125  cm/s,  as  well 
as  experimental  temperature  measurements  taken  along  the  axial  length  of  a  reactor  channel 
while  operating  at  these  conditions,  normalized  by  axial  position  and  overlaid  onto  the  analytical 
domain.  Figure  D.3  shows  the  results  of  the  model  and  experimental  measurements  for  varying 
inlet  velocity  from  75  to  175  cm/s  with  equivalence  ratio  held  constant  at  0.44. 
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Figure  D.3:  Analytical  model  (AM)  gas  (Tg)  and  wall  (Tw)  temperatures  and  experimental 
(EXP)  temperatures  measured  along  axial  position  in  channel  for  varying  u  and  <|)  =  0.44 

In  both  cases,  the  inlet  temperature  was  specified  at  577  °C  to  account  for  preheating  that  occurs 
in  the  reactor  heads.  The  combustion  zone  is  identified  by  a  sharp  increase  in  gas  temperature, 
while  the  wall  temperature  shows  a  broader  profile.  The  combustion  zone  adjusts  its  location 
towards  the  center  of  the  reactor  following  a  decrease  of  (|>  or  an  increase  in  u  in  Figures  D.2  and 
D.3,  respectively.  For  increasing  <j)  (at  fixed  u),  the  temperature  of  the  flame  increases  and  less 
preheat  is  required  for  the  mixture  to  react  and  so  the  flame  stabilizes  nearer  the  entrance  than  for 
the  cases  at  lower  (j)  which  require  greater  preheating  (Fig.  D.2).  For  increasing  u  (at  fixed  (j)),  the 
reactants  travel  further  into  the  channel  before  reaching  the  critical  temperature  for  combustion 
(Fig.  D.3).  Through  this  adjustment  in  combustion  zone  location  a  burning  velocity  is  attained 
that  equals  the  inlet  velocity. 

Wall  temperatures  also  respond  to  changes  in  (j)  or  u,  both  in  magnitude  and  broadness  of 
temperature  profile.  As  (j)  is  increased  (Fig.  D.2)  or  u  is  decreased  (Fig.  D.3)  and  the  combustion 
zone  shifts  upstream,  hot  products  flow  over  the  remaining  length  of  the  channel.  This  condition, 
as  well  as  conduction  through  the  wall,  results  in  the  broad  profile  predicted  by  the  model. 
Similarly,  as  (j)  is  decreased  or  u  is  increased  and  the  reaction  zone  shifts  towards  the  center  of 
the  reactor,  the  central  portion  of  the  wall  receives  more  heat  and  the  peak  temperature  zone  in 
the  wall  is  narrower.  Experimental  temperature  measurements  are  in  qualitative  agreement  with 
the  results  of  the  analytical  model  in  terms  of  relative  values  at  the  tested  operating  points. 

While  the  analytical  results  can  only  be  interpreted  qualitatively  due  to  simplifications  in  the 
model,  they  do  accurately  capture  temperature  trends  and  the  general  conditions  under  which 
flashback  and  blowoff  occur,  which  are  at  high  (j)-low  u  and  low  (j)-high  u,  respectively.  The 
results  of  this  study  can  also  be  compared  to  a  previous  study  in  which  the  counterflow  reactor 
was  operated  as  a  fuel  reformer,  and  the  analytical  model  was  used  to  predict  its  behavior  under 
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fuel  rich  conditions  [10].  In  that  study  the  model  similarly  predicted  changes  in  combustion  zone 
location  in  response  to  changes  in  (j)  and  u.  In  the  case  of  fuel  rich  conditions  the  combustion 
zone  migrated  downstream  for  increases  in  u,  as  in  this  study.  However,  the  zone  migrated 
downstream  for  increases  in  (j)  because  more  preheating  is  required  to  attain  the  same  peak 
reaction  temperature  as  a  rich  mixture  becomes  progressively  richer,  analogous  to  the  case  for  a 
lean  mixture  becoming  progressively  leaner.  Therefore  the  model  indicates  that  the  stabilization 
method  for  rich  and  lean  mixtures  in  the  counterflow  reactor  is  similar. 

D.3.2  Experimental  results 

Operating  Range 

Methane 

Figure  D.4  shows  the  stable  operating  range  of  the  reactor  operating  on  methane/air  reactants  in 
terms  of  inlet  velocity  and  both  standard  and  normalized  equivalence  ratio  [23].  The  combustion 
process  becomes  unstable  due  to  either  flashback,  in  which  the  reaction  front  propagates 
upstream  past  a  section  of  porous  flow  straightener,  blowoff,  in  which  the  reaction  front 
propagates  out  of  the  outlet  of  the  reactor  channel,  or  extinction.  In  addition  to  the  data  presented 
in  this  paper  for  operation  in  the  lean  regime,  data  from  an  earlier  study  [16]  on  rich  methane/air 
mixtures  is  also  shown.  The  stable  operating  points  on  both  the  rich  and  lean  sides  are  beyond 
the  conventional  flammability  limits. 

In  terms  of  standard  equivalence  ratio,  the  stable  operating  range  on  the  lean  side  is  considerably 
narrower  than  that  on  the  rich  side.  Stability  is  strongly  influenced  by  the  temperature  of  the 
products  and  the  resultant  preheating  of  the  reactants  in  the  neighboring  channel.  It  is  well 
known  that  in  comparison  to  the  rich  side  of  stoichiometric,  the  temperature  of  the  products 
changes  more  significantly  with  standard  equivalence  ratio  on  the  lean  side  resulting  in  a  smaller 
range  of  stable  conditions.  Due  to  the  asymmetry  inherent  in  the  standard  equivalence  ratio 
definition,  the  stability  range  is  also  presented  in  terms  of  normalized  equivalence  ratio.  In 
normalized  form,  the  width  of  the  stable  range  in  the  rich  regime  is  approximately  60%  broader 
than  the  lean,  which  is  attributed  to  differences  in  chemical  kinetics  between  the  two  regimes. 
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Figure  D.4:  Stability  map  for  counterflow  reactor  operating  on  lean  (this  study)  and  rich 
([16])  methane/air  reactants,  with  both  stable  (ST)  and  unstable  points  due  to  flashback 
(FB),  blowoff  (BO)  and  extinction  (EX)  shown  as  a  function  of  inlet  velocity  and  both 
standard  and  normalized  equivalence  ratios 

Propane 

Figure  D.5  is  a  stability  plot  for  counterflow  reactor  operation  on  propane/air  reactants,  where 
operating  conditions  are  defined  in  terms  of  inlet  velocity  and  both  standard  and  normalized 
equivalence  ratio.  Normalized  equivalence  ratio,  defined  as  cj)/(  1  +cj>),  permits  direct  comparison 
of  rich  and  lean  operating  regimes  because  it  compensates  for  the  inherent  asymmetry  in  the 
definition  of  standard  equivalence  ratio  [23]. 
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Figure  D.5:  Stability  map  for  counterflow  reactor  operation  on  lean  and  rich  [17] 
propane/air  reactants,  with  stable  (ST)  and  unstable  points  due  to  flashback  (FB),  blow-off 
(BO)  and  extinction  (EX)  shown  as  a  function  of  inlet  velocity  and  both  standard  and 
normalized  equivalence  ratios.  Lean  and  rich  flammability  limits  for  propane/air  reactants 

at  inlet  conditions  are  indicated. 
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The  conditions  where  the  reactor  operates  stably,  as  well  as  conditions  where  operation  becomes 
unstable,  are  shown.  Combustor  instability  can  occur  as  flashback,  in  which  the  reaction  front 
propagates  upstream  past  a  section  of  porous  flow  straightener;  blow-off,  in  which  the  reaction 
front  propagates  out  of  the  outlet  of  the  reactor  channel;  or  extinction,  in  which  the  rate  of  heat 
loss  exceeds  the  rate  of  heat  release  and  combustion  is  not  sustained.  In  addition  to  stability  data 
acquired  for  lean  combustor  operation,  previously  obtained  data  for  reactor  operation  on  fuel- 
rich  mixtures  of  propane  [17]  are  included  in  the  plot.  Taken  together,  a  broader  understanding 
of  stable  reactor  operation  is  attained. 

Some  operating  points  were  not  tested  at  high  inlet  velocities  in  order  to  avoid  damage  to  the 
reactor  due  to  high  temperatures.  Lean  operation  is  achieved  at  minimum  and  maximum  propane 
standard  equivalence  ratios  of  0.36  and  0.44,  over  a  range  of  inlet  velocities  from  60  to  225  cm/s. 
These  conditions  correspond  to  Reynolds  numbers  of  approximately  100-600  based  upon  inlet 
conditions,  putting  reactor  operation  in  the  laminar  regime.  Stable  operation  in  the  lean  regime  is 
entirely  below  the  conventional  lean  flammability  limit  of  propane,  <f>  =  0.51,  while  the  stable 
rich  range  spans  from  below  to  above  the  conventional  rich  flammability  limit,  <j)  =  2.5  [24],  In 
normalized  form,  the  width  of  the  stable  rich  regime  is  approximately  35%  broader  than  the  lean. 

Heptane 

Figure  D.6  is  a  stability  plot  for  counterflow  reactor  operation  on  heptane/air  reactants,  where 
operating  conditions  are  defined  in  terms  of  inlet  velocity  and  both  standard  and  normalized 
equivalence  ratio.  The  operating  conditions  where  the  reactor  operates  stably,  as  well  as 
conditions  where  operation  becomes  unstable  due  to  flashback,  blow-off  or  extinction,  are 
shown.  In  addition  to  stability  data  acquired  for  lean  combustor  operation,  previously  obtained 
data  for  reactor  operation  on  fuel -rich  mixtures  of  heptane  [18]  are  included  in  the  plot. 
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Figure  D.6:  Stability  map  for  counterflow  reactor  operation  on  lean  and  rich  [18] 
heptane/air  reactants,  with  stable  (ST)  and  unstable  points  due  to  flashback  (FB),  blow-off 
(BO)  and  extinction  (EX)  shown  as  a  function  of  inlet  velocity  and  both  standard  and 
normalized  equivalence  ratios.  Lean  and  rich  flammability  limits  for  heptane/air  reactants 

at  inlet  conditions  are  indicated. 

Lean  operation  is  achieved  at  minimum  and  maximum  heptane  standard  equivalence  ratios  of 
0.37  and  0.46,  over  a  range  of  inlet  velocities  from  100  to  225  cm/s.  The  minimum  stable  inlet 
velocity,  100  cm/s,  is  significantly  higher  than  previously  tested  fuels,  methane  and  propane, 
which  stabilized  at  60  cm/s.  Destabilization  by  extinction  at  the  lower  velocity  limit,  which 
occurs  at  75  cm/s  for  heptane  combustion  at  (j)  =0.41  in  the  counterflow  reactor,  has  been 
observed  in  previous  studies  as  well  [16-18].  This  mode  of  destabilization  results  from  the  low 
rate  of  enthalpy  input  at  low  inlet  velocities,  while  the  rate  of  heat  loss  is  a  function  of  reactor 
temperature,  which  does  not  change  proportionately.  Stable  operation  in  the  lean  regime  is 
significantly  below  the  conventional  lean  flammability  limit  of  heptane,  <j>  =  0.56,  while  the 
stable  rich  range  spans  from  below  to  slightly  above  the  conventional  rich  flammability  limit,  <f>  = 
3.8  [24].  In  both  standard  and  normalized  form,  the  width  of  the  stable  rich  regime  is 
significantly  broader  than  the  lean. 

Emissions 


Exhaust  emissions  of  CO,  NOx  and  UHCs,  and  peak  wall  temperature  measurements  are 
reported  for  lean  reactor  operation  on  methane,  propane  and  heptane.  Results  are  presented  for 
two  sets  of  tests:  variation  of  inlet  velocity  with  equivalence  ratio  held  constant,  and  variation  of 
equivalence  ratio  with  inlet  velocity  held  constant. 
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Effect  of  Inlet  Velocity 


For  a  specified  equivalence  ratio,  the  equilibrium  temperature  is  fixed  and  is  not  affected  by 
volumetric  flow  rate.  In  a  practical  reactor,  however,  increasing  the  inlet  velocity  at  a  specified 
equivalence  ratio  increases  the  firing  rate  while  altering  the  ratio  of  volumetric  heat  release  to 
external  heat  loss,  which  is  primarily  affected  by  convection  from  the  outside  surface  area.  In 
addition,  the  residence  time  in  the  reactor  decreases  with  increasing  velocity.  These  two  factors 
of  temperature  and  residence  time  are  important  in  the  formation  of  emissions. 

Figure  D.7a  shows  CO  exhaust  concentrations  and  peak  measured  wall  temperatures  for  lean 
reactor  operation  on  heptane,  propane  and  methane  as  a  function  of  inlet  velocity.  Figure  D.7b 
shows  UHC  exhaust  concentrations  and  peak  wall  temperatures  for  propane  and  methane.  For 
heptane,  UHCs  were  undetected  at  all  points  and  are  not  included  in  Figure  D.7b.  In  order  to 
facilitate  direct  comparison  of  the  data  presented  in  this  study  to  other  bodies  of  work,  these 
results  are  presented  in  terms  of  standard  equivalence  ratio.  Inlet  velocity  varies  from  60  to  225 
cm/s  with  equivalence  ratio  held  constant  at  (f>  =  0.44  for  methane  and  <f>  =  0.41  for  propane,  and 
from  100  to  200  cm/s  for  heptane  at  (j>  =  0.41.  Calculated  adiabatic  equilibrium  temperatures  for 
each  fuel  at  the  tested  equivalence  ratios  are  shown.  Propane  and  heptane  data  were  acquired  at  (f> 
=0.41,  as  highlighted  in  the  stability  maps  shown  in  Figs.  D.5  and  D.6,  where  broad  ranges  of 
stable  inlet  velocities  were  achieved  for  both  fuels.  Methane  flames  stabilized  over  a  broad  range 
of  inlet  velocities  at  an  equivalence  ratio  of  (f>  =0.44,  as  highlighted  in  Fig.  D.4.  Although  this 
small  difference  in  equivalence  ratio  complicates  the  comparison  of  the  different  fuels,  it  is  still 
instructive  to  examine  the  trends  associated  with  changes  in  inlet  velocity  which  significantly 
affect  reactor  temperatures.  This  effect  is  attributed  to  an  increase  in  the  ratio  of  volumetric  heat 
release  within  the  reactor  to  heat  losses  from  the  surface  of  the  reactor  as  inlet  velocity  increases. 
The  change  in  reactor  temperatures  strongly  impacts  reaction  rates,  and  therefore  emissions. 

Additionally,  regions  of  combustor  operation  are  defined  based  upon  peak  wall  temperatures, 
and  are  indicated  in  Figures  D.7a  and  D.7b.  While  peak  gas  temperatures  are  expected  to  be 
superadiabatic  at  all  operating  points  based  upon  previous  theoretical  analyses  [10],  peak  wall 
temperatures  are  a  readily  measured  indicator  of  reactor  conditions.  Peak  wall  temperatures  are 
the  lowest  temperatures  experienced  by  the  reaction  zone;  this  has  significant  implications  for 
reaction  chemistry  and  quenching,  and  therefore  emissions.  The  operating  regions  are  designated 
as  subadiabatic,  near-adiabatic,  and  superadiabatic  depending  on  the  value  of  the  peak  reactor 
wall  temperature  relative  to  the  adiabatic  temperature,  where  near-adiabatic  operation  is 
identified  where  peak  reactor  temperatures  are  within  approximately  50°C. 

Propane  combustion  produces  peak  reactor  wall  temperatures  that  vary  from  890  to  1233°C  at 
the  lowest  and  highest  tested  inlet  velocities,  respectively,  and  are  superadiabatic  at  inlet 
velocities  above  100  cm/s.  The  positive  correlation  of  peak  wall  temperature  with  inlet  velocity 
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is  also  shown  for  methane  as  well  as  in  previous  studies  [16-18].  This  is  a  result  of  increased 
rate  of  heat  release  within  the  reactor  with  increasing  velocity  compared  to  rate  of  heat  loss  to  the 
surroundings.  Heptane  combustion  produces  a  similar  temperature  trend,  with  peak  wall 
temperatures  increasing  from  1046  to  1149°C  with  increasing  inlet  velocity.  Heptane  has  a 
narrower  stable  operating  range  than  propane  and  methane  at  the  tested  equivalence  ratios; 
instability  occurs  at  inlet  velocities  below  100  cm/s  and  above  200  cm/s. 

The  CO  concentration  from  heptane  combustion  decreases  from  259  ppm  to  133  ppm  with 
increasing  inlet  velocity.  CO  concentration  from  propane  also  decreases  with  increasing  inlet 
velocity  above  75  cm/s,  with  minor  variation  at  the  highest  tested  inlet  velocities  that  is  within 
calculated  uncertainty.  At  the  lowest  tested  inlet  velocity  of  50  cm/s,  the  concentration  of  CO  is 
significantly  less  than  the  peak  value  of  123  ppm  that  occurs  at  75  cm/s.  UHCs  from  propane  are 
detected  at  inlet  velocities  of  50  and  75  cm/s. 

The  results  of  propane  combustion  show  temperature  and  CO  trends  that  are  similar  to  methane. 
Detection  of  UHCs  at  the  lowest  tested  inlet  velocities  indicates  incomplete  fuel  breakdown  for 
both  fuels.  The  absence  of  measured  UHCs  for  heptane  is  attributed  to  the  narrower  stable 
operating  range  of  heptane  within  the  counterflow  reactor,  which  is  limited  to  a  minimum  inlet 
velocity  of  100  cm/s.  It  is  anticipated  that  UHC  would  be  detected  if  heptane  flames  could  be 
stabilized  in  the  reactor  at  lower  inlet  velocities,  following  trends  observed  for  methane  and 
propane,  however  the  current  reactor  geometry  does  not  support  heptane  flames  at  inlet  velocities 
below  100  cm/s  at  <f>  =  0.41. 
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Figure  D.7:  (a)  CO  exhaust  concentration  and  peak  wall  temperature  (Tw)  for  lean 
counterflow  combustor  operation  on  heptane,  propane  and  methane,  and  (b)  UHC  exhaust 
concentration  and  peak  wall  temperature  measured  for  lean  counterflow  combustor 
operation  on  propane  and  methane  are  shown.  UHC  are  undetected  at  all  tested  heptane 
operating  points  and  are  not  shown.  A  range  of  inlet  velocities  are  tested  with  equivalence 
ratio  held  constant  at  ^  =  0.41  for  heptane  and  propane,  and  </>  =  0.44  for  methane. 
Calculated  adiabatic  equilibrium  temperatures  (Tad)  for  the  tested  equivalence  ratios  of 
each  fuel  are  shown.  Subadiabatic,  near-adiabatic  and  superadiabatic  operating  regions, 
defined  in  terms  of  peak  wall  temperature,  are  also  shown. 
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Figure  D.8  shows  NOx  exhaust  concentrations  for  combustion  of  lean  heptane,  propane  and 
methane.  Results  are  presented  for  inlet  velocities  from  60  to  225  cm/s  with  equivalence  ratio 
held  constant  at  (j>  =  0.41  for  propane  and  <j>  =  0.44  for  methane,  and  from  100  to  200  cm/s  for 
heptane  at  <f>  =  0.41.  Peak  measured  wall  temperatures  and  calculated  adiabatic  equilibrium 
temperatures  are  shown.  Additionally,  regions  of  combustor  operation,  defined  by  subadiabatic, 
near-adiabatic  and  superadiabatic  peak  wall  temperatures,  are  indicated. 
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Figure  D.8:  NOx  exhaust  concentration  and  peak  wall  temperature  (Tw)  measured  for  lean 
counterflow  combustor  operation  on  heptane,  propane  and  methane.  Inlet  velocity  is  varied 
with  equivalence  ratio  held  constant  at  <j>  =  0.41  for  heptane  and  propane,  and  (/>  =  0.44  for 
methane.  Adiabatic  equilibrium  temperatures  (Tad)  for  the  range  of  tested  equivalence 
ratios  of  each  fuel  are  shown.  Subadiabatic,  near-adiabatic  and  superadiabatic  operating 
regions,  defined  in  terms  of  peak  wall  temperature,  are  also  shown. 


NOx  concentration  is  below  13  ppm  at  all  tested  points.  The  highest  NOx  measurements  occur  at 
125  cm/s  for  propane,  and  the  concentrations  decrease  with  increasing  and  decreasing  inlet 
velocities  from  that  point.  The  same  trend  is  observed  for  methane  combustion.  The  initial 
increase  in  NOx  with  inlet  velocity  is  attributed  to  temperature  dependence,  whereas  the 
variation  at  higher  inlet  velocities  is  largely  within  the  estimated  uncertainty.  A  slight  decrease  at 
the  highest  tested  inlet  velocities  is  attributed  to  decreased  residence  times  at  these  conditions. 
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Heptane  combustion  produces  the  highest  NOx  concentration  at  the  lowest  tested  inlet  velocity. 
Concentration  decreases  with  an  increase  in  inlet  velocity,  and  remains  constant  within  the 
estimated  uncertainty  over  the  remainder  of  the  tested  inlet  velocities. 

Effect  of  Equivalence  Ratio 

For  a  fixed  inlet  velocity,  a  change  in  equivalence  ratio  changes  the  reactor  firing  rate  while 
maintaining  an  approximately  constant  residence  time.  The  necessary  firing  rate  is  determined 
by  process  requirements,  while  environmental  considerations  dictate  acceptable  emissions  levels. 

Figure  D.9  shows  CO  exhaust  concentrations  and  peak  wall  temperatures  for  lean  reactor 
operation  on  heptane,  propane  and  methane.  Equivalence  ratios  are  varied  from  0.38  to  0.45  for 
heptane,  0.37  to  0.44  for  propane  and  0.41  to  0.48  for  methane,  with  inlet  velocity  held  constant 
at  u  =  125  cm/s.  Calculated  adiabatic  equilibrium  temperatures  for  each  fuel  at  the  tested 
equivalence  ratios  are  also  shown.  The  widths  of  stable  equivalence  ratio  ranges  are  the  same  for 
each  fuel  at  an  inlet  velocity  of  125  cm/s,  where  all  emissions  were  measured  for  equivalence 
ratio  variation.  However,  as  can  be  seen  in  the  stability  maps  shown  in  Figs.  D.4,  D.5  and  D.6, 
the  upper  and  lower  limits  of  these  stable  ranges  vary  between  the  different  fuels.  Methane  in 
particular  exhibits  a  stable  range  that  is  shifted  more  significantly  towards  stoichiometric  as 
compared  to  the  other  two  tested  fuels.  Because  there  are  no  flameholders  in  the  counterflow 
reactor,  stable  operation  is  obtained  when  the  inlet  flow  is  balanced  by  the  burning  rate,  which  is 
a  function  of  equivalence  ratio,  preheat  temperature,  heat  losses  and  recirculation,  and  fuel 
chemistry. 
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Figure  D.9:  CO  exhaust  concentration  and  peak  wall  temperature  (Tw)  measured  for  lean 
counterflow  combustor  operation  on  heptane,  propane  and  methane  with  varying  </>  and 
inlet  velocity  held  constant  at  u  =  125  cm/s  are  shown,  as  well  as  calculated  adiabatic 
equilibrium  temperatures  (Tad)  for  the  range  of  tested  equivalence  ratios  of  each  fuel 
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Methane,  propane  and  heptane  combustion  result  in  modest  increases  of  peak  wall  temperatures 
with  equivalence  ratio.  For  example,  for  heptane  and  propane,  temperatures  increase  from  1046 
to  1149°C  and  1071  to  1146°C,  respectively.  The  positive  correlation  of  temperature  with 
equivalence  ratio  for  both  fuels  is  also  in  agreement  with  trends  predicted  by  chemical 
equilibrium.  Comparison  of  measured  and  equilibrium  temperatures  shows  that  peak  wall 
temperatures  are  superadiabatic  at  all  tested  methane  and  propane  equivalence  ratios.  Peak  wall 
temperatures  for  heptane  are  superadiabatic  at  equivalence  ratios  below  0.42,  and  are 
subadiabatic  at  the  highest  tested  equivalence  ratios.  Peak  wall  temperatures  as  a  function  of 
equivalence  ratio  are  generally  within  50°C  of  the  adiabatic  equilibrium  temperature  for  all  fuels, 
and  are  considered  to  be  in  the  near-adiabatic  regime,  thus  the  three  operating  regimes  identified 
in  Figures  D.7  and  D.8  are  not  identified  in  Figures  D.9  and  D.10. 

The  measured  CO  concentration  for  propane  combustion  initially  increases  slightly  with 
increasing  equivalence  ratio  from  <f>  =  0.37  to  (j>  =  0.39,  however  the  variation  in  this  range  is 
within  calculated  uncertainty.  At  higher  equivalence  ratios,  CO  concentration  decreases  with 
increasing  equivalence  ratio  from  a  peak  value  of  106  ppm  to  50  ppm  at  the  highest  tested 
equivalence  ratio  of  0.44.  Likewise,  the  measured  CO  concentration  for  heptane  combustion 
decreases  from  a  peak  value  of  373  ppm  at  the  lowest  tested  equivalence  ratio  of  0.38  to  99  ppm 
at  the  highest  tested  equivalence  ratio  of  0.45.  This  trend  is  in  agreement  with  the  results  for 
methane,  which  show  a  decrease  in  CO  concentration  with  increasing  equivalence  ratio  across 
the  tested  range.  As  with  methane,  this  negative  correlation  is  in  contrast  to  chemical 
equilibrium,  which  indicates  an  increase  in  CO  concentration  as  equivalence  ratio  approaches 
stoichiometry.  However,  the  measured  values  for  heptane  and  propane  combustion  are 
significantly  in  excess  of  equilibrium,  which  predicts  less  than  1  ppm  of  CO  at  all  tested 
equivalence  ratios.  Mesoscale  reactors  are  particularly  prone  to  reaction  quenching  at  the  channel 
walls  as  compared  to  larger  reactors  due  to  increased  surface  area-to-volume  ratios,  which  leads 
to  incomplete  combustion.  An  additional  challenge  of  high  area-to-volume  ratios  in  mesoscale 
reactors  is  the  increased  heat  loss  at  channel  walls,  which  produces  lower  reaction  temperatures 
and  therefore  decreased  reaction  rates.  Reactor  temperatures  and  heat  loss  are  also  highly 
dependent  on  inlet  conditions,  including  inlet  velocity  and  equivalence  ratio,  which  significantly 
impact  the  rate  of  volumetric  heat  release.  The  observed  CO  trend  is  attributed  to  a  dependence 
on  reactor  temperature,  which  increases  with  equivalence  ratio  and  drives  the  oxidation  of  CO 
toward  C02  and  equilibrium. 

Figure  D.10  shows  NOx  exhaust  concentrations  for  lean  reactor  operation  on  heptane,  propane 
and  methane.  Results  are  presented  for  varying  equivalence  ratio  with  inlet  velocity  held  constant 
at  u  =  125  cm/s.  Again,  peak  wall  temperatures  and  calculated  adiabatic  equilibrium 
temperatures  are  shown.  The  measured  NOx  concentration  for  heptane  and  propane  combustion 
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increases  with  increasing  equivalence  ratio  from  5  ppm  at  (f>  =  0.37  to  17  ppm  at  </;  =  0.44.  The 
same  trend  is  observed  for  methane  combustion. 


Unbumed  hydrocarbons  are  undetected  at  all  tested  equivalence  ratios  for  lean  heptane,  propane 
and  methane  combustion,  indicating  complete  fuel  breakdown. 


0.36  0.38  0.4  0.42  0.44  0.46  0.48  0.5 

Equivalence  ratio 

Figure  D.10:  NOx  exhaust  concentration  and  peak  wall  temperature  (Tw)  measured  for 
lean  counterflow  combustor  operation  on  heptane,  propane  and  methane  with  varying  (f 
and  inlet  velocity  held  constant  at  u  -  125  cm/s  are  shown,  as  well  as  adiabatic  equilibrium 
temperatures  (Tad)  calculated  for  the  range  of  tested  equivalence  ratios  of  each  fuel 


Effect  of  Superadiabatic  Operation 


The  peak  reactor  wall  temperatures  measured  during  experiments  are  categorized  as 
subadiabatic,  near-adiabatic  and  superadiabatic  relative  to  the  calculated  adiabatic  equilibrium 
temperature  of  the  unburned  fuel  mixture,  and  are  labeled  accordingly  in  Figs.  D.7  and  D.8  In  a 
reactor  that  operates  without  heat  recirculation,  temperatures  will  necessarily  be  subadiabatic 
everywhere  because  of  unavoidable  heat  losses  from  the  system.  In  a  heat-recirculating  reactor, 
however,  temperatures  in  excess  of  the  adiabatic  equilibrium  temperature  can  be  achieved 
through  conduction  and  radiation  by  the  solid  structure  of  the  reactor,  and  through  convection  by 
the  gas.  Previous  theoretical  analyses  of  the  counterflow  reactor  indicate  that  peak  gas 
temperatures  within  the  reactor  channels  are  significantly  in  excess  of  wall  temperatures  [10]. 
Furthermore,  experiments  and  computations  of  stationary  flames  in  ducts  constructed  of  low 
thermal  conductivity  materials  show  that  wall  temperatures  in  such  systems  are  significantly  less 
than  centerline  peak  temperatures  and  adiabatic  temperatures  [25-27].  While  peak  gas 
temperatures  at  the  centerlines  of  the  reactor  channels  are  therefore  expected  to  be  significantly 
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higher  than  peak  wall  temperatures,  the  latter  are  more  readily  measured  without  interfering  with 
reactor  dynamics,  such  as  by  flameholding.  Furthermore,  reactor  wall  temperatures  are  quite 
important  in  determining  quenching  and  heat  loss  rates.  Therefore,  peak  wall  temperatures 
provide  an  important  and  practical  metric  by  which  reactor  performance  can  be  evaluated  and 
monitored.  The  results  of  this  study  reveal  common  trends  in  fuel  emissions  as  the  peak  wall 
temperatures  change  with  varying  inlet  conditions  of  equivalence  ratio  or  inlet  velocity.  Regions 
of  operation  can  be  distinguished  in  terms  of  peak  wall  temperature,  where  distinct  trends  are 
noted  depending  on  the  value  of  the  peak  wall  temperature  relative  to  the  calculated  adiabatic 
equilibrium  temperature.  The  results  highlight  the  importance  of  achieving  significantly 
superadiabatic  peak  temperatures  within  the  reactor  channels,  which  will  result  in  near-adiabatic 
or  superadiabatic  wall  temperatures  due  to  heat  losses  at  external  surfaces,  in  achieving  favorable 
emissions. 

The  current  analysis  of  peak  reactor  wall  temperature  relative  to  calculated  adiabatic  equilibrium 
temperature  is  distinguished  from  a  global  definition  of  superadiabatic  performance  that  is  used 
in  the  discussion  of  heat  recirculating  reactors,  which  refers  to  the  peak  temperature  in  the  gas 
phase  of  the  reactor  relative  to  the  calculated  adiabatic  equilibrium  temperature.  The  peak  wall 
temperatures  observed  in  this  study  approach  or  are  above  the  calculated  adiabatic  equilibrium 
temperatures  for  each  fuel.  Therefore,  based  upon  previous  theoretical  and  experimental  analyses 
of  heat  recirculating  reactors  which  account  for  heat  losses  at  the  wall  at  all  conditions,  all 
reported  cases  in  this  study  most  likely  exhibit  superadiabatic  temperatures  at  the  centerline  of 
the  reactor  channel,  and  operation  would  therefore  be  deemed  superadiabatic  by  global 
definition.  The  current  analysis,  however,  focuses  on  the  categorization  of  peak  wall  temperature 
relative  to  calculated  equilibrium  temperature,  and  uses  the  designation  of  subadiabatic,  near- 
adiabatic  and  superadiabatic  to  convey  the  relative  magnitudes  of  these  values. 

Inlet  velocity  has  a  significant  effect  on  the  peak  wall  temperature,  which  is  as  much  as  200°C 
greater  than  the  adiabatic  temperature  at  the  conditions  tested  in  this  study.  In  Figures  D.7  and 
D.8,  regions  of  subadiatic,  near-adiabatic,  and  superadiabatic  peak  wall  temperature  conditions 
have  been  identified.  Stable  operation  for  heptane  was  not  obtained  at  significantly  subadiabatic 
wall  temperature  conditions.  For  methane  and  propane,  high  levels  of  UHCs  are  detected  at 
subadiabatic  wall  temperature  conditions,  and  NOx  concentrations  are  low.  These  results  are 
consistent  with  incomplete  combustion. 

At  conditions  where  near-adiabatic  wall  temperatures  were  measured,  CO  peaks  and  then 
decreases  with  increasing  inlet  velocity  for  all  fuels  and  UHCs  are  not  detected  at  the  highest 
tested  inlet  velocities  in  this  regime.  The  variation  of  NOx  among  the  different  fuels  in  this 
region  is  not  consistent:  concentrations  increase  for  propane  and  methane  and  decrease  for 
heptane.  It  should  be  noted,  however,  that  all  NOx  concentrations  are  less  than  15  ppm. 
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At  conditions  where  superadiabatic  peak  wall  temperatures  were  measured,  the  emissions  of  all 
reported  species  are  fairly  constant  with  the  exception  that  CO  for  heptane/air  continues  to 
decrease. 

A  negative  correlation  of  CO  with  temperature  is  observed  in  all  data,  except  those  points  where 
UHCs  are  detected  for  methane  and  propane  (Figures  D.7  and  D.9).  These  products  are  far  from 
equilibrium  values,  which  predict  less  than  1  ppm  for  both  CO  and  UHCs.  In  contrast,  no  UHCs 
are  detected  for  heptane  combustion.  Likewise,  no  significantly  subadiabatic  wall  temperatures 
are  measured  at  any  stable  heptane  operating  points.  The  stable  range  of  inlet  velocities  is 
narrower  for  heptane  as  compared  to  methane  and  propane.  Flashback  occurs  at  a  higher  inlet 
velocity  for  heptane  compared  to  the  other  two  fuels,  and  the  lowest  stable  velocity  is  limited  to 
approximately  where  the  measured  peak  wall  temperatures  become  subadiabatic.  In  contrast, 
methane  and  propane  are  able  to  stabilize  at  significantly  subadiabatic  wall  temperatures,  but 
UHCs  are  measured  at  these  points  (Figure  D.7). 

Heptane  is  likewise  limited  at  the  high  end  of  the  tested  inlet  velocity  range  by  blow-off  at  225 
cm/s,  whereas  methane  and  propane  are  stable  at  this  condition.  Peak  wall  temperature  values  are 
lower  for  heptane  than  methane  and  propane  across  the  tested  range  of  inlet  velocities  (Figures 
D.7  and  D.8).  Flame  temperature  and  speed  are  strongly  dependent  on  reactant  temperature; 
therefore  peak  wall  temperature  and  burning  rate  within  the  reactor  are  expected  to  be  highly 
dependent  on  the  extent  of  counterflow  heat  exchange  from  combustion  products  in  order  to 
preheat  counterflowing  reactants  in  adjacent  channels.  Since  the  counterflow  reactor  requires 
burning  rate  to  equal  reactant  velocity  in  order  to  operate  stably,  insufficient  preheating  of 
reactants  leads  to  heptane  blow-off  at  high  inlet  velocities. 

The  results  of  lean  methane,  propane  and  heptane  combustion  suggest  significant  changes  in 
stability  and  emissions  behavior  when  the  peak  wall  temperatures  become  subadiabatic. 
Emissions  concentrations  indicate  that  low  wall  temperatures  in  the  subadiabatic  peak  wall 
temperature  regime  do  not  promote  sufficient  reaction  rates  for  complete  fuel  breakdown,  and 
increasing  wall  temperatures  into  the  adiabatic  and  superadiabatic  regimes  drive  reactions  toward 
equilibrium.  Emissions  and  stability  range  results  suggest  that  combustor  operation  with  peak 
wall  temperatures  in  the  near-adiabatic  and  superadiabatic  regimes  is  desirable  for  allowing 
increased  firing  rates  and  decreased  emissions.  The  most  significant  decline  in  emissions  is  seen 
in  the  transition  from  subadiabatic  to  near-adiabatic  peak  wall  temperatures,  while  further 
temperature  increase  into  the  superadiabatic  regime  produces  marginal  additional  advantage.  The 
more  significant  advantage  of  operation  in  the  superadiabatic  wall  temperature  regime  is  found 
in  extension  of  the  stable  combustor  operating  range,  as  highlighted  by  the  absence  of  heptane 
stability  in  the  subadiabatic  wall  temperature  regime.  This  finding  is  in  agreement  with  a 
previous  investigation  of  preheated  flames  [28]. 
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In  contrast  to  the  previous  cases  in  which  peak  wall  temperature  was  significantly  affected  by 
changes  in  velocity,  the  peak  wall  temperature  as  a  function  of  equivalence  ratio  (Figures  D.9 
and  D.10)  is  generally  within  50°C  of  the  adiabatic  equilibrium  temperature  at  the  tested  inlet 
velocity  of  125  cm/s  for  each  fuel,  and  is  considered  to  be  in  the  near-adiabatic  regime. 
Emissions  of  CO  decrease  and  NOx  increase  with  increasing  equivalence  ratio  for  all  three  fuels. 
Measured  concentrations  are  significantly  above  equilibrium  values,  and  the  observed  trends  are 
attributed  to  temperature  dependence. 

While  the  data  for  equivalence  ratio  variation  with  an  inlet  velocity  of  125  cm/s  show  peak  wall 
temperatures  in  the  near-adiabatic  regime,  Figures  D.7  and  D.8  indicate  that  an  increase  or 
decrease  in  the  inlet  velocity  may  result  in  a  shift  of  the  peak  wall  temperatures  into  the 
superadiabatic  or  subadiabatic  regimes.  The  effects  of  such  a  shift  on  stability  and  emissions  can 
be  anticipated  based  upon  the  results  for  inlet  velocity  variation. 

D.4  Conclusion 

Lean  premixed  combustion  of  methane,  propane  and  heptane  was  studied  in  a  heat  recirculating 
mesoscale  counterflow  reactor.  This  reactor  design  has  previously  been  used  to  study  the 
conversion  of  rich  methane,  propane  and  heptane  reactant  mixtures  to  hydrogen-rich  synthesis 
gas.  The  current  study  expands  the  established  capability  of  the  counterflow  reactor  to  operate  in 
the  lean  combustion  regime  on  increasingly  complex  fuels.  In  particular,  liquid  fuels  are 
generally  easier  and  safer  to  transport  than  gaseous  fuels,  and  are  studied  because  of  their 
logistical  importance. 

Results  included  reactor  stability  maps  for  methane,  propane  and  heptane  combustion  that  show 
stable  operating  points,  as  well  as  points  where  operation  became  unstable  due  to  flashback, 
blow-off  or  extinction,  for  rich  and  lean  reactant  mixtures.  Emissions  of  CO,  NOx  and  UHC 
were  presented  for  lean  methane,  propane  and  heptane  operation.  Additionally,  operating  regions 
were  defined  in  terms  of  peak  reactor  wall  temperatures.  Similarities  in  trends  between  the  three 
tested  fuels  highlight  the  correlation  of  emissions  with  peak  reactor  wall  temperature,  and 
support  previous  findings  that  a  minimum  firing  rate  is  necessary  for  optimal  performance. 
Further  analysis  highlights  the  importance  of  achieving  peak  wall  temperatures  near,  or  in  excess 
of,  the  calculated  adiabatic  equilibrium  temperature  of  the  fuel  at  inlet  conditions  in  heat¬ 
recirculating  reactors  for  emissions  minimization  and  stable  range  extension. 
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Appendix  E:  Syngas  production  from  burner-stablized  methane/air  flames: 
the  effect  of  preheat 

The  funding  from  ARO  was  leveraged  with  support  from  supplemental  sources  which  permitted 
the  investigation  of  other  questions  related  to  the  conversion  of  a  hydrocarbon  fuel/air  mixture  to 
syngas.  The  work  below  is  published  in  C.H.  Smith,  D.I.  Pineda,  and  J.L.  Ellzey,  “Syngas 
Production  From  Burner- Stabilized  Methane/Air  Flames:  The  Effect  of  Preheated  Reactants”, 
Comb.  Flame  (160)  pp.  557-564,  2013. 

E.l  Introduction 

Heat  recirculating  reactors,  in  which  heat  is  transferred  from  hot  products  to  cold  reactants,  allow 
the  combustion  of  mixtures  with  equivalence  ratios  (4>)  outside  the  conventional  flammability 
limits  and  with  burning  velocities  far  in  excess  of  laminar  flame  speeds  [1-3].  In  addition,  local 
peak  temperatures  can  be  superadiabatic,  i.e  greater  than  those  predicted  by  equilibrium.  Because 
of  these  unique  characteristics,  heat  recirculating  reactors  have  been  proposed  as  low  emissions 
radiant  heaters  and  combustors  [4-6]  and  fuel  reformers.  A  significant  amount  of  work  has  been 
done  on  fuel  reforming  in  particular,  where  a  rich  mixture  of  fuel  and  oxidizer  is  reacted  to 
produce  syngas — a  mixture  of  hydrogen  (H2),  carbon  monoxide  (CO)  and  other  species  [3-11]. 

Generally,  in  these  reactors  both  conduction  and  radiation  are  enhanced  over  gaseous  values  by 
the  presence  of  a  solid  surface,  and  in  some  designs,  such  as  a  packed  bed  or  filtration  reactor, 
mass  transfer  is  enhanced  over  the  gas  phase  levels  [12,  13].  The  complexity  of  the  transport 
processes  in  these  reactors  makes  analysis  and  the  determination  of  the  critical  design  parameters 
difficult. 

Previous  research  has  shown  that  the  stable  operating  range,  defined  as  the  sets  of  (|)  and  inlet 
velocity  (V)  where  stable  combustion  is  attained,  varies  widely  depending  on  the  reactor  design 
even  when  operating  on  the  same  fuel.  For  example,  when  operating  on  methane  stable  operation 
in  a  packed  bed  reactor  was  demonstrated  up  to  <])  =  5.0  [14]  while  for  a  counterflow  reactor  the 
limit  was  at  4>  =  2.5  [15].  When  the  reaction  front  was  stabilized  at  the  interface  of  two  sections 
of  porous  media  the  limit  was  (j)  =  1.9  [10]. 

The  conversion  of  the  reactants  to  syngas  also  depends  on  the  reactor  design.  Experiments 
reforming  methane  to  syngas  with  filtration  reactors  [14,  16],  a  counterflow  reactor  [15]  and 
another  type  of  porous  media  reactor  [7]  showed  a  variety  of  syngas  yields  with  maximum  H2 
yields  ranging  from  40%  to  more  than  70%  amongst  the  different  reactors.  Furthermore, 
experiments  reforming  heptane  to  syngas  with  two  different  reactor  types  [17,  18]  showed 
different  reforming  efficiencies.  The  filtration  reactors  typically  produce  the  highest  syngas 
yields,  but  since  the  stable  operating  conditions  of  the  various  reactors  typically  do  not  overlap,  a 
direct  comparison  is  difficult. 
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Several  computational  studies  have  shown  the  importance  of  conduction,  radiation,  and  solid-to¬ 
gas  heat  transfer  in  these  reactors  [19-22].  The  significant  outcome  of  these  enhanced  processes 
is  that  the  reactants  are  preheated  and  will  therefore  react  even  at  extreme  equivalence  ratios.  In 
these  devices,  the  preheat  temperature  is  a  function  of  both  the  properties  of  the  reactor  and  the 
combustion  process.  In  one  study  [6],  the  inlet  mixture  to  a  porous  reactor  was  preheated 
externally  to  a  controlled  temperature  but  since  additional  heat  transfer  takes  place  within  the 
reactor,  the  effective  preheat  temperature  is  difficult  to  discern.  The  development  of  an  optimal 
design  is  challenging  because  the  level  of  preheat  required  for  conversion  and  the  specific  effect 
of  preheat  temperature  on  the  conversion  are  unclear.  An  additional  complication  is  that 
preheating  sometimes  occurs  under  conditions  where  mass  transfer  is  substantially  enhanced. 

The  purpose  of  our  work  was  to  understand  the  effect  of  preheat  temperature  on  the  conversion 
of  rich  mixtures  of  methane/air  to  syngas.  In  order  to  examine  the  specific  effect  of  preheat 
temperature,  we  investigated  flames  stabilized  on  flat  flame  burners.  In  contrast  to  heat 
recirculating  reactors,  flat  flame  burners  allow  the  control  of  preheat  temperature  and  do  not 
have  enhanced  transport  processes  that  complicate  the  analysis. 

There  has  been  very  substantial  work  on  premixed  methane/air  laminar  flames  stabilized  by 
various  methods  including  on  flat  flame  burners.  Most  studies  have  focused  on  lean  or  near- 
stoichiometric  conditions,  though  there  are  many  published  investigations  of  rich  burner- 
stabilized  flames  as  well.  Most  of  the  experiments  with  rich  flames  were  studies  of  detailed 
flame  structure  conducted  at  low-pressure  [23].  It  is  known  that  kinetics  can  differ  between  low 
and  high  pressure  [23,  24]  so  studies  at  high  pressure  are  important  to  supplement  the  detailed 
low-pressure  flame  investigations.  The  studies  at  atmospheric  pressure  [25-28],  have  typically 
focused  on  flame  structure  or  measurements  of  flame  speed  rather  than  the  concentrations  of 
species  over  a  wide  range  of  conditions  as  is  reported  in  this  work.  In  one  exception  [26],  CO 
concentrations  over  a  range  of  equivalence  ratios  were  measured.  Previous  studies  of  rich 
premixed  methane/air  flames  with  preheated  reactants  have  not  included  exhaust  product 
measurements  with  the  exception  of  NOx  [29-33]. 

In  this  study,  we  examined  experimental  measurements  and  computational  predictions  of  exhaust 
products  for  rich  methane/air  flames  with  and  without  preheated  reactants.  Chemical  kinetic  and 
flame  models  have  not  been  validated  extensively  for  very  rich  equivalence  ratios  [34],  therefore 
we  include  comparisons  of  various  global  flame  characteristics  such  as  standoff  distance,  flame 
temperature,  and  stable  operating  range  to  establish  confidence  in  these  models  to  make  adequate 
predictions.  The  models  are  then  used  to  compare  to  the  experimental  data  from  rich  flames  with 
and  without  preheated  reactants  and  also  to  extend  the  range  of  temperature  conditions  that  are 
investigated. 
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E.2  Numerical  and  analytical  approach 


One  dimensional  simulations  of  burner-stabilized  flames  and  laminar  freely-propagating  flames 
were  performed  using  Cantera  [35]  with  GRI3.0  [36].  GRI3.0  was  initially  modified  to  include  a 
CH*  reaction  using  a  set  of  reactions  by  Hai  Wang  (personal  communication),  though  we  found 
that  the  location  of  the  CH*  and  CH  were  nearly  identical.  Because  of  this  similarity,  the  original 
GRI3.0  mechanism  was  used  for  all  simulations. 

For  a  burner-stabilized  flame,  the  energy  equation  was  solved  with  an  upstream  boundary 
condition  of  constant  temperature  and  a  downstream  boundary  condition  of  zero  temperature 
gradient.  The  species  conservation  equation  was  solved  with  multicomponent  transport.  The 
upstream  species  boundary  condition  is  zero  concentration,  and  the  downstream  condition  was 
zero  gradient.  The  domain  was  1  cm  for  all  simulations. 

For  a  freely-propagating  flame,  the  energy  equation  was  solved  with  upstream  and  downstream 
boundary  conditions  of  zero  temperature  gradient.  The  species  conservation  equation  was  solved 
with  multicomponent  transport.  The  upstream  and  downstream  species  boundary  conditions  were 
zero  gradient.  The  domain  was  30  cm.  The  details  of  the  burner-stabilized  flame  model  and  the 
freely-propagating  flame  model  can  be  found  in  the  Cantera  documentation  [35]. 


In  addition  to  the  detailed  computational  models,  the  analytical  solution  of  Law  [37]  was  used  to 
compute  the  standoff  distance  as  a  function  of  operating  conditions  and  to  understand  the  basic 
behavior  of  burner-stabilized  flames.  The  nondimensional  governing  equations  for  flame 
temperature,  standoff  distance  and  heat  loss  are: 


f2  =  P1) 

J  \TadJ 


~  /  1  1  \ 

J  exp 

.  Ta  \f}  W. 

(E.l) 


f[Tad  +  (Ts  -  fu)\  -Ls=  ftf 


(E.2) 


7  _  fjff-Ts) 
s  e~rxf- 1 


(E.3) 


where  the  nondimensonal  variables  for  mass  flux,  temperature,  standoff  distance  and  heat  loss 
are  given  by: 


f  T  xf°  _  dT 

1  r  Tad  —  Tu  ’  PD’  s  dx 
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and  /  is  the  mass  flux,  f°  is  the  laminar  burning  flux,  T  is  the  temperature,  Tad  is  the  adiabatic 
flame  temperature  based  on  the  unbumed  reactant  temperature,  Tu,  x  is  the  standoff  distance,  p  is 
the  density,  and  D  is  the  mass  diffusion  coefficient. 


To  obtain  dimensional  values  for  the  quantities  of  interest,  laminar  burning  fluxes  and  adiabatic 
flame  temperatures  obtained  from  Cantera  solutions  were  used.  Two  other  parameters,  the 
activation  temperature,  Ta,  and  the  product,  pD ,  required  dimensional  values.  Both  of  these 
parameters  were  set  by  calibration  to  experimental  data.  The  activation  temperature  was  set  to 
15098  K  such  that  the  flame  temperature  given  by  the  analytical  model  matched  the  flame 
temperature,  1850  K,  for  methane/air  at  (j)  =  1.2,  V  =  11  cm/s  as  reported  by  [27].  Setting  pD  = 
5.5E-5  yielded  the  best  fit  to  our  experimental  data  for  standoff  distance  at  (j)  =  1.25.  This  value 
of  pD  corresponded  to  the  mass  average  value  of  pD  for  a  mixture  of  methane  and  air  with  ((>  = 
1.25  at  approximately  1025  K. 


We  also  performed  equilibrium  calculations  for  varying  (j)  and  initial  reactant  temperature  to 
compare  with  data  from  combustion  processes  that  approach  equilibrium  in  the  limit  of  infinitely 
fast  kinetics  and  adiabatic  conditions  [11,  17].  In  these  calculations,  the  initial  temperature  and  (|> 
were  selected  and  enthalpy  and  pressure  were  held  constant  for  equilibration.  The  equilibrium 
calculations  provide  insight  into  the  thermodynamic  properties  of  the  mixtures,  and  equilibrium 
has  been  assumed  to  represent  the  species  concentrations  in  the  post-flame  region  of  burner- 
stabilized  flames  [27,  38].  These  calculations  were  performed  in  Cantera  with  thermodynamic 
data  from  GRI3.0.  Solid  carbon  formation  was  included  in  all  of  the  presented  equilibrium 
calculations.  These  calculations  show  that  no  solid  carbon  is  formed  for  methane/air  mixtures 
until  c[>  exceeds  approximately  3.3. 

The  models  are  used  to  predict  flame  standoff  distance  in  terms  of  the  location  of  maximum  CH 
concentration,  flame  temperature,  stability  limits  and  species  yields.  Experimentally,  we  define  a 
stability  limit  as  (]>,  V  pair  where  if  V  is  increased  the  flame  ceases  to  be  flat.  For  the  burner- 
stabilized  flame  model,  the  limit  is  defined  as  the  maximum  V  for  a  given  <])  that  produces  a 
converged  solution.  An  additional  limit  comparison  is  based  on  the  laminar  flame  speed  of  a 
mixture  where  the  predicted  laminar  flame  speeds  are  compared  to  the  limiting  velocities 
determined  experimentally  at  each  (j).  Lastly,  the  models  and  experiments  were  compared  in 
terms  of  exhaust  species  yields.  A  summary  of  the  conditions  for  comparison  is  shown  in  Table 
E. 


95 


Measured  Valued 

Calculated  Value 

Burner  Stabilized  Flame 

Calculated  Value 

Free  Flame 

Given  <|),  V 

Standoff  distance 

Standoff  distance 

NA 

Given  <|),  V 

Flame  temperature 
from  [27] 

Flame  temperature 

NA 

Given  <|) 

Maximum  V  for  flat 

flame 

Maximum  V  at  which  solution 
converged 

Laminar  flame  speed  for 
given  c[> 

Given  V 

Maximum  cf>  for  flat 
flame 

Maximum  c[)  at  which  solution 
converged 

cf>  at  which  laminar  flame 
speed  equals  given  V 

Given  <|),  V 

Species  yields 

Species  yields 

Species  yield  if  laminar 
flame  speed  =  V  at  given  d 

Table  E.l:  Conditions  for  comparing  models  and  experiment.  In  each  row  the  quantity  of 
interest  for  comparison  is  listed  for  operating  condition(s)  given  in  the  far  left  column. 


E.3  Experimental  approach 
Experimental  Apparatus 

Experiments  were  performed  with  two  different  burners.  Experiments  without  preheat  were 
conducted  on  a  water-cooled  McKenna  burner  (Holthius  &  Associates),  which  has  a  sintered 
porous  bronze  plate  as  the  burner-surface.  The  water  flow  was  ~13  g/s  for  all  experiments, 
efficiently  cooling  the  burner  surface.  The  McKenna  burner  was  used  to  establish  a  baseline  with 
which  to  compare  the  experiments  with  preheated  reactants  and  to  determine  the  accuracy  of  the 
models  since  the  McKenna  burner  is  well-established  for  this  comparison.  We  constructed  a 
second  burner,  similar  to  that  used  in  other  studies  of  flat  flames  [33,  39],  for  experiments  with 
preheated  reactants.  This  burner  was  a  square  channel  mullite  ceramic  matrix  with  1.34  mm 
channels  and  0.35  mm  walls.  The  matrix  was  51  mm  long  and  55  mm  in  diameter.  A  diagram  of 
this  burner  is  shown  in  Fig.  E.l. 
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Figure  E.l:  Diagram  of  ceramic  burner 


As  described  by  Law  [37],  when  the  mass  flux  and  mixture  properties  (cj),  unbumed  gas 
temperature)  are  specified  in  Equation  E.l,  the  flame  temperature  is  determined  and  is 
independent  from  the  surface  temperature  and  standoff  distance.  Since  the  reaction  zone  is 
confined  to  a  thin  region  at  or  very  near  the  flame  temperature,  this  flame  temperature  can  be 
considered  the  controlling  temperature  for  the  chemical  reactions  and  product  species.  Therefore 
the  product  species  distributions  for  the  ceramic  burner,  which  does  not  have  a  controlled  surface 
temperature,  the  McKenna  burner  and  the  burner- stabilized  flame  model  may  be  directly 
compared  when  the  mass  flux  and  mixture  properties  are  identical.  In  addition  they  may  be 
compared  in  terms  of  limit  conditions  because  under  limit  conditions  there  is  no  thermal 
interaction  between  the  flame  and  the  burner  surface. 

The  standoff  distances  are  not  compared  to  the  modeling  results  for  preheated  reactants  because 
the  standoff  distance  is  a  function  of  the  surface  temperature,  which  is  not  controlled  in  these 
experiments  using  the  ceramic  burner.  In  contrast  to  the  McKenna  burner  where  the  surface 
temperature  is  effectively  controlled  by  the  temperature  of  the  cooling  fluid,  in  the  ceramic 
burner  the  surface  temperature  is  a  function  of  operating  conditions.  Stabilization  is  achieved 
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through  heat  loss  to  the  environment  primarily  by  radiation  from  the  burner  surface.  The  required 
heat  loss  and  burner  surface  temperature,  then,  determine  the  standoff  distance. 

In  all  experiments,  a  quartz  tube  of  55  mm  ID  was  used  to  shield  the  flames  from  external  air. 
The  quartz  tubes  were  used  instead  of  inert  gas  shielding  because  the  exhaust  gas  samples  were 
extracted  from  the  center  of  the  flame,  where  the  quartz  tube  boundary  had  no  effect.  We  found 
no  observable  difference  in  any  measurements  by  running  experiments  with  and  without  the 
quartz  tube  as  a  shield. 

The  flows  were  controlled  with  calibrated  Hastings  mass  flow  controllers.  The  air  was 
compressed,  dry  laboratory  air,  and  the  methane  was  CP  grade  from  Air  Liquide.  Preheating  was 
accomplished  by  running  the  premixed  methane/air  through  a  helical  coil  of  stainless  steel  heated 
by  a  Bunsen-bumer  before  entering  the  reactor  housing.  The  temperature  of  the  reactants  was 
monitored  by  a  fine  wire  K-type  thermocouple  placed  within  the  ceramic  plug  upstream  of  the 
preheat  zone  of  a  burner-stabilized  flame.  It  was  confirmed  that  fuel  breakdown  did  not  occur 
due  to  preheating  before  reaching  the  flame  by  sampling  the  gases  after  heating  but  without  a 
flame. 

We  also  built  a  system  to  measure  the  spatial  location  and  intensity  of  CH*  chemiluminescence 
for  the  measurement  of  flame  standoff  distance  and  flatness.  This  system  was  a  digital  webcam 
(Logitech™  C910  HD),  a  vertical  traverse  with  a  micrometer  head,  and  a  Lab  VIEW™  control 
system,  along  with  a  430+10  nm  optical  bandpass  filter. 

Species  measurement 

Measurements  of  species  concentrations  in  the  exhaust  were  taken  with  a  Varian  (Agilent)  CP- 
4900  gas  chromatograph,  which  measures  H2,  CO,  nitrogen  (N2),  carbon  dioxide  (C02), 
methane  (CH4),  ethylene  (C2H4),  acetylene  (C2H2)  and  select  higher  hydrocarbons.  No 
compounds  greater  than  C2  were  detected  in  any  experiments.  The  samples  were  all  drawn  from 
1  cm  above  the  burner  surface  with  a  water-cooled  sample  probe.  The  probe  had  a  1  mm  OD 
injection  point  to  minimize  system  disturbance.  The  sources  of  uncertainty  were  in  the  flow  rates 
of  methane  and  air,  the  calibration  standards,  and  the  gas  chromatograph  calibration. 
Uncertainties  in  species  measurements  were  calculated  using  a  Student-t  distribution,  and 
uncertainties  for  calculations  were  based  on  sequential  perturbation.  The  species  results  are 
reported  in  terms  of  dry  species  concentrations. 

Determination  of  flame  flatness 

When  a  flame  approaches  a  blow-off  condition,  the  flame  loses  its  flatness,  and  when  a  flame  is 
in  a  flashback  condition,  it  is  not  visible.  Since  we  wanted  to  compare  our  experimental  results  to 
one-dimensional  computational  models  of  burner-stabilized  flames,  it  was  important  to  establish 
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a  criterion  for  a  “flat”  flame  that  was  also  not  in  a  flashback  condition.  In  addition,  we  define  the 
stability  limit  as  the  as  the  (|>,  V  pair  beyond  which  a  flame  would  not  be  flat  if  the  inlet  velocity 
was  increased  at  the  given  <|>.  We  found  these  limiting  conditions  by  perturbing  (j)  and  V  while 
monitoring  the  flame  flatness  though  CH*  chemiluminescence  imaging.  A  representative 
grayscale  CH*  chemiliuminescence  image  is  shown  in  Fig.  E.2. 


Figure  E.2:  Pictorial  diagram  describing  the  flatness  metric  used  in  this  study  on  a  highly 

non-flat  flame 

The  flame  flatness  was  quantified  by  finding  the  variance  of  the  full  width  (in  y  direction)  at  half 
maximum  (FWHM)  grayscale  intensity  along  the  x  axis.  By  this  definition,  a  perfectly  flat  flame 
has  zero  FWHM  intensity  variance. 

To  determine  the  experimental  limiting  conditions  (the  conditions  where  increases  in  (j)  or  V 
would  produce  a  nonflat  flame),  velocity  was  held  constant  while  (j)  was  adjusted  by  increments 
of  0.005  for  five  different  values  of  (j)  surrounding  a  suspected  limit  and  images  were  taken  to 
determine  the  FWHM  variance  for  each  condition.  A  limit  value  of  (j)  was  defined  as  a  tested  (j) 
limit  where  the  FWHM  variance  found  in  the  set  of  five  conditions  at  least  doubled.  Flames  with 
operating  conditions  within  the  bounds  of  the  limiting  conditions  were  considered  flat,  and  in 
this  paper  we  only  report  measurements  from  flat  flames. 

E.4  Results  and  discussion 

Previous  research  on  the  conversion  of  rich  fuel/air  mixtures  has  shown  that  optimal  production 
of  syngas  may  occur  beyond  the  conventional  flammability  limits.  Optimization,  however,  of 
reactor  performance  remains  challenging  because  the  level  of  preheat  necessary  is  unclear.  In 
the  following  sections,  we  discuss  experimental,  computational,  and  analytical  results  which 
provide  some  insight  into  the  importance  of  preheating  for  conversion  of  rich  methane/air 
mixtures. 

Comparison  of  experimental  and  numerical  modeling  data 

Since  combustion  models  are  not  extensively  compared  to  experimental  data  in  the  rich  regime 
[34],  we  present  comparisons  of  standoff  distance,  stability,  and  temperature  in  order  to  gain 
confidence  in  the  predictions  of  the  models,  which  are  subsequently  used  to  investigate  a  broader 
range  of  conditions  than  was  investigated  experimentally.  In  this  section  the  computational 
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models  were  compared  to  the  results  from  experiments  on  the  McKenna  burner  and  the 
analytical  model  under  non-preheated  conditions. 

Standoff  distance 


For  a  flat  flame  to  stabilize  on  the  surface  of  a  burner,  the  burning  velocity  of  the  flame  must 
match  the  inlet  velocity.  Matching  of  these  velocities  occurs  because  heat  loss  to  the  burner 
decreases  the  flame  temperature  and  therefore  decreases  the  burning  velocity  as  described  in 
detail  in  the  literature  [37,  40,  41].  The  standoff  distance,  which  we  use  for  comparison  between 
the  model  and  experiment,  is  a  fundamental  characteristic  of  burner-stabilized  flames  and  has 
been  used  by  other  researchers  for  such  comparison  [42,  43]. 


Experimental,  computational,  and  analytical  results  show  that  the  flame  standoff  distance 
decreases  with  increasing  inlet  velocity,  as  shown  in  Fig.  E.3.  In  the  range  of  inlet  velocities 
tested,  the  model  and  experiment  agreed  very  well  in  both  trend  and  magnitude,  and  the 
analytical  model  fits  the  data  very  well.  This  trend  has  been  observed  by  other  researchers  and  is 
explained  in  detail  in  the  literature  [32,  41]. 
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Figure  E.3:  Standoff  distance  vs.  inlet  velocity  at  equivalence  ratio  =  1.25  ±  0.05 


We  also  compared  the  modeling  predictions  with  the  experimental  measurements  of  standoff 
distance  as  a  function  of  <])  while  V  was  held  constant  at  15  cm/s  (Fig.  E.4).  Again,  the  agreement 
between  the  experiments,  the  burner-stabilized  flame  model  and  the  analytical  model  is  very 
good.  The  same  trend  was  observed  by  Furguson  and  Keck  [40],  and  their  standoff  distance 
measurements  were  within  about  1  mm  of  those  reported  here  even  though  the  reactant 
temperature  in  their  experiments  was  50  K  greater  than  the  temperature  of  the  reactants  in  our 
experiments. 
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The  standoff  distance  is  one  of  the  fundamental  global  parameters  that  characterize  burner- 
stabilized  flames.  Over  a  range  of  V  and  (j),  the  burner-stabilized  flame  model  was  a  very  good 
predictor  of  standoff  distance.  Law’s  analytical  model  also  fits  the  experimental  and  modeling 
data  very  well  in  trend  and  in  magnitude  except  for  a  slight  overprediction  near  stoichiometric. 


Figure  E.4:  Standoff  distance  vs.  equivalence  ratio  at  inlet  velocity  =  15  +  0.2  cm/s 

Flame  Temperature 

Simulations  of  burner-stabilized  flames  are  often  conducted  using  experimental  temperature 
profiles  as  inputs  rather  than  solving  the  energy  equation  for  the  temperature  profile.  In  our  case, 
the  energy  equation  was  solved  and  the  temperature  profile  was  a  simulation  output.  The  model 
results  for  flame  temperature  were  compared  to  experimental  data  from  published  CARS 
measurements  [27].  The  maximum  difference  between  the  model  prediction  of  flame 
temperature  was  80  K  and  the  average  difference  was  26  K  at  all  conditions,  from  (j)  =  1  to  1 .4 
and  inlet  velocities  from  7  to  15  cm/s. 

Limit  conditions 


In  practice  a  flat  flame  can  be  stabilized  on  the  burner,  and  if  the  velocity  is  increased  the 
standoff  distance  decreases.  At  some  velocity,  the  standoff  distance  is  minimum,  and  a  further 
small  (~1  cm/s)  increase  of  the  inlet  velocity  causes  the  flame  to  lose  its  flatness.  Further 
increase  always  results  in  the  flame  becoming  less  flat  either  forming  corrugated  or  cellular 
shapes  as  described  in  the  literature  [40]  or  under  very  rich  conditions  forming  more  of  a  cone  or 
dome  shape  (Figure  E.2  is  an  example)  sometimes  with  an  inversion  in  the  center. 


101 


As  described  previously,  the  intention  of  this  study  is  to  understand  the  behavior  of  flat  flames 
and  1-D  numerical  models  only.  In  the  burner-stabilized  flame  simulations,  as  the  velocity  is 
increased  there  is  a  velocity  beyond  which  a  converged  solution  is  not  obtained,  and  converged 
solutions  should  not  be  obtained  under  conditions  where  stable,  flat  flames  are  not  observed 
experimentally.  The  comparison  between  the  experimental  limiting  conditions  and  the  maximum 
velocities  for  which  converged  solutions  to  the  burner  stabilized  flame  simulation  is  shown  in 
Fig.  E.5.  Additionally,  as  described  in  the  literature,  flames  stabilized  on  a  McKenna-type  burner 
should  become  non-flat  when  the  flame  is  nearly  adiabatic  [31,  40].  Therefore  the  experimental 
limiting  conditions  are  also  compared  to  the  calculated  laminar  flame  speeds  in  Fig.  E.5. 
Throughout  the  range  of  equivalence  ratios,  the  match  between  the  experiments  and  burner- 
stabilized  flame  computations  is  very  good.  From  (|>  =  0.6  to  (()  ~  1.2  the  flame  loses  its  flatness 
below  the  laminar  flame  speed  especially  near  stoichiometric.  When  (j)  is  moderately  rich,  the 
flame  loses  its  flatness  at  an  inlet  velocity  near  the  laminar  flame  speed,  but  then  at  lower  inlet 
velocities  in  comparison  to  the  laminar  flame  speed  as  the  equivalence  ratio  approaches  the 
maximum  value  tested  (c[)  ~  1.6). 
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Figure  E.5:  Limit  inlet  velocity  and  laminar  flame  speed  vs.  equivalence  ratio 


These  comparisons  indicate  that  the  computational  model  predicts  flame  standoff  distance  and 
temperature  across  a  broad  range  of  4>.  In  addition,  the  limit  velocities  for  the  experiments  show 
excellent  agreement  with  the  velocities  of  maximum  convergence  for  the  burner-stabilized  flame 
model.  The  laminar  flame  speeds  exceed  the  limit  velocities  determined  experimentally  under 
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lean  and  near  stoichiometric  conditions,  but  not  as  substantially  under  richer  conditions,  which 
are  the  focus  of  this  paper.  These  comparisons  establish  confidence  that  the  models  can  be  used 
to  predict  the  behavior  of  burner-stabilized  flames  especially  at  rich  conditions. 

Effect  of  preheating  on  the  conversion  of  rich  mixtures  to  syngas 

To  provide  a  baseline  for  comparison  to  the  results  for  flames  with  preheated  reactants, 
experiments  and  simulations  were  conducted  with  reactants  at  standard  conditions.  These  data 
are  important  also  to  complete  the  comparison  between  the  models  and  experiments  and  to 
establish  the  similarity  between  the  two  burners  used  experimentally.  Then  data  are  presented  for 
experiments  and  modeling  of  burner-stabilized  flames  with  preheated  reactants  to  show  the  effect 
of  preheat  on  syngas  production. 

Comparison  of  McKenna  and  ceramic  burner  at  baseline  ( unheated )  conditions 

As  was  discussed  previously,  the  McKenna  burner  cannot  withstand  the  preheat  temperatures  of 
interest  in  this  study  so  we  constructed  a  ceramic  burner  for  this  purpose.  In  Fig.  E.6,  the 
concentrations  of  H2,  CO  and  C02  are  shown  for  experiments  with  both  burners  as  well  as 
predictions  from  burner- stabilized  flame  simulations,  freely-propagating  flame  simulations,  and 
equilibrium  calculations.  The  error  bars  shown  in  Fig.  E.6  are  the  maximum  uncertainties  in  the 
presented  experimental  data.  The  average  uncertainty  in  the  concentration  measurements  are  as 
follows:  ±0.29  for  [H2],  ±0.32  for  [CO],  ±0.4  for  [C02]. 


Figure  E.6:  Concentration  vs.  equivalence  ratio  at  inlet  velocity  =  15  ±  0.2  cm/s  (circled 
points  indicate  ceramic  burner  data,  all  computational  data  are  from  burner-stabilized 
flame  model  except  boxed  points  indicating  free  flame  computations  at  ()>  =  1.38) 
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As  Fig.  E.6  shows,  the  experimental  data  from  both  reactors  and  the  models  match  in 
concentration  extremely  well  over  the  tested  range  of  (j>.  The  ceramic  burner  has  a  flashback 
limit,  so  data  are  presented  for  (j)  from  1.25  to  1.35  where  flat  flames  were  stabilized  on  the 
surface  of  the  burner.  At  stoichiometric  conditions,  very  nearly  all  of  the  carbon  in  methane  is 
converted  to  C02  and  little  CO  is  observed.  As  (j)  increases,  the  amount  of  carbon  converted  to 
C02  decreases  as  the  amount  converted  to  CO  increases.  Negligible  amounts  of  small 
hydrocarbons  were  measured  by  the  GC.  Therefore,  as  (j)  is  increased,  more  hydrogen  in  methane 
is  converted  to  H2  at  the  expense  of  H20.  Both  computations  and  experiments  showed  very 
nearly  equilibrium  levels  of  C02,  CO  and  H2.  Lastly,  the  concentration  of  syngas  increases 
monotonically  with  4>  suggesting  that  even  richer  (j)  might  improve  syngas  production  if  stable 
operation  could  be  accomplished. 

Concentrations  for  CO,  C02  and  H2  are  almost  identical  to  the  equilibrium  values  and  are 
consistent  with  the  trend  of  species  concentrations  seen  in  the  burner-stabilized  flame 
experiments  and  models.  In  the  experiments  if  (j)  is  increased,  the  standoff  distance  increases,  but 
the  flame  remains  flat  (burner-stabilized)  until  the  (j)  ~  1.32  on  the  McKenna  burner  and  (|)  ~  1.35 
on  the  ceramic  burner.  When  (j)  is  increased  further,  the  flame  becomes  unstable  because  the 
laminar  flame  speed  of  the  mixture  approaches  V;  under  this  condition  the  flame  can  be 
considered  nearly  adiabatic,  and  the  predictions  of  the  freely-propagating  flame  model  become 
appropriate  for  comparison.  At  300  K,  a  free  flame  of  methane  and  air  has  a  laminar  flame  speed 
of  15  cm/s  when  (j)  ~  1.38. 

We  measured  species  concentrations  while  varying  V  and  holding  (j)  constant  at  1.25  to 
determine  if  the  two  burners  and  the  model  results  were  consistent  across  a  range  of  V.  These 
data  are  shown  in  Fig.  E.7. 
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Figure  E.7:  Concentration  vs.  inlet  velocity  at  equivalence  ratio  =  1.25  ±  0.05  (circled  points 
indicate  ceramic  burner  data,  all  computational  data  are  from  burner-stabilized  flame 
model  except  boxed  points  indicating  free  flame  computations  at  V  =  29  cm/s) 


The  computations  predict  almost  no  change  in  species  concentration  as  a  function  of  V  in  the 
range  of  inlet  velocities  that  were  tested  experimentally.  However,  at  very  low  inlet  velocities  the 
computations  predict  that  more  methane  is  converted  to  C02  at  the  expense  of  CO.  The 
experiments  show  almost  no  change  in  concentration  for  both  the  McKenna  burner  and  the 
ceramic  burner.  The  equivalence  ratio  strongly  affects  species  concentrations,  as  shown  in  Fig. 
E.6,  and  there  is  uncertainty  in  (j)  for  all  experimental  data  points  when  V  is  varied  and  (j>  is  held 
constant.  The  uncertainties  in  species  concentration  due  to  uncertainty  in  (j)  around  (j)  =  1.25  are 
estimated  as  follows:  ±0.8  for  [H2],  ±0.9  for  [CO]  and  ±0.8  [C02].  Because  of  these 
uncertainties,  which  are  shown  in  Fig.  E.7  because  they  are  larger  than  the  raw  measurement 
uncertainty,  no  variation  in  the  experimental  data  as  a  function  of  inlet  velocity  is  implied  in  Fig. 
E.7. 


The  major  conclusions  from  the  previously  presented  data  are  that  the  ceramic  burner  and  the 
McKenna  burner  produce  nearly  identical  species  concentrations  across  a  range  of  (j)  and  V.  In 
addition,  the  burner-stabilized  flame  model  and  the  free  flame  model  were  very  accurate  in 
predictions  of  standoff  distance,  limit  conditions  and  [CO],  [C02]  and  [H2]  at  all  (j)  and  V  tested. 
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Results  from  experiments  and  modeling  with  preheated  reactants 


Using  the  ceramic  burner,  we  performed  experiments  with  preheated  reactants.  Fig.  E.8  shows 
concentration  as  a  function  of  (j)  with  V  held  constant  at  25  cm/s  and  the  reactants  preheated  to 
617  K.  Under  unheated  conditions  on  the  ceramic  burner  for  V  =  25  cm/s  the  limit  (j)  is  1.26. 


Figure  E.8:  Concentration  vs.  equivalence  ratio  at  inlet  velocity  =  25  +  0.05  cm/s  and 
reactant  temperature  =  617  ±  10  K  (all  computational  data  are  from  burner-stabilized 
flame  model  except  boxed  points  indicating  free  flame  computations  at  <)>  =  1.9,  uncertainty 

in  temperature  is  ±  10  K) 

The  burner-stabilized  flame  model  matches  the  experimental  data  very  closely.  Under  this 
preheat  condition,  a  stable  flame  was  obtained  up  to  (j)  =  1.75,  which  is  greater  than  the 
conventional  flammability  limit  of  1.67  and  significantly  greater  than  the  operating  limit 
observed  when  reactants  were  unheated  ((j)  =  1.26).  The  maximum  [H2]  and  [CO]  achievable 
under  unheated  conditions  at  V  =  15  cm/s  were  6.5%  and  8.3%,  respectively,  both  occurring  at  (}) 
=  1.35  as  shown  in  Fig.  E.5.  With  reactants  preheated,  the  maximum  [H2]  was  12.7%  and  the 
maximum  [CO]  was  12%.  These  data  show  that  preheating  significantly  extends  the  operating 
range  of  burner-stabilized  flames  such  that  the  stable  burning  of  mixtures  that  produce  high 
syngas  yields  can  be  achieved.  The  data  also  show  that  the  burner- stabilized  flame  model 
produces  very  good  predictions  of  species  at  (j)  even  beyond  the  rich  flammability  limit. 

It  has  been  shown  that  the  limiting  conditions  for  the  experiments,  the  burner  stabilized  flame 
model  and  the  laminar  flame  speed/equivalence  ratio  pairs  were  nearly  coincident  for  flames 
with  unheated  reactants.  As  shown  in  Fig.  E.8  with  reactants  preheated  to  617  K  the  limiting 
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equivalence  ratio  (4>  =  1.75)  for  a  converged  solution  with  the  burner- stabilized  flame  model  is 
nearly  identical  to  the  experimental  value  as  well.  However,  the  mixture  equivalence  ratio  (cj)  = 
1.9)  having  a  flame  speed  of  25  cm/s  was  significantly  greater  than  the  experimental  and  model 
limit  for  the  burner-stabilized  flame. 

A  difference  between  these  data  and  the  data  with  unheated  reactants  is  that  the  experimental 
results  and  the  computational  results  do  not  match  the  equilibrium  values  at  the  highest  values  of 
(j).  Up  to  the  maximum  achievable  cj)  for  the  burner-stabilized  flame,  the  experiments  and  the 
model  concentrations  are  very  close  to  equilibrium  for  [C02],  however  [CO]  begins  to  diverge 
from  equilibrium  slightly  and  [H2]  begins  to  diverge  significantly  when  cj)  is  greater  than  about 
1.4.  This  divergence  is  somewhat  unexpected  since  preheating  should  increase  the  rate  at  which 
the  mixture  reaches  equilibrium  relative  to  unheated  conditions.  However,  under  the  preheated 
conditions,  the  cj)  that  were  accessible  were  all  relatively  high  (cj)  >  1.4);  at  high  cj)  the  heat  release 
is  low,  thus  slowing  the  kinetic  processes,  and  the  relatively  slow  water-gas  shift  reaction 
becomes  important.  Even  though  preheating  increases  the  range  of  accessible  cj)  by  increasing 
reaction  rates,  the  reaction  rate  increase  is  not  great  enough  to  allow  full  equilibration  of  the 
mixture.  At  very  rich  equivalence  ratios,  the  relatively  slow  steam  reforming  reaction  becomes 
important  Lastly,  the  freely-propagating  flame  model  shows  even  more  substantial  divergence 
from  equilibrium  for  both  [CO]  and  [H2]  while  [C02]  remains  very  close  to  the  equilibrium 
value.  The  concentrations  for  the  freely-propagating  flame  are  consistent  with  the  trend  in 
concentration  for  the  burner-stabilized  flame.  However,  the  difference  between  the  limit  cj)  for  the 
burner-stabilized  flame  model  and  experiment  and  the  cj)  having  a  flame  speed  of  25  cm/s  casts 
doubt  on  the  accuracy  of  the  freely-propagating  flame  model  under  this  preheated  condition. 

The  effect  of  reactant  temperature  was  investigated  by  preheating  reactants  while  the  equivalence 
ratio  was  held  at  1.4  and  the  inlet  velocity  was  held  at  24  cm/s  (Fig.  E.9).  The  minimum 
temperature  at  which  a  stable  flame  was  observed  for  these  conditions  was  365  K,  and  the 
maximum  temperature  to  which  reactants  were  preheated  in  experiments  was  630  K. 
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Reactant  Temperature  (K) 

Figure  E.9:  Concentration  vs.  reactant  temperature  at  inlet  velocity  =  24  +  2  cm/s  and 
equivalence  ratio  =  1.4  +  0.05  (all  computational  data  are  from  burner-stabilized  flame 
model  except  boxed  points  indicating  free  flame  computations  at  temperature  =  385  K) 

In  the  range  of  temperature  where  there  was  overlap  between  the  experiments  and  computations, 
the  computations  predicted  [H2]  and  [CO]  very  well  while  under-predicting  [C02]  slightly.  As 
was  shown  in  all  of  the  previously  presented  data,  the  experimental  results,  burner-stabilized 
flame  model  and  freely-propagating  flame  model  had  species  concentrations  nearly  identical  to 
the  equilibrium  values  at  unheated  conditions.  As  the  preheat  temperature  increases  beyond  700 
K,  the  difference  between  equilibrium  and  the  burner-stabilized  flame  model  increases  for  both 
[CO]  and  [C02],  however  the  difference  between  the  model  and  equilibrium  in  terms  of  [H2] 
does  not  increase  significantly  with  temperature.  Lastly,  the  overall  yield  of  syngas  ([H2]  and 
[CO])  remains  constant  with  preheat  temperature  indicating  that  preheating  is  important  for 
stabilizing  a  reaction,  but  not  for  improving  syngas  yield.  The  concentration  uncertainties  shown 
in  this  plot  are  due  to  uncertainty  in  equivalence  ratio  because  these  uncertainties  are  larger  than 
the  raw  measurement  uncertainties. 

The  minimum  temperature  at  which  a  solution  converged  with  the  burner-stabilized  flame  model 
was  within  90  K  of  the  limit  temperature,  which  is  the  minimum  temperature  at  which  a  stable, 
flat  flame  was  observed  experimentally.  The  preheat  temperature  at  which  the  freely-propagating 
flame  model  reported  a  flame  speed  of  25  cm/s  was  within  20  K  of  this  limit  temperature.  The 
coincidence  of  this  temperature  amongst  the  models  is  another  indication  that  the  models  are 
accurate  in  their  prediction  of  flame  characteristics  at  low  reactant  temperature  conditions. 
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E.5  Concluding  Remarks 


An  experimental  and  modeling  study  of  burner- stabilized  flames  was  conducted  to  understand 
the  effect  of  preheat  on  syngas  production.  First  experiments  were  conducted  with  unheated 
reactants  so  that  confidence  in  the  models  could  be  gained  for  their  prediction  of  combustion 
characteristics,  especially  species  yields,  at  high  equivalence  ratios.  The  effect  of  preheat  was 
investigated  by  preheating  reactants  up  to  630  K  while  holding  the  inlet  velocity  and  equivalence 
ratio  constant. 

The  results  showed  that  the  burner  stabilized  flame  model  accurately  predicts  standoff  distance, 
limit  conditions  and  species  yields  for  experiments  with  flat  flame  burners.  Experiments  with 
preheated  reactants  demonstrated  that  the  operating  limits  for  a  flat  flame  burner  increase 
substantially  with  increased  reactant  temperature,  allowing  the  stable  operation  of  the  burner  at 
richer  equivalence  ratios — where  syngas  yields  are  greater — than  those  achievable  with  unheated 
reactants.  The  specific  effect  of  inlet  temperature  was  investigated,  and  it  was  found  that  inlet 
temperature  did  not  affect  the  yield  of  syngas  under  the  tested  range  of  temperatures.  Burner- 
stabilized  flames  with  preheated  reactants  showed  deviations  from  equilibrium  results  in  terms  of 
species  yields,  in  contrast  with  unheated  conditions  where  equilibrium  predicts  species  yields  for 
burner-stabilized  flames  very  well.  The  burner-stabilized  flame  model  was  able  to  predict  the 
deviation  from  equilibrium  when  reactants  were  preheated  to  -630  K.  At  high  preheat 
temperatures,  the  limiting  operating  conditions  of  the  experiment  were  well-predicted  by  the 
burner-stabilized  flame  model  but  not  by  the  laminar  flame  speed. 

These  results  have  important  implications  with  regards  to  heat  recirculating  reactors.  First,  the 
models,  and  therefore  the  kinetics  mechanism,  performed  very  well  up  to  richest  equivalence 
ratio  examined  of  1.75,  which  exceeds  the  conventional  rich  flammability  limit.  The  only 
significant  difference  between  the  models  and  the  experiments  was  the  difference  between  the 
limit  equivalence  ratio  and  the  equivalence  ratio  where  the  laminar  flame  speed  was  25  cm/s 
when  the  reactants  were  preheated  to  617  K.  This  means  that  those  working  on  reactor  models 
for  rich  combustion  should  have  confidence  in  GRI3.0  for  predicting  syngas  production  even 
beyond  the  rich  flammability  limit.  In  heat-recirculating  reactors,  preheat  is  accomplished  by 
transport  enhancement,  which  occurs  in  proximity  to  the  reaction  zone,  as  part  of  the  reactor 
design.  In  the  experiments  described  here,  the  preheat  occurred  far  upstream  of  the  reaction  zone 
and  was  input  from  an  external  source.  This  experimental  design  provided  a  means  to  study  the 
effect  of  preheat  temperature  independent  of  enhanced  transport  mechanisms  and  interaction 
with  the  reaction  zone.  The  results  show  that  preheat  is  necessary  to  attain  stable  operation  at 
extreme  equivalence  ratios.  This  was  the  expected  result  since  conditions  (such  as  downstream 
propagation  of  the  reaction  zone  in  a  filtration  reactor)  that  enhance  preheat  have  been  observed 
in  heat-recirculating  reactors  when  mixtures  expected  to  be  unstable  are  burned.  It  was  also 
found  that  further  preheat  does  not  enhance  syngas  production.  This  result  aids  in  the  design  of 
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heat  recirculating  reactors  because  it  indicates  that  maximal  heat  recirculation  may  not  be 

necessary  for  conversion  of  fuel  to  syngas. 

E.6  References 

1.  Weinberg,  F.J.,  Combustion  Temperatures:  The  Future?  Nature,  1971.  233:  p.  239. 

2.  Jones,  A.R.,  S.A.  Lloyd,  and  F.J.  Weinberg,  Combustion  in  Heat  Exchangers.  Proc.  R. 
Soc.  Lond.  A,  1978.  360:  p.  97-115. 

3.  Brenner,  G.,  K.  Pickenacker,  O.  Pickenacker,  D.  Trimis,  K.  Wawrzinek,  and  T.  Weber, 
Numerical  and  Experimented  Investigation  of  Matrix-Stabilized  Methane/Air  Combustion 
in  Porous  Inert  Media.  Combust.  Flame,  2000.  123:  p.  201-213. 

4.  Sathe,  S.B.,  M.R.  Kulkarni,  R.E.  Peck,  and  T.W.  Tong,  An  Experimental  and  Theoretical 
Study  of  Porous  Radiant  Burner  Performance.  Proc.  Combust.  Inst.,  1990.  23:  p.  1011- 
1018. 

5.  Wood,  S.  and  A.T.  Harris,  Porous  burners  for  lean-burn  applications.  Progress  in  Energy 
and  Combustion  Science,  2008.  34(5):  p.  667-684. 

6.  Kim,  N.I.,  S.  Kato,  T.  Kataoka,  T.  Yokomori,  S.  Maruyama,  T.  Fujimori,  and  K.  Maruta, 
Flame  Stabilization  and  Emission  of  Smcdl  Swiss -Roll  Combustors  as  Heaters.  Combust. 
Flame,  2005.  141:  p.  229-240. 

7.  Al-Hamamre,  Z.,  S.  VoB,  and  D.  Trimis,  Hydrogen  production  by  thermal  partied 
oxidation  of  hydrocarbon  fuels  in  porous  media  based  reformer.  International  Journal  of 
Hydrogen  Energy,  2009.  34(2):  p.  827-832. 

8.  Kennedy,  L.A.,  J.P.  Bingue,  A.V.  Saveliev,  A.A.  Fridman,  and  S.I.  Foutko,  Chemical 
structures  of  methane-air  filtration  combustion  waves  for  fuel-lean  and  fuel-rich 
conditions.  Proceedings  of  the  Combustion  Institute,  2000.  28(1):  p.  1431-1438. 

9.  P.  Bingue,  J.,  A.  V.  Saveliev,  A.  A.  Fridman,  and  F.A.  Kennedy,  Hydrogen  production  in 
ultra-rich  filtration  combustion  of  methane  and  hydrogen  sulfide.  International  Journal  of 
Hydrogen  Energy,  2002.  27(6):  p.  643-649. 

10.  Pedersen-Mjaanes,  H.,  F.  Chan,  and  E.  Mastorakos,  Hydrogen  production  from  rich 
combustion  in  porous  media.  International  Journal  of  Hydrogen  Energy,  2005.  30(6):  p. 
579-592. 

11.  Pastore,  A.  and  E.  Mastorakos,  Rich  n-heptane  and  eliesel  combustion  in  porous  media. 
Experimental  Thermal  and  Fluid  Science,  2010.  34(3):  p.  359-365. 

12.  Koch,  D.L.,  R.G.  Cox,  H.  Brenner,  and  J.F.  Brady,  The  Effect  of  Order  on  Dispersion  in 
Porous  Media.  J.  Fluid  Mech,  1989.  200:  p.  173-188. 

13.  Hackert,  C.L.,  J.L.  Ellzey,  O.A.  Ezekoye,  and  M.J.  Hall,  Transverse  Dispersion  at  High 
Peclet  Numbers  in  Short  Porous  Media.  Exp.  Fluids,  1996.  21:  p.  286-290. 

14.  Dhamrat,  R.S.  and  J.L.  Ellzey,  Numerical  and  Experimental  Study  of  the  Conversion  of 
Methane  to  Hydrogen  in  a  Porous  Media  Reactor.  Combustion  and  Flame,  2006.  144:  p. 
698-709. 


110 


15.  Schoegl,  I.  and  J.L.  Ellzey,  A  mesoscale  fuel  reformer  to  produce  syngas  in  portable 
power  systems.  Proceedings  of  the  Combustion  Institute,  2009.  32(2):  p.  3223-3230. 

16.  Drayton,  M.K.,  A.V.  Saveliev,  L.A.  Kennedy,  A. A.  Fridman,  and  Y.  Li,  Syngas 
Production  Using  Superadiabatic  Combustion  of  Ultra-Rich  Methane -Air  Mixtures.  Proc. 
Combust.  Inst.,  1998.  27:  p.  1361-1367. 

17.  Dixon,  M.J.,  I.  Schoegl,  C.B.  Hull,  and  J.L.  Ellzey,  Experimental  and  numerical 
conversion  of  liquid  heptane  to  syngas  through  combustion  in  porous  media.  Combustion 
and  Flame,  2008.  154(1-2):  p.  217-231. 

18.  Belmont,  E.,  S.M.  Solomon,  and  J.L.  Ellzey,  Syngas  Production  from  Heptane  in  a  N on- 
catalytic  Counterflow  Reactor.  Combustion  and  Flame,  2012.  In  Press. 

19.  Barra,  A.J.  and  J.L.  Ellzey,  Heat  Recirculation  and  Heat  Transfer  in  Porous  Burners. 
Combust.  Flame,  2004.  137:  p.  230-241. 

20.  Ronney,  P.D.,  Analysis  of  Non-Adiabatic  Heat-Recirculating  Combustors.  Combust. 
Flame,  2003.  135:  p.  421-439. 

21.  Chen,  C.-H.  and  P.D.  Ronney,  Three-dimensional  effects  in  counterflow  heat¬ 
recirculating  combustors.  Proceedings  of  the  Combustion  Institute,  2010.  33(2):  p.  3285- 
3291. 

22.  Mishra,  S.C.,  M.  Steven,  S.  Nemoda,  P.  Talukdar,  D.  Trimis,  and  F.  Durst,  Heat  transfer 
analysis  of  a  two-dimensional  rectangular  porous  radiant  burner.  International 
Communications  in  Heat  and  Mass  Transfer,  2006.  33(4):  p.  467-474. 

23.  Khose-Hoinghaus,  K.  and  J.B.  Jeffries,  eds.  Applied  Combustion  Diagnostics. 
Combustion:  An  International  Series.  2002,  Taylor  &  Francis:  New  York. 

24.  Hansen,  N.,  T.A.  Cool,  P.R.  Westmoreland,  and  K.  Kohse-Hoinghaus,  Recent 
contributions  of  flame -sampling  molecular-beam  mass  spectrometry  to  a  fundamental 
understanding  of  combustion  chemistry.  Progress  in  Energy  and  Combustion  Science, 
2009.  35(2):  p.  168-191. 

25.  Marinov,  N.M.,  W.J.  Pitz,  C.K.  Westbrook,  M.J.  Castaldi,  and  S.M.  Senkan,  Modeling  of 
Aromatic  and  Polycyclic  Aromatic  Hydrocarbon  Formation  in  Premixed  Methane  and 
Ethane  Flames.  Combustion  Science  and  Technology,  1996.  116-117(1-6):  p.  211-287. 

26.  Schoenung,  M.  and  R.K.  Hanson,  CO  and  Temperature  Measurements  in  a  Flat  Flame 
by  Faser  Ab  orption  Spectroscopy  and  Probe  Techniques.  Combustion  Science  and 
Technology,  1980.  24(5-6):  p.  227-237. 

27.  Weigand,  P.,  R.  Luckerath,  and  W.  Meier,  Documentation  of  Flat  Premixed  Faminar 
CH4/Air  Standard  Flames:  Temperatures  and  Species  Concentrations ,  D.I.o.C.T. 
(http://www.dlr.de/vt/PortalData/29/Resources/dokumente/ch4  air  flames.pdf).  Editor. 
2003:  Stuttgart. 

28.  Gersen,  S.,  A.V.  Mokhov,  and  H.B.  Levinsky,  Extractive  probe/TDLAS  measurements  of 
acetylene  in  atmospheric-pressure  fuel-rich  premixed  methane/air  flames.  Combustion 
and  Flame,  2005.  143(3):  p.  333-336. 


Ill 


29.  Kumar,  K.  and  C.-J.  Sung,  Laminar  flame  speeds  and  extinction  limits  of  preheated  n- 
decane/02/N2  and  n-dodecane/02/N2  mixtures.  Combustion  and  Flame,  2007.  151:  p. 
209-224. 

30.  Liu,  F.S.  and  O.L.  Gulder,  Effects  of  pressure  and  preheat  on  super-adiabatic  flame 
temperatures  in  rich  premixed  methane/air  flames.  Combustion  Science  and  Technology, 
2008.  180(3):  p.  437-452. 

31.  Bosschaart,  K.J.  and  L.P.H.  de  Goey,  The  laminar  burning  velocity  of  flames  propagating 
in  mixtures  of  hydrocarbons  and  air  measured  with  the  heat  flux  method.  Combustion 
and  Flame,  2004.  136(3):  p.  261-269. 

32.  Gibson,  J.,  M.  Ayoobi,  and  I.  Schoegl,  Behavior  of  Preheated  Premixed  Flames  at  Rich 
Conditions.  Proc.  Combust.  Inst.,  2012.  34(In  Press). 

33.  Sepman,  A.,  A.V.  Mokhov,  and  H.B.  Levinsky,  A  laser-induced  fluorescence  and 
coherent  anti-Stokes  Raman  scattering  study  of  NO  formation  in  preheated,  laminar,  rich 
premixed,  methane/air  flames.  Proceedings  of  the  Combustion  Institute,  2002.  29:  p. 
2187-2194. 

34.  Hidaka,  Y.,  K.  Sato,  Y.  Henmi,  H.  Tanaka,  and  K.  Inami,  Shock-tube  and  modeling  study 
of  methane  pyrolysis  and  oxidation.  Combustion  and  Flame,  1999.  118(3):  p.  340-358. 

35.  Goodwin,  D.G.  (2012)  Cantera  -  An  object-oriented  software  toolkit  for  chemical 
kinetics,  thermodynamics,  and  transport  processes,  http://code.google.eom/p/cantera/ 

36.  Smith,  G.P.,  D.M.  Golden,  M.  Frenklach,  N.W.  Moriarty,  B.  Eiteneer,  M.  Goldenberg, 
C.T.  Bowman,  R.K.  Hanson,  S.  Song,  J.W.C.  Gardiner,  V.V.  Lissianski,  and  Z.  Qin, 
GRI-Mech  3.0  URL  http://www.me.berkeley.edu/gri-mech/version30/.  1999. 

37.  Law,  C.K.,  Combustion  Physics.  2006:  Cambridge  University  Press. 

38.  Sutton,  G.,  A.  Levick,  G.  Edwards,  and  D.  Greenhalgh,  A  combustion  temperature  and 
species  standard  for  the  calibration  of  laser  diagnostic  techniques.  Combustion  and 
Flame,  2006.  147:  p.  39-48. 

39.  Mokhov,  A.V.  and  H.B.  Levinsky,  A  L1F  and  cars  investigation  of  upstream  heat  loss 
and  flue-gas  recirculation  as  NOx  control  strategies  for  laminar,  premixed  natural- 
gas/air  flames.  Proceedings  of  the  Combustion  Institute,  2000.  28(2):  p.  2467-2474. 

40.  Ferguson,  C.R.  and  J.C.  Keck,  Stand-off  distances  on  a  flat  flame  burner.  Combustion 
and  Flame,  1979.  34:  p.  85-98. 

41.  Kurdyumov,  V.N.  and  M.  Matalon,  The  porous-plug  burner:  Flame  stabilization,  onset 
of  oscillation,  and  restabilization.  Combustion  and  Flame,  2008.  153:  p.  105-118. 

42.  Williams,  B.A.  and  J.W.  Fleming,  Comparison  of  species  profiles  between  02  and  N02 
oxidizers  in  premixed  methane  flames.  Combustion  and  Flame,  1995.  100(4):  p.  571-590. 

43.  Evertsen,  R.,  J.A.  Van  Oijen,  R.T.E.  Hermanns,  L.P.H.  De  Goey,  and  J.J.  Ter  Meulen, 
Measurements  of  absolute  concentrations  of  CH  in  a  premixed  atmospheric  flat  flame  by 
cavity  ring-down  spectroscopy.  Combustion  and  Flame,  2003.  132:  p.  34-42. 


112 


